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FLOOD- FREQUENCY RELATIONSHIPS IN THE PACIFIC NORTHWEST 


By George L. Bodhaine,! M, ASCE 


SYNOPSIS 


This paper describes, in general terms, a method of determining the mag- 
nitude and frequency of expected floods that is applicable to any area in the 
Pacific Northwest of the United States. The development of a composite fre- 
quency curve applicable to any site in a large homogeneous region is ex- 
amined, and the effect of mean elevation and size of drainage area on the 
shape of the frequency curves in designated parts of the Pacific Northwest 
.§ demonstrated. A formula for determining the mean annual flood for any 
site on the west side of the Cascade Range in Washington is shown, and the 
rainage-basin characteristics used in its derivation are listed, A flood- 
frequency curve can be drawn for any given site in the Pacific Northwest 
y combining the results obtained from the mean annual flood formula and 
he composite frequency curve that are both applicable to that site. Com- 
ments about limitations on the use of results obtained by this method are 


INTRODUCTION 


The purpose of this paper is to present a method of correlating stream- 
‘low data and basin characteristics to obtain information relative to the 


method reflects the latest developments based on a continuing study of the 
ubject by engineers of the United States Geological Survey, 


 Note.—Discussion open until April 1, 1961. To extend the closing date one month, a 
ritten request must be filed with the Executive Secretary, ASCE. This paper is part 
the copyrighted Journal of the Hydraulics Division, Proceedings of the American So- 
iety of Civil Engineers, Vol. 86, No. HY 9, November, 1960. 
1 Flood Specialist, U. S. Geological Survey, Tacoma, Washington. 
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Interior (USGS). It is expected that the method will be revised from time t 
time as new ideas are conceived and tested. 

The discharge records collected by the USGS and others are the data 
on which flood-frequency studies are based. The records from about 800 
gaging stations in the Pacific Northwest were of sufficient length to furnish 
data for deriving the flood-frequency curves. A curve based on a gaging- 
station record, however, applies only to that one particular site, whereas, 
generally, the information is desired at an ungaged point, frequently on an 
ungaged stream. The two parts to the problem are (1) computation of flood 
frequency at a point, and (2) adaption of the point data to apply over abasin 
or a region. 

In general, the needs of most engineering studies are met by information 
concerning floods with average recurrence intervals of 10 yr, 25 yr, 50 
yr and 100 yrs. For example, some states’ Highway Departments design 
their primary highway drainage structures to withstand floods that are 
expected to occur, on the average, once in 50 yrs whereas their secondary 
highway structures are designed to withstand floods that are expected to 
occur, on the average, once in 25 yr. 

The Pacific Northwest is defined as the Columbia River basin and the 
western slope of the Cascade Range in Washington and Oregon. This area is 
designated as Parts 12, 13, and 14 by the USGS. 


FLOOD FREQUENCY AT A GAGING STATION 


The first step in a flood compilation is to select the gaging stations to be 
included. There should be at least 10 yr of record at each station, although 
records as short as 5 yr are used when no others are available. Adjust- 
ments should be made for storage, or other artificial factors, that, tend to 
modify flood discharges significantly, or the record should not be used, 

Flood Series.—There are two kinds of flood series to be considered in 
plotting frequency curves, the annual-flood series and the partial-duration 
series. The two methods give essentially the same results for recurrence 
intervals greater than 10 yr. A relationship between the two series has been 
derived by Langbein2, of the USGS, so that, when necessary, the annual-flood 
curve can be converted to the partial-duration series, For example, a flood 
having a recurrence interval of 2.0 yr in the annual-flood series would have 
a recurrence interval of 1.45 yr in the partial-duration series, At a recurr- 
ence interval of 10.5 yr the comparative value would be 10.0 yr and at 50.5 
yr, 50.0 yr. The annual-flood method is used in these studies because of its 
simplicity. 4 

Plotting Procedure.—Recurrence intervals are computed from the formula 


Ties (N's 10/Mi rears sel dtc oeogsens () 


in which T is recurrence interval in years, N is the number of years of rec- 
ord, and M is the relative magnitude of the event beginning with the highest 
as 1. : 

Annual floods are plotted on a special form developed for analysis of flood 
frequencies by the Gumbel8 method. Discharge is plotted as ordinate and time 


2 « annual Floods and the Partial-Duration Flood Series,” by W. B, Langbein, Ame: 
ican Geophysical Union Transactions, Vol. 30, 1949, pp. 879-881, ay 
3 “Floods Estimated by the Probability Method,” by E. J. Gumbel, Eng > News 
Record, June 14, 1945, 
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Ss abcissa. Both linear and logarithmic ordinate scales are used, depending 
nm which plot more nearly approximates a straight line. An example of each 
ind of plot is shown in Fig. 1. The same data are plotted here to show the 
haracteristic difference in the shape of the final curves, 

Fitting Frequency Graphs.— After the annuai floods are plotted on a fre- 
uency diagram, a curve is fitted to the data. Considering the fact that most 
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7 FIG. 1.—MAGNITUDE AND FREQUENCY OF ANNUAL FLOODS, 


freamflow records are less than 35 yr inlength, and analytically fitted func- 
on has little value. Therefore, graphical treatment only is applied, and a 
imple curve of visual best-fit is drawn. 


REGIONAL DISTRIBUTION 


— 


-Flood-frequency compilations, prepared in accordance with procedures 
escribed, result in a curve of flood magnitude versus recurrence intervals 
rr each station studied. The curve is not applicable directly to ungaged 
In most instances, the flood records are short and time-sampling 
rrors may be large. That is, the flows sampled by the period of record 
na not be representative of conditions for long periods. Also, the records 
re, generally, for different periods of time. By combining the records for 
1 stations in a drainage basin ora larger region, time-sampling and areal- 
1 pling errors are reduced, and the resulting composite curves are ap- 
ficable to ungaged areas. However, these curves simply show floods of 
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various recurrence intervals in terms of mean annual flood, so the latte 
must be computed before the composite curve can be applied to an ungage 
area. That can be accomplished by relating the mean annual flood to meas 
urable characteristics of the drainage basin. 

Combining Records.—If the flood data for gaging stations at which onl 
records of short length have been collected could be added together to mak 
long records, than five records each 20 yr long would provide a 100-yr re 
cord, Unfortunately, this cannot be done because the data are not independent 
In a homogeneous region, many of the gaging stations record the same floo 
event so we simply have a20-yr record at each of the five stations. Howevel 
since the stations are measuring the same events, they may be averaged to 
gether to provide a better measure of the frequency characteristics. 

The stations are grouped into homogeneous regions, generally on the basi 
of the slopes of the individual frequency curves. The slope usedis the ratio o 
the 10-yr to 2.33-yr flood. According to Gumbel’s theory of extreme values 
the mean of a large number of annual floods should equal the flood with re 
currence interval of 2.33 yr. The slopes of the upper ends of the curves ar 
considered also, especially if the curves break sharply to the left or right 

Time Base.—Experience has demonstrated that floods for one station can 
not be compared with floods for another station with a different length o 
period of record. In studies where mean annual floods are used as a correla 
tion factor, the records must be comparable so they must be adjusted to th 
same period of time. 

The base time period may be chosen equal to any period desired for whic 
records are available, but, ingeneral, itis best to use the period of the longes 
records. A37-yr base, 1921-57, was used for Idaho, and a 46-yr base, 1912-5" 
was used for Oregon, Washington, and western Montana. 

Mean Annual Flood.—The mean used for each frequency distribution is th 
graphical mean determined by the intersection of the visually best fitting fre 
quency line with the vertical line corresponding tothe 2.33-yr recurrence 
interval, The arithmetic mean for relatively short periods of record can b 
greatly influenced by the chance inclusion of a major flood. The graphica 
mean is preferred because it is more stable, with greater weight given to th 
medium floods than to the extreme floods. 

Composite Frequency Curve.—All homogeneous stations are grouped fo 
the purpose of computing an average frequency curve applicable to that regior 
For each station, the flood discharges for several recurrence intervals ar 
determined from the individual frequency curves. These discharges are divide 
by the mean annual flood inorder toplace them on a dimensionless basis. Th 
medium value of the ratios is thendeterminedfor each of the selected recur 
rence intervals. Each median flood ratio is plotted to its corresponding recur 
rence interval on a frequency chart, and a frequency curve is drawn. Thi 
curve, showing discharge in ratio to the mean annual flood versus recurrence 
interval in years, is based on all signifieant discharge records available an 
may be considered as representing the most likely flood-frequency values fo 
all places in the region. (See Fig. 2). ; 

Adjustments to Composite Curves.— Recent studies of streams in Ida 
have shown that, in some areas, the shape of the curve changes with elevatic 
Also, some investigators have found that the size of the drainage area has” 
effect on the shape of the frequency curve, and thatthe effect increases \ 
increase in recurrence interval. 
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Elevation Adjustments.—In certain areas, such as the Clearwater River 
basin in Idaho, there is a definite relationship between flood discharge and 
elevation. There, the slope of the frequency curve decreases as the elevation 
increases, For example, a 50-yr flood for a stream, with an average basin 
elevation of 2,000 ft is 36% larger than for one of the same size area with an 
average basin elevation of 4,000 ft, andis about twice as large as for a stream 
with this area at an elevation of 6,000 ft. Fig. 3 is a family of curves showing 
the variation of flood ratio with recurrence interval and elevation, 

Area Adjustments.—The use of the dimensionless flood-frequency curve 
assumes that the size of the drainage area has no effect on the slope of the 
frequency curve, However, an effect has been found to exist in several areas 
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FIG, 4.—AREA CORRECTION FACTORS TO REGIONAL CURVE. F 


in the Pacific Northwest. For example, in the Spokane River—Pend Oreille 
River region, the frequency curve is not affected by variation in size of drain- 
age areas up to about 200 sq miles, but for areas greater than 200 sq miles 
the ratios of the larger floods to the mean annual flood become smaller as the 
area increases, The effect is not as pronounced in the Spokane River—Pend 
Oreille River area as in some other regions, Even so, the 50-yr ratio of 
2.05 determined for an area of 200 sq miles is reduced to 1.62 for 20,000 sq 
miles, Fig, 4 is a family of curves showing the amount of correction with 
respect to both drainage area and recurrence interval. The values taken 
from these curves are subtracted from the flood ratios of the composite 


cage to obtain flood ratios that are applicable to the size area being con- 
sidered. ; 


————— oe 
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Derivation of Formula for Mean Annual Flood.— A plot of drainage area 
against mean annual flood in any region usually results in a graph in which 
the standard error of the plotted points about the regression line is large, 
sometimes 100% or more. Sucha large standard error indicates that drainage- 
yasin characteristics other thandrainage area may effect the size of the mean 
annual flood. Sometimes, it is possible to obtain a satisfactory correlation 
xy dividing the region into hydrologic zones in which the runoff per square 
mile is similar, If the standard error of the plotted points for a hydrologic 
zone still seems excessive, other drainage basin characteristics must be 
considered. 

In this study, the standard error of drainage area plotted against mean 
annual flood was excessive, about 90%, so other drainage-basin factors 
were considered. Mean elevation, mean annual runoff, area of lages and 
yonds, precipitation, and channel slope were some of the additional factors 
hat were found to improve the correlations. Graphical-linear-multiple cor- 
relations were made to determine those parameters that were most signif- 
cant before entering into mathematical correlations. 

In the area designated as “the west side of the Cascade Range in Wash- 
ngton” drainage area, mean annual runoff and area of lakes and ponds were 
‘ound to affect the mean annual flood most significantly and were used as 
dependent variables in a correlation with mean annual flood. Residual de- 
yiation of about equal magnitude remaining in these correlations were found 
© group areally, so they were assumed to represent the effect of some un- 
lefined basin characteristic such as geology. A map was prepared on which 
“esidual deviations were shown. The map was divided into areas by drawing 
ines enveloping residuals of about equal size. The median residual for each 
rea was used to represent a geographical factor for that area. 

: The standard error of the multiple correlation, before introducing geo- 
rraphical factors, was about 46%, and the coefficient of correlation was 
about 0.97. After assigning values for geographical location and reducing 
Ihe degrees of freedom by one for each zone used, the standard error was 
“educed to about 24%, and the coefficient of correlation increased to 0.99. 
\ standard error of 24% indicates that about two-thirds of the time a mean 
innual flood computed by the method described should be within 24% of the 
‘orrect value. 

Formulas for other areas in the Pacific Northwest, such as western Ore- 
fon and the Snake River basin, were developed in a similar manner, 

| Formula and Nomograph.—The formula derived for western Washington 
Ss 


Qioess ae Rtldp--MG 0... ea) 


‘which Q refers to the meanannual flood in cubic feet per second, A denotes 
e drainage area in square miles, R is the mean annual runoff in inches for 
eriod 1930-57, L refers to the area of lakes and ponds in percent of drainage 


rea, and G is the geographical factor. 
Eq. 2 may be solved by use of logarithms or by slide rule, but it is also 


es igned for the above formula, one for drainage areas of 0.1 to 10 sq miles 
ud another for areas of 10 to 1,000 sq miles. The nomograph for small 


as is shown in Fig. 5. 
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Derivation of Frequency Curve for Any Site.—The procedure for deriving 
a flood-frequency curve for any site in a region is as follows: (1) Obtain the 
drainage area, mean annual runoff, area of lakes and ponds, and other char- 
acteristics, as necessary, from appropriate maps (Mean annual runoff is ob- 
tained from a map showing isopleths based on gaging-station records.); (2) use 
the nomograph or formula to determine the mean annual flood; (3) select the 
proper flood-ratio curve and determine ratios to mean annual flood for se- 
lected recurrence intervals up to 50 yr; and (4) multiply these ratios by the 
mean annual flood computed in step 2 and plot against recurrence intervals 
to define the frequency curve. 


LIMITATIONS 


The methods described may be used to predict the magnitude of floods of 
givenreturnintervals on streams in the Pacific Northwest. However, they 
must be used with full knowledge of their limitations. 

1. The term “return interval”, when considering annual floods, represents 
the average length of time, in years, between the individual occurrences of 
flood peaks of a given magnitude over a very long period. It does not indicate 
that the flood will occur on a specific time schedule, 50-yr flood is one that 
has a 2% chance of occurring in any year. 

2. The relationships shown in this report were developed from data 
gathered on the natural flows of representative streams, and should not be 
applied to streams on which the flood-flow potential has been altered by the 
works of man, or to streams with unique drainage-basin characteristics. 
Many of the smaller streams in the semiarid areas of eastern Washington 
and Oregon are essentially dry except when they experience flash floods 
caused by local thunderstorms, This type of stream does not lend itself to the 

ethods of analysis presented herein but will require special studies. 

3. Extrapolation beyond the areal limits of the basic data is not advised. 
| Se mean annual flood formulas should not be applied to drainage areas 
outside the regions for which they were derived. 

4, The length of existing gaging-station records does not justify extending 
the frequency curves beyond the 50-yr recurrence interval, Extrapolation of 
he curves beyond this limit may introduce considerable error. However, 
uch extrapolation may provide the best available estimate. 

5, The confidence that may be placed inthe accuracy of the results obtained 
by using these data is related to the areal density of the basic data used. 

‘he denser the gaging sites, the more confidently the frequency curves can 
L defined, and the more accurately the geographical and runoff factors for 


‘he flood formulas can: be determined. 

6. Use of the methods outlined in this paper will facilitate design of 
tructures commonly built in or along stream channels. Where design re- 
yuires consideration of recurrence intervals greater than 50 yr, or the de- 
Bign is for a site otherwise outside the limitations stated, a special study 


should be made. 


BL 


. 


' SUMMARY AND CONCLUSIONS 


A method has been presented for determining the magnitude and frequency 
f floods to be expected at any site in the Pacific Northwest, with exceptions 
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as noted in item 2 of the limitations. Flood-frequency curves for individual 
stations are combined into several regional flood-frequency curves. The 
regions include stations for which slopes of individual flood-frequency 
curves are similar. The regional curves show flood magnitudes in terms of 
mean annual floods. To obtain the frequency curve for an ungaged site, it is 
only necessary to compute the mean annual flood for that site and convert 
the appropriate regional frequency curve to flow in cubic feet per second by 
multiplying the ordinates of that curve by the mean annual flood. 

Magnitudes of mean annual floods have significant correlations with drain- 
age area, mean annual runoff, area of lakes and ponds and in some places with 
elevation. After all these factors are taken into account in multiple correla- 
tions, large residual standard errors still exist. These residuals tend to be of 
about equal magnitude in given geographical areas, suggesting that a factor 
such as geology of the areahasa significant effect. So-called geographic fact- 
ors based on residual errors are developed. 

The formula for mean annual floods in western Washington is expressed in 
terms of drainage area, area of lakes and ponds, mean annual runoff, and the 
geographical factor. Similar equations have been developedfor other areas in 
the Pacific Northwest. Drainage area and area of lakes and ponds are easily 
obtained from topographic maps. Mean-annual runoff is obtained from an 
isopleth map based on gaging-station records. Geographical factors are ob- 
tained from a map that has been divided into areas having residual errors of 
similar magnitude. The geographical factor is based on the median residual 
error for that area. Because all of these factors are readily available, the 
formula is easily applied. 

The methods presented herein are subject to limitations, They cannot be 
applied to streams in which the regimen has been so altered by the works of 
man that satisfactory adjustments for such cannot be applied. Nor can such 
adjustments be applied to certain small streams in semi-arid eastern Wash- 
ington and Oregon in which practically the sole sources of runoff are local 
thunderstorms, Extensions of frequency curves beyond the 50-yr recurrence 
interval are not justified, and extrapolation beyond that limit should be a 
computed risk taken only when no better estimate is available. 
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PERIODICAL GRAVITY WAVE ON A DISCONTINUITY 


By Bernard Le Méhauté! 


SYNOPSIS 


A general theory is presented to study the transmission and reflection of 
sravity waves on a discontinuity to the first order of approximation. 

Particular attention is given to the cases of an obstruction, a change of 
width, and a change of depth. In the latter case only, the gravity wave arriving 
at an angle is studied. 


INTRODUCTION 
| 
: In harbors, the study of wave motion, caused by the action of incident waves 
‘rom the open sea, is generally a difficult theoretical problem. The complexity 
of wave patterns is increased by successive reflections within the harbor. 
Some authors2,3,4 gave linear solutions for some simple cases, assuming 
hat there were no energy losses in the harbor, that is no beach or wave-trap. 
hus, the solutions for circular,2 square2)3 and sector harbors connected to 
he open sea by a small opening are known. The general method consisted in 
| ie a potential function ¢, satisfying the usual free-surface, bottom and 
limit conditions for the harbor, and combining this with the potential function 
t sea. 
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Unfortunately such a potential function ¢ is always complicatedif there are 
singularities in the limits. (The reader can see, for example, the very simple 
case, studied by Ursell,’ of the wave transmission in infinite depth under a 
vertical wall.) 

Therefore, this kind of solution would not be expected to be of practical 
importance in harbor engineering. Furthermore, the solutions of many prac- 
tical cases are still quantitatively unknown. In particular, the effects of the 
harbor neighborhood shapes, and of the beaches inside the basins, are usually 
ignored. 

As far as the motion can be considered approximately two dimensional, 
there now exist linear computation methods more simple than the potential 
function method. 

F. Biesel and the author used the complex number computation method to 
study», 6,7 wave transmission and reflection on symmetrical obstacles and 
applied these results to study agitation between a combination of these ob- 
stacles, that is, in simple shape basins. In the Sogreah laboratory in France, 
more complicated shape basins that had not been Seca eee! analyzed were 
investigated. The results of these experiments are available. /,8 

The purpose of this paper is to present a more general wavetheory for 
non-symmetrical discontinuities using a similar method of computation. These 
theoretical results can be applied to establish a general formula of the agita- 
tion value, in a basin subjected to ocean waves, that will verify all the previ- 
ously published experimental results. Furthermore, this theory allows one 
to generalize the experimental results and to analyze, quantitatively, the effect 
of a beach in a harbor. 

Perhaps the calculus hypotheses may seem too simple for a rigourous logi- 
cal mind, but the relative simplicity of established formulas allows engineers 
to use them in practice. Above all, the assumptions are justified by close 
agreement between theoretical and experimental values. 

Notation.—The letter symbols adopted for use in this paper are defined 
where they first appear, in the illustrations or in the text, and are arranged 
alphabetically, for convenience of reference, in Appendix II. 


LIMITS OF STUDIED PROBLEMS 


The study is limited to periodical irrotational gravity waves. It is assumed 
that the fluid is ideal (except when specified otherwise), and any terms of the 
second or higher power of the wave steepness are neglected; that is, the square 
of the velocity (that is valid if the product of the steepness ne the wave and the 
cube of the inverse of the relative depth is small— (=) (# ) P 

The motion is mainly two dimensionaland will be treated as such. If there 
exists a third dimensional component, it must have secondary effects, and only 
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the main two dimensional components are analyzed. Finally, some details o 
three dimensional motion will be given. 

When a wave arrives on a two or three dimensional discontinuity, as show: 
in Fig. 1, the incident wave energy is partially transmitted, partially reflecte: 
(and sometimes partially destroyed by friction). In the close neighborhood o 
discontinuities, the waves become deformed, but the deformations disappea 
exponentially with the distance from the obstacle. 

In two-dimensional cases, motion is defined by a potential function that ha 
the form of a linearized sum of constant coefficients (generalized constants o 
Fourier) of which some have an exponential form and tend rapidly to zero 

In three-dimensional cases, the motion becomes two-dimensional at a shor 
distance from a discontinuity if the width of the considered domain is smalle 
than a half wave length. If the width of the domain is greater than this value 
Jeffreys’ waves are theoretically superimposed on the two-dimensional mo 
tion. Only in cases in which the width is large with respect to the wave lengtt 
is it necessary to consider their effect. This, then, becomes a pure three 
dimensional diffraction phenomenon. 

On the other hand, if three-dimensional discontinuities are notSymmetrica 
about an axis in the direction of the incident waves, resonant two-dimensionz 
motion can be excited in a direction perpendicular to the incident waves. 

In this paper, it is assumed that all these three-dimensional effects do nc 
occur, or are negligible, and the deformation of the waves disappears ex 
ponentially with the distance from the discontinuity (that is within one or tw 
times the depth). The following is a consideration of the two-dimensioné 
periodical wave motion. 


DEFINITIONS 


Consider a discontinuity in a wave flume (such as one of those in Fig. % 
separating the two domains (1) and (2). Incident waves from domain (1) ar 
transmitted to domain (2) and reflected in domain (1). 

Sufficiently far from the discontinuity for the waves to have the character 
istics of a two-dimensional periodical wave, it may be assumed that the rat 
of the transmitted wave height to the incident wave height has a value a that i 
independent of the absolute value of the wave heights. 

Similarly, it may be assumed that the ratio of the reflected wave height 1 
the incident wave height has a value 8. These hypotheses follow from the as 
sumption of the linear theory. - 

On the other hand, the transmitted waves are out of phase with the incide: 
waves by an amount @ Similarly, the reflected waves have a certain phas 
difference from the waves that would be reflected by a perfectly vertical wa 
placed in the same position as the discontinuity. It is assumed that this pre 
vious result is applicable for a discontinuity that is not strictly located. Bu 
in this case, it is necessary to introduce an imaginary reflecting referent 
plane related to the discontinuity. It is then possible to designate the phase | 
the waves as being that of their vertical amplitude in front of this referen 
plane, whatever their direction of propagation. For a theoretical plane 
rier, it is natural to choose it as the reference plane. When the discontin 
possesses a plane of vertical symmetry perpendicular to the wave direc 
this plane would be chosen as the reference plane. 

Consequently, the motion of the free surface atany point can be represe! i 
by a vector running inthe plane of imaginaries (Fig. 3), analogous to that us 
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n the study of alternating currents in electricity. Therefore, we can repre- 
sent the coefficients of transmission and reflection by vectors, that is by the 
somplex numbers @ and 8, representing these coefficients both by their argu- 
nents: (@and @) and by their moduli ( @ and 8). 

' When the phase and amplitude of the incident waves are represented by the 
somplex number: i A, the transmitted and reflected waves will be represented 
y @ Aand BA, respectively. 

If the discontinuity is symmetrical about a plane that is perpendicular to the 
vave direction and separates two identical domains, it is clear that the value 


Pl 


FIG. 3.—VECTORIAL REPRESENTATION OF A PERIODICAL GRAVITY WAVE. 


‘the coefficients @ and B are independent of the side from which the incident 


However, in the general case, @ and B are probably different for the two 
rections of incident waves. Therefore, we have to consider two transmission 
efficients @, and @, and two reflection coefficients 1 and £2. 


Consequently, a discontinuity is characterized byeight parameters: Trans- 

‘ a 

4 - 1 eed i ; 

ion coefficient in domain (2)—a, ia ; Reflection coefficient in domain 
B 1 i a> ‘ 

-By { A ; Transmission coefficient in domain (1)—a, d, ; Reflection co- 


eee 
ient in domain (2)—Bo { ‘ i 
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From the continuity and energy conservation principles, relationships be 
tween these eight parameters will be established. The reader who wishes t 
use directly the results of this study will refer to Tables 1 through 5 in whic 
the practical results of this study are given for a number of common cases 


THEORETICAL RELATIONSHIPS 


All the mathematical developments and demonstrations are given in Appen 
dix I. 

Let Ly and Ly be the wave lengths in the domains (2) and (1) respectivel, 
and ly and 1, the widths of these domains. Also, 


L 12 
Day Fp tee os theron gten ene ( 
From the continuity principle it is found that: 
Z Qy Bs By MR 5 gig od 6 bees o) 2) wl ec ee ee (2: 
and ws = 
7, Og, Bg = hes ta ens oa (2 


in which i is a complex number of modulus: i = 1 and argument: te 0 
Let hg and hy, be the depths of the domains (2) and (1), respectively, an 


m = 2a and m = 27 ( 
2 Lo 1, Ly ¢ &. €-0.0 20. 0 oh ae 
and, to simplify the writing 
2m9h 
i) Se eS 22 
sinh 2 mg ho 
A = | ———~————_ |... ..... ... eee ( 
2 my hy 
sinh 2 my, hy 


It is pertinent to note from Eq. 4 that A = 1 in the cases in which (1) t 
two domains, (1) and (2), have the same depth, (2) there exists shallow wate 


wa = : : 2mh 
since: sinh 2 m h*= 2 mh, and (3 = | 
, and (3) deep water exists, since sinh mn 
to zero. : 
The conservation of energy principle leads successively to the equaliti 
2 208 
AZ ay + By Sh ea) sashe! ass een Aen ee q 

1 #) Pie 

iz Ap + Bo = Levee rere cree ee ee 
and e 
A Atal ; 
Q - By = -(@) - BL) +7 2nz) OO" i eae : 
in which n=0,1,2... . 


The two vectorial equalities, Eqs. 2a and 2b, equivalent to four algeb ‘2 
equalities, form with Eqs. 5, 6, and 7 a system of seven independent equatic 
with eight unknowns. : 3 

Hence, we are able to express the values of seven parameters as functll 
of only one, a1, and the factors Z and A, that characterize the two domai 


— 


iy 9 PERIODICAL GRAVITY aL 


md do not depend on the exact shape of the considered discontinuity. Conse- 
juently, it is sufficient toknow the exact value of @, aS afunction of this shape, 
n order to know the others from the following theoretical formulas Ge Zya)): 


Qs = AZ a, Sree SOS, Ey oS aes (8) 
1/2 
s i 2 
B= By (1 AZ 2 CaS. aa aa (9) 
Ne " - S 
cos a = cos @ = > (AME). ee oh dws (10) 
2, 
2- Zaf(A+ Z) 
NN 
cos 8) = ee ae ee (11) 
2 1/2 
2(1 “ A Za?) 
nd 2 
2- a; (A+Z) 
cos By = ee a ae ee (12) 
ofidize 2/2 
(1 - AZaz) 
‘urthermore, it may be shown Pah ~ 
By +B 
Te ec 2 “ile 
Q, = & = 5 +5 Ge A a sa (13) 


GRAPHICAL REPRESENTATION 


_ Assuming a, is known, we can easily derive Eq. 10. Now consider the half 
ircle of diameter AB (equal to unity) shown in Fig. 4. On this circle, we 


raw BC = cos @), or the angle Kee = . On the segment BC, we plot the 


pint D, so that BD = A Z a; and draw the circle of center A and radius AD, 
itersecting CD at the point E. It is easy to verify that this sketch satisfies 
ll the above conditions. Thus, we have 
Shira Hee Hil a Siac SP hen We Pa he (14a) 
BAD = Be re eet cee (14b) 
ea gic ena (14¢) 
Yeap ss 
I yyceeae oe Aang eM Denies (14a) 
UES SAS Se eee Oe eee: (14e) 
AE EB cis ns es wns os (14g) 
DB = Za, = = ER ye ee eS (14h) 
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FIG, 4.—VECTORIAL GRAPHICAL REPRESENTATION 
WITHOUT LOSS OF ENERGY, 
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1 1 
ap UTA % meteys tary? ag Qy 
SS RE TT eT tRNA) oe (LAD) 
r more simply 
A = 8 RACES OR GRR Pens Se (15a) 
D = Z dy =.05 RE ae a en Se aR (15b) 
—— 1 — —> 
AE= > RO c,h sieeve og Md nae (15c) 
ee 8 
PE ag et eR ee a ee (15d) 
2 
Aa; +—a Za, + 
— rey i oe 
CB =——,"— = —__4 ~- 2 (a+2).... (15e) 


me can see also, without further computation, the following noticeable rela- 
onships: 


A a + Bo Biri SE Ces A RE} PELE 2 ee (16) 
a5 2 — | 
A Piet sitcihen dies ovo opus aoe (17) 
nd 2 


§ SOME PARTICULAR CASES 


The relative variations of each coefficient are easily studied by means of 


ig. 5. When A Z is greater than 1, D and E, fy and fg, Ay and Ay ...are 
sversed. (a) If A, Z tend to 0, either by a reduction of domain (2), or by an 
crease of domain (1), as shown in Fig. 5, then, according to the sketch for 
1e wave moving from the domain (1) to the domain (2) 


A 
: A 


7 
t is, there is total reflection without change of sign (as in a closed acoustic 
pe); for the wave moving from the domain (2) to the domain (1), the sketch 


ows 


is 


& A — 0 a — 0 
re is also total reflection with change of sign, (as in an open acoustic pipe). 


(b) If AZ tends to 1, that is, if the domains (1) and (2) are identical, the re- 
ts previously published are found, that is 


By = Bo 


20 November, 1960 HY | 


Even if the discontinuity has an unsymmetrical shape (Fig. 6), the effects 
onwaves are symmetrical. (This is due to the assumption that the discontinu- 
ity is local, and that the linear theory is used. ) 

However, this result is also true when domains (1) and (2) are different, 
with A, Z = 1, (for example, as in the case shown in Fig. 7. But A, Z can be 
equal to unity for only one value of the wave period. 

(c) If a, tends to 0 (by obstruction), A, Z being constant 


a5 == 0 
A AN as 
a, and a ——s 


3 By and By — 1 
AN AN 
By and By — 0 


There is total reflection in the two domains without change of sign, that is 
there is total obstruction. 


VALUES OF A AND B AS FUNCTIONS OF THE SHAPE 
OF THE DISCONTINUITIES: EMPIRICAL FORMULA 


It has been seen that all the coefficients @ and B are known when we ar 
able to express only one (for instance a1) as a function of the shape of the dis: 
continuity. This value has been obtained with the help of theory, physical con: 
siderations, and experimental results, some of which have been previoush 
published.7,8 (However, the formula obtained can be modified after othe: 
complementary tests or further theory.) 

Consider the most general possible shape of discontinuity frequently me 
with at the mouth of harbors (Fig. 8). : 

The coefficient a; is given by the empirical formula (N) 


1/2 /1,\1/4 : 
-(& adel pebeees | 
ay (i) Ce) Sa iy A.M Seb Rot (19 


Lf 
in which b is the opening between domains (1) and (2). However, this empirica 
equation could be justified by a number of theoretical and physical consideay 
tions. 
(a) It is easy to verify that, if we apply the theoretical formula of Eq. 8 w 
find 3 


rf, 
the function a» is perfectly symmetrical with respect to Qy. This is physicall 
consistent. 


(b) If 1y = lp = 1, and hy = hy = h, (Lg = Ly, A Z = 1) we obtain 


oe 
; 


ay 9 


ae 
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FIG. 7.—A CASE WHERE AZ COULD BE EQUAL TO UNITY DE- 
SPITE THE CHANGE OF DEPTH AND WIDTH. 


FIG. 8.—-OBSTRUCTION, CHANGE OF DEPTH, CHANGE 
OF WIDTH; NOTATIONS. 
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1 - p\1/2 
piss Bess (ee Sp oe ee ee (22) 
cos a = ay @, = Q, ee Ne ete oes (23) 
Antal’ 
cos By = Bo By, s Bo = B Pit ee ee ee eek (24) 


and so on. This satisfies the hypothesis that the transmitted and reflected 
energies are proportional to the opening and obstruction, respectively. This 
is also consistent. 

(c) If 1 = 1, = 1g =), and hy # hp (Ly # Ly) we obtain 


re 2 
a = Ege Mp ee (25) 
1+A— 
Ly 
and 
Lo 
DP sass wither adoeh lace (26) 
Ly 
On the other hand, if we consider the formula 
a 
cos ay = ~ Cava Zeta: BSS hehe eane (27) 


and the value a; in the case of a simple change of depth: 


2 cos a 
dee NE EI ro (28) 
L2 ‘ 

Pha taper 

Ly 
Equating these two values leads to A and Z = 1, or hy = hg when ay = 0 that is 
consistent. 5 
In the case of a discontinuity with deepening and without obstruction, as 
shown in Fig. 9, when A Z varies from 0 to 1, that is for instance, when the 
dimensions of domain (2) vary from 0 to those of domain (1), a, varies from 
1 to 2, By and By from 1 to 0, and a from 0 to 1, at the same time. Conse- 
quently, the wave in the domain (2) cannot have a height greater than twice the 
height in domain (1), and this occurs when the ratio of the dimensions of the 
two domains tends to zero; that is when i — 0. 
It is interesting to note that the wave height always increases when the 
depth decreases suddenly. That can explain the great agitation of shallow 

water close to deep water. 


Finally, it should be noted that in shallow water L = TV gh,and: sinh 2m P 
=2mh. Hence, A =1 and 4 
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FIG. 9.—A CASE OF MAXIMUM TRANSMISSION, 


FIG. 10.—VARIATION OF THE DISTANCE BETWEEN OR- 
THOGONALS FOR A WAVE AT AN ANGLE, 
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9 and later 


The formulas of this particular case were obtained by Poincare 
by Lamb. 10 

(d) In the simple case of a change of section (Z # 1), without obstruction, 
the opening (b) is equal to the width of the smaller domain. For example, if 


the smaller domain is (2), b = lg, and one obtains 


1,\1/4 
+e Bek oe ok ee (30a) 
1+A cco 
and 1 
1 3/4 
ee ye SSS ee (300 
ao \h 1 ly 
TPE Set 
A Lo 


Some authors, after Lamb, 10 have adopted a formula that, with the present 
notation, leads to 


= j 
ay ie WPS Ae ae age PP ee (31a) 
and 2 
ay Totti tt tees (31b, 
is; zm 


These two last formulas do not verify the theoretical relationship demon- 
strated above as Eq. 8 (with A # 1) holds true. 

On the other hand, it is difficult to introduce the effect of an obstructior 
between the domains (1) and (2), in this formula, without verifying both the 
theoretical relationship such as Kq. 8, and giving a correct formula for any 
particular case. 

Indeed, if for example, one admits, as it seems logical that 


b 1/2 2 
ay - & Tp eres: pee oko (32a, 
and 
Sh 1/2 2 
eo) aoe 1 +s © oe + ie tele’ is) ‘ew Metce (32b, 
2 ie arms 


the fundamental theoretical relationship of Eq. 8 is not verified. In practice. 
the difference between the two formulas for (a1) (and (aq) ) is not very im- 
portant, but it would seem better to keep the proposed formulas that Satisfy} 
both the theoretical relationship and the physical considerations for any par- 
ticular case. 

Hence, the general formula of Eq.19is valid for any case, that is: obstruc- 
tion without change of characteristics of domains (1) and (2) or change of char- 
acteristics without obstruction. 

Introducing the value of a, given by Eq. 19 in the theoretical formulas (Eqs 
8 through 12) permits us to compute the value of the other coefficients; a, By 


9 “Méchanique Céleste,” by Poincaré, Torne, 131, p. 87. 
10 “Hydrodynamics,” by Lamb, 6th ed., p. 263, 
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\, Qo, a, Bo, Ay, in the general case and, in a number of particular cases, 
is it is shown in Tables 1 through 5. 


WAVE ON A DISCONTINUITY AT AN ANGLE 


When a periodical gravity wave arrives obliquely at a sudden change of 
lepth, a part of its energy is reflected. The angle of the reflected wave with 
he line of discontinuity is the same as the angle of the incident wave, as is 
shown by Fig. 10. The distance between orthogonals does not change. 

Because of the change of depth, the velocity of the transmitted wave changes 
suddenly. Hence, the transmitted wave is refracted, and the distance between 
yrthogonals varies. For example, if the distance between orthogonals ly is the 
leeper domain (1), and 15 is the shallower domain (2), we obtain 


m which 6; and 59 are the angles of the waves with the line of discontinuity 
m the domains (1) and (2), respectively. On the other hand, because of the 
thange of velocity of the waves 

i 


Lo sin 69 


lp 
Combining these relationships, the value of ae is obtained as functions of 
1 
=— } and one of the angles 5; or 69, that could be chosen as the angle of the 


L 


1 
= wave 


Leh ee lik 

i at Ly soe cos 69 

8 4 ae ae EE ae ss (35) 

1! cos 5, ij 1/2 
2 


If the continuity and conservation of energy principles are applied to this 
henomenon, all the previous results are valid, provided that the ratio of dis- 


neces between orthogonals|—} takes place, instead of the ratio of width of 


1 1 : j 7 
domains (1) and (2), (3) Hence, without further computation, introduc- 


s the new value i 


aR ee ee (36) 
L, L 
_ the previously demonstrated theoretical relationships, are valid 
AZ ob +, pi-= 1 oO ee ee oe (37) 
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and soon... giving 


and so on. 
In order to express completely the value of a1, it will be noticed that a re 
duction of depth causes an increase of distance between orthogonals. Hence 
this case is similar to the case shown by Fig. 7. 
Hence, if the domain (2) is assumed to be the shallower domain, the value 
ay and a9 can be written 


1, \3/4 2 
a) =| = —  eeaed? i eeeee (40: 
1 ' Jc 
1 2 
2 t+A — 
Ly 
and 
"\1/4 
15 \K/ 2 
a <=) ————_ ........... wa. (40) 
2 i 16 
1 1 - LS a 
Lo 
1 
Introducing the values (*) leads to 
1 
cos 54 3/4 2 


ay = i fhm hate R i iw rate bees (41: 


and 


R 
i] 


D) f(a allo, Yee Lhe « We OW aS (4 1] 


eee 
7 La eae 


Since the domain (2) is the shallower domain, Lg is always smaller than L 
and a, always has a definite value. But Q@_ tends to be infinite if the a 
arrives at an increase of ie at an angle 69 such that (sin 6 2) tends to (G 
In the same time, 6» tends to 5 amy At (sin 52) is greater than oe , that is sin 6 


> rot there is total Sette at the discontinuity. In AB such a war 


diffracts inthe deeper domain. The anology of this phenomena with light arrit 
ing at a medium of a different refraction index is evident. 

This particular solution has been considered to protect the mouth of a ha 
bor from agitation, but it does not offer a very practical solution because 
the cost of maintenance by dredging. 


t 
a 
iS 
E 
The most complete verification of this theory has been made by studyin 
the agitation in a basin, that is, between two discontinuities. Indeed, a theo 


EXPERIMENTAL VERIFICATION 


O_o 
11 “Etude de la houle en profondeur finie-Note sur la réflexion totale” by Carlott 
La Houille Blanche, Mars and Avril, 1947, pp. 112-116, 7 
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FIG. 11.—-EXPERIMENTAL RESULTS IN THE CASE OF A CHANGE OF WIDTH. 
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retical formula for the agitation in a harbor will be established. This theo 
retical formula is based on the theory for waves at a discontinuity presente: 
here. If such a theoretical formula is well verified by experiments, it is con 
sequently true that the theory for waves at a discontinuity, the basis for thi 
further computation, is also verified. 

The results of some tests made at Queen’s University in Kingston, Ontario 
are presented in Fig.11. It is important to recall that the theory is valid onl 
if the quadratic convective inertiais small, that is if (2) (z) is small. Hene 
all the tests have been carried out with small steepness and large relativ 
depth. It must also be remembered that, except in the case of a simple chang 
of depth, the width of the channel must be smaller than L/2. 

In addition, it has been assumed that the loss of energy due to friction i 
negligible. The main way to achieve this, is to maintain conditions prevailin 
for negligible convective inertia, that is, a small wave steepness and a larg 
relative depth. Unfortunately, since the discontinuities used for these test 
were made of simple concrete blocks placed in a wave flume, some energy wa 
lost because of their roughness and because ofthe leaks occurring between th 
blocks. This loss of energy is evident in the case of a complete closure of th 
wave flume, because the reflection coefficient is only between 0.8 and 0.9 in 
stead of 1.0. However, in any case the loss of energy, computed by the form 


ulas 
2 i 2 
- (AZat - 3)| i Seer a+ 2) 

always stays very small, as long as there is no breaking of the wave at an 
point. : 

Fig. 11 shows the results obtained in the case of a simple change of width 
The theoretical curves proposed by Lamb are compared with the theoreticz 
curves derived here and the experimental results obtained. Although not bein 
definitely conclusive, it seems, according to these results, better to adopt tk 
values proposed in this paper. Indeed, the fact that a certain amount of energ 
was lost by friction against the rough concrete blocks and by leaks in betwee 
them must be considered. This especially effected the reflection coefficien 

Another occurrence of discrepancies was caused by the fact that the trans 
mitted wave was partially reflected on the terminal beach of the wave flum 
This effect is not negligible if the steepness of the transmitted wave is smal 
as in the most frequent case. Hence, it may be concluded that this theory he 
been verified within the accuracy of the tests. However, it would seem relz 
tively simple to obtain some better experimental results. 3 

Some other points such as the effect of an obstruction, Eq. 21 or 


erraaenrees ae 


or a change of depth, Eq. 25 have been verified with the same degree of at 
curacy. : 
No tests have yet been made to verify the value of the theory of obliqt 
waves at a discontinuity. The author would be happy to know if there exist al 
experimental results that demonstrate the validity of his theory and furth 
tests to verify or to improve his empirical formula (Eq. 19) giving the va 
of ay. a 


, 
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It seems that many theoretical results are valid, even when the width is 

a : 11 \ / 12 

; dimensionless factors such as: |— |,{—“-] could be taken 
Ly} \Lg 

into account, with the condition that there are no dysymmetrical effects that 

may cause a transverse resonance. 


greater than = 


LOSS OF ENERGY 


In the case of a loss of energy at the discontinuity, the relationships ob- 
tained from the continuity principle (Eqs. 2a and 2b) are still valid. However, 
because of the energy loss 


2 
AZ of ot 


ah 2 
eats 


In this case, the graphical representation has to be modified as is shown by 
Fig. 12. The theoretical limits are obtained when a, and ag tend to zero, and 


La Le 


FIG. 12,—VECTORIAL GRAPHICAL REPRESENTATION WITH LOSS OF ENERGY. 


the theoretical maximum loss of energy on a located discontinuity is equal to 


1/2, when 
1 1 
: VAZ ay =—5- “BL = jis eplammee a3) 
L 
ao 1 


1 
a = stelibehe GLa cee MSS 
Yaz : B27 = (44) 
Physically, this case could be conceived if the discontinuity is a pervious 
obstruction such as a rockfill breakwater submitted to long waves. 
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TABLE 1 
Domain 2 | Domain | 
Z . Incident 
a, ) —— 
(ee THEORETICAL GENERAL 
I ATION 
Wibidant 7 i) F pater FORMULAE 
a eee Ot a 
Be | f( A, 23 o;) 
Coefficient in case of Lots 
Change of Depth & Width. Ly {; 


ae : 2Moh 
Coefficient in case of i+ Sink Ohta 
2m, h 
| * sinh 2myh 


Change of Depth. 


Coef of Transmission 


in Domain 2. 


Coef. of Transmission 


in Domain |. AZa, 
Phase of Transmission COS a,(=cos hal 
ad > (A + 2) 
3 V2 
Coef. of Reflection B, (=g,) C= SA a) 


Phase of Reflection 2 - Za® (A +Z) 


cos B, 
in Domain |. 2(1- AZ a*)'72 


Phase of Reflectio 
a 2.— rae (Ae Z) 


cos P. PPE SE AE: 
in Domain 2. Be 2(1 - AZa® )! 


GRAPHICAL =e 
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TABLE 2 


Obstruction Change of Width Obstruction 
Change of Width 


hy = he=ah 


taatal 
falta)” 
ale) GTA 
tc 
2 Leste) ~) 


2B, 


2 BIT 


2B 


Same as_ General 
Case. | 
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TABLE 3 


Change of Depth Obstruction and Change of Width 


and Depth 


[sy 2(A+Z) 
te ne 
Seah 


/2 
az [at 24 AZ 


[i+ ate] 


2 —Za?(A+Z) 


A, = 0 
2 A 
)— 2 BZ(A +2) 2 - atA +Z) | 
cos Ba Lo 7 —— et ¢ 
[1+ atz] 2B ‘ 


Same as General Case. 


Same as General Case. 


(1/z)a 
Ae 
Sateeey 
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TABLE 4 


Obstruction Limit Case Limit Cases Pervious Wall 
Change of Depth f.—- 0 {, =f. -f 
and Width ho Oo hy=he-h 


iA 2 


rae 


1+ Ate 


e2 

: 

50s a, 

: ay 

=cos @ paras a) 


— | 


Reflection 


Total 


eee | 22 Zawasz) 2, 0 Ay = a+ B=7 
Ai 2B, loop with a loop a? + B7<\ 


Nn 
Bo=zm 
with a node 


2 -9°(A4Z) 


28, 


with a loop 


Same as General 


Case. 
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TABLE 5 


Pervious Wall & Wave at an Angle 


Change of Depth 


Change of Depth & Widin 


1+ — 3. 
sinh 2Mo 


Same as Previous Case. 


cos 82 i 


2 
L, : a} iy \ cy 
[1 ie {i sin [eee ihuais : 


3A +2’) 


3/4 
cos 8, 2 
| | 
Lice 24/2 L2 
fone [1 - {<2 sin a} a ion 


cos a) 
(=c0s a 


L 
i ee 
If sin3.> GC 


By =| | 
Total Reflection 2 


B; = Bo 


Za, + B= 


(1/Z) a+ Bo =7 2 
cos 8, ar Fe Bgl Aten 
AZa’ + Bi <I 2B 
2 
— LP bey 2- Zaj(A+Z) 
2 Az 22 marae ty hoa 
: ; 
c 
‘ 
Same as General Case. E 
" 
: 
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APPENDIX I. MATHEMATICAL DEVELOPMENT 


Continuity Principle.—We will assume that the discharge of incident and 
flected waves in domain (1) is equal to the discharge of the transmitted 
yvave in domain (2) on the discontinuity at any time. It is natural to make 
his assumption, for if this were not the case, a quantity of water would suc- 
essively be accumulated near the discontinuity and would issue from it. 

The discharge is the product of the velocity times the cross sectional area 
depth times width). The velocity is a function of the vertical abcissa y. 
lence, the discharge could be written 


=h 
l F VEY VST Ce Se oetane pag aia ore ee mcs pene (45) 
Vee 


n which 1 is a constant width and h the depth of the considered cross section. 

Now assuming that the domains (1) and (2) have constant widths ly and lg, 
nd constant depths hj and hg, respectively, the continuity principle is given 
y equating the value of the discharge in the two domains at the discontinuity. 
“hat is 


hy hg 
lif u(y) dy = ly f Ho (yey. sae snes (46) 
Oo 


.ssuming that the motions in domains (1) and (2) are defined by the potential 
unctions $1 and $9, respectively, and that the horizontal x-axis origin coin- 
ides with the discontinuity, uj (y) and ug(y) could be expressed as a function 
f $1 or $9 as follows: 


o¢ 
uj(y) = ea At hig eg iegranteoa) 


a¢ 
ug(y) = ( a A la .  (46b) 


lence, the continuity relationship becomes: 


O¢4 cites fee 
i ly ales dy =p f Lee patty, Deni (47) 


“The potential function $1 is given by the sum of the potential function of an 
cident periodical gravity wave of amplitude 2a travelling towards the dis- 
ntinuity and the potential function of a reflected wave, of amplitude 2ax BI 
-avelling seawards and having a difference of phase By at the discontinuity. 


Beer nL 20 
ce, $1 may be written |h = +7 m= TL 
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h cosh mj (hit y) 
id neha mj. sinh mj hy 


h A 
S pingnic se cosh in (014 1) cae ty eens 


my sinh my hy 


sin (ht - mj @) 


The potential function ¢9 is that of a wave of amplitude 2a X Q4, travelling ir 
the same direction as the incident wave and having a phase difference by wit 
the incident wave at the discontinuity; 


h cosh m9 (ho + y) 
m9 sinh m2 hg 


04 cosh mj (hi + y A 
eae 2 poosh mai | + 2 my hy [cos ht - 81 cos (nt+8})| -. . (5m 


ech cosh mg (hg + y) | 
ce Ae =a sinh m9 ho [24 cos (ht + 2) depe 6 take (51, 


Introducing these values in Eq. 46 since 


d9 = -aal sin (ht - mg a+4y) rs 


Hence, 


and 


h 
cosh m(h+y). .. tL | 
, Binh ath ae We eae (52 


fe) 
we obtain 


ly Ly [cos hte-924 cos (ht + B1)| = lo Lg aq cos (ht + a1) 


Assuming that ht = 0 andht = > successively, leads to 


11 Li [1 - 81 cos By] = lg Lg a1 cos ay oo 2 (53a 
A 
ly Ly [o - By sin By | = Ip Lg ay sin ay Paar A (53b 
or 
A TS ES scalp 9 nd sirefeo te nen ime (54 
in which 
lg Lg 


Similarly, considering an incident wave from the opposite direction we fin 


Soe hey p 
Z, ag + Bo i qe nalts bee bisagg Sal 


Already, one can see the eight parameters will be functions of the pecil 
(included in Z), except in shallow 


lo 2 ‘ 
water, Z Mis hy , and in infinite depth, or, if the depths of the two domair 
‘ 
4 
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lg 
ure equal Z yy 


Conservation of Energy.—The transmitted energy of a periodical gravity 
wave of amplitude 2 a and with a length of crest equal to 1 is 


1 LA9 2mh 
5 Pe aol ( ee.) Ascoe OS 255 (57) 
er unit of time. 


On the other hand, the most important loss of energy at the discontinuity is 
jue to the convective inertia. Since this force is negligible, one may neglect 
he loss of energy and equate the incident wave energy to the sum of the re- 
‘lected energy and the transmitted energy. To simplify the writing we intro- 
juce the new symbols 

2 my hy 
Ay= 1+ oh2my hy 

2mg9 hg 
Ag=1+* inh 2 mp ho 


[he conservation of energy principle may be written as follows 


Incident energy = reflected energy + transmitted energy 
2 
a2 1, Ly Ay = (a 61)? 1, Ly Ay + (2.21) "1g Ly Ag 
Dividing by a2 1, Ly Ay and introducing Z and A = Ky yields 
| NE tea 1) ee nnn: (59) 


Similarly, for a wave from the opposite direction 
1 2 Vee 

Kz%2 + oD ede Se (60) 

Now, if we consider two waves of the same period but arrivingfrom the two 

pppostie directions, with height 2 a, and 2 ag and withany phase difference, the 

conservation of energy principle takes the form 


ao ae 
Ay ly Ly [81|“+ Azle be [2] "= 


i? 


enna LA be eaten 
Ar ly Ly [3 By + Bo % [7+ Aale be [81% +828 |” .--~ (60) 
hat is 
‘ 


a4? (Ay ty by 61? + Ag te be 2) + ag? (Aq lt Ly Ba” + Ag le Le a9”) 


s ey ae 
Ay ly Ly |@y| 2+ Ag lp Le |@2|~ = 


AS AYE ctASe IN 
; ~ 2 Ay ly Ly ay ag By 2g cos (24 + By - 22 - @2) 


DS a IN 
= 2 Ay Ig Lig ay 29.04 Bg 08 (4 +44 - & - Ba) Aer aeer eet n (Pd) 


38 November, 1960 HY § 
that can be written 
2 a 


AZ 


AGE WAT A’, 
- 2a, a9 [Ay ly Ly %p By cos (ay - 89 + Bo 
NN 


1 2 
Ay ly Ly ay? (A Z ay? + By) + Ag lp Lg ap” (az 22 +B9 ) 
AN 
- &) 
LP PN 
+ Ag ly Ly a, Bg cos (2 - a9 +Q - Bo) eine plas eee (63) 
Simplifying with respect to Eqs. 58 and 59, we obtain 
4. ght aA aE 4. 
Ay By cos (4) - a + By -4))+AZay Bg cos (aj - a9 +@, - Bg) =0... (64) 
A 
This relationship must be verified for any value, (4) - a). 
Assuming that 4\ - a, = 0, 4\ - a> = 3 successively, we find 
gig? & eee a Oe 


Gg Py = Bo - act. a. 20 fo vs ee wbalee vs (65a) 
or 

ay - By = -(& - by) +0 (+207) Pe er (65b) 

Akar ss 


(The solution a» - By =B9 -O1(+ 2 n 7) must be eliminated, a1, 81, a9, 62 al- 
ways being positive.) We also obtain 


as Pye AD ay Bo oP et Cue ee (66a) 
or P 
2 1 

Rae AD ke IAG UES A SE ARS Sa re 66b, 

Bo a5 (6b) 


The two vectorial equalities, Eqs. 53 and 55 being identical tofour algebraic 
equations, the system formed by Eqs. 53, 55, 58, 59, and 65b gives seven in- 
dependent equations with eight unknowns. (Indeed, it is easy to demonstrate 
that Eq. 64b depends on the others.) Therefore, it is possible to express the 
eight coefficients: a1, @y, By, B, a9, @, Bo, is as a function of only one (for 
instance @1), and there exist some other noticeable relations. From the pre- 
ceding equalities involving energy conservation d 


7 


2 ay ? 
=gey os Z cae & 2. @ 6S) se wae 6 @ 6.10) 6 eeee j 

Bo 4 By a 

and : 
; ‘ 

A Z ay? su By2 = ECA a” - Bo? = uf Cy eee (68) 

we deduce easily: ‘ 


a =AZay ee Ly eee ae ee (69) 
and to ; 


Bt = By =(1- AZ 042)? piety ans ae wate, Tes 
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Continuity Principle and Energy Conservation Principle.—From the conti- 
nuity relationship we find, seccessively, for the first incident wave (Fig. (13) 


Z2 a2 = 1+ p12 - 2 By cos fy chal ihe vee ethos dak (71) 


On inserting Eq. 69, deduced from the conservation of energy principle, this 
yields 


Aun BeBe? (A-+Z) 


EG Belo a oe a a ae ee (72) 
2 (1 -AZ ay”) 
On the other hand 
By2 = 1+ 22 ay2-2Za4 cos ay ... Sette. exors a Ata) 
yields 
cos AY = (A+ 2) aera eR ae (74) 


_ Similarly, for the second incident wave, from Eq. 55 and from the Eqs. 68 


FIG. 13.-VECTORIAL REPRESENTATION OF THE CONTINUTIY PRINCIPLE. 


and 69, we find 
A 2-042 (A+ Z) 
cos Bg = 3 1 ° ry # 440 4.6 8 (75) 
2 
7 2 (1- AZ a2) 
and , 
cos & = (A+ Z) Ia RIGS, BO te are Tec Le me x. (76) 


cos 4) = cos dy ihabiodcbte debbaia te 


iS 5) Peder te no een tO) 
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The solution ay =- & is not possible. Indeed, in this case, Eq. 64b would 


yield By =- bo +m (+2n 7) and cos Ay = - cos by. Now, from Eqs. 71 and 74, 
we could establish a relation between a; and Z and A, and all the parameters 
would be known only as functions of Zand A. This is physically impossible. 
Indeed, these parameters depend not only on domain characteristics, but also 
on discontinuity shapes, 


APPENDIX II. NOTATION 


a coefficient of transmission at the discontinuity: ratio of transmitted 
wave height to incident wave height; 

ay : coefficient of transmission from the domain (1) to thedomain (2); 

A : coefficient of transmission from the domain (2) to thedomain (1); 

B coefficient of reflection at the discontinuity: ratio of the reflected 
wave height to the incident wave height; 

By : coefficient of reflection in the domain (1); 

Bo : coefficient of reflection in the domain (2); 

@ : phase of the transmitted wave by comparison with the incident wave 
at the discontinuity; 

ay : change of phase relative to the coefficient a1; 

a, : change of phase relative to the coefficient a9} 

® : phase of the reflected wave by comparison with the incident way 
at the discontinuity; 

ay : relative to the coefficient By3 

BS : relative to the coefficient £9; 

a : complex number of modulus a and argument Qs 

B : complex number of modulus @ and argument ». 

dy : complex number of modulus: i = 1, and argument: ee 0; 

2a : incident wave-height in domain (1); 

4\ : phase of the incident wave in domain (1); : 

2a : incident wave-height in domain (2); 

a : phase of the incident wave in domain (2); 

$4 : potential function of the motion in the domain (1); ¥ 

bo : potential function of the motion in the domain ( 2): ; 

Ox : horizontal axis; : 

Oy : vertical axis, y=0 at the still water level of the free surface; ) 


subscripts: 


em: |e 
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: period; 


2 
i er 27h\, 
wave length: (. on tanh Ben ) 


: wave-height of the incident wave; 

: depth of the water; 

: width of the domains; 

: width of the opening of the discontinuity; 
: distance between orthogonals; 


: angle of oblique wave crest with the line of discontinuity; and 


Lol 
z= 22 2% 
Lily 
2 mo ho 
Yi ee 
kar he: sinh 2 mg ho 
Ay ‘ 2my, hy 


sinh 2 mq hy 


(1) relative to domain (1) 


(2) relative to domain (2) 


/ 
+h m 4 * Ree . 
r 4 
ph 
¢ bat 4 a Cok ve 
: ; ef era t 
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WATER EDDY FORCES ON OSCILLATING CYLINDERS 


By Alan D. K, Laird!, Charles A. Johnson2, and Robert W. Walker3 


SYNOPSIS 


Neighboring cylinders may increase drag force on vertical cylinders by 
75% or decrease it to mins 20%, and may triple lift forces, for spacing to 3 
jiameters and cause significant deviations to 10 diameters. Neighboring cyl- 
inder eddies cause large force fluctuations. 


INTRODUCTION 


: Recent investigations 4,5 have suggested that eddies formed along piles in 
»cean waves may make important contributions to wave forces on piling. The 
sresent laboratory investigation of these eddy forces is part of a continuing 
esearch project being conducted at the University of California in Berkeley. 


4 


a ———$—— 
Note.—Discussion open until April 1, 1961. To extend the closing date one month, a 

ritten request must be filed with the Executive Secretary, ASCE, This paper is part 

‘the copyrighted Journal of the Hydraulics Division, Proceedings of the American So- 

ety of Civil Engineers, Vol. 86, No. HY 9, November, 1960. 

I Assoc. Prof,, Mech. Engrg., Univ, of California, Berkeley, Calif. 

2 Research Asst., Inst. of Engrg. Research, Univ, of California, Berkeley, Calif. 

3 Research Asst., Inst. of Engrg. Research, Univ. of California, Berkeley, Calif. Pre- 
tly, Design Engr., Catepillar Tractor Co., Peoria, Ill, 

4 “Ocean Wave Forces onCircular Cylindrical Piles,” by W. L. Wiegel, K. E. Beebe, 

nd James Moon, Proceedings, ASCE, Vol, 83, April, 1951, 

5 “waterForces on Accelerated Cylinders,” by A. D. K. Laird, C, A. Johnson, and R, 


alker, Proceedings, ASCE, Vol. 85, No, WW 1, March, 1959. 
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TEST CYLINDER WITHOUT NEIGHBORS,# 


TABLE 1,—RATIOS OF AVERAGED MAXIMUM DRAG FORCES ON TEST CYLINDER 
WITH PARALLEL NEIGHBORING CYLINDERS, TO MAXIMUM DRAG FORCE ON 


to left and right of test cylinder 


b 


(628 runs used).9»° 


(A) Two neighbors symmetrically spaced 


To each side, diameter 


Ahead 
diameter 


0.50 0.75 1.00 1.25 1.50 2.00 2.50 3.00 3.50 4.00 4.50 5.00 


AE HEFDADHROHMNOTRHAAHODNMOE NDE NOWNH oOo 
NHN HAHCOOSCHOSOHRHSCORDHRARORRWROOH 
e @G &, Bae See oe so. «¢ * B.S BP a Se eee 8 ee ce eee 
HHH AHHH HA HAA HOH AHA OHOOCOCOHOONHNHS 
KE ODOONHOEFANAEANDOALFADANYWDWDOAMYROON 
NHN AAANHDODCORGOSCDSCSORSOSSCSARQAAN 
LWA AH AAAS SSCORSSSSSROGGSOARARR 
HHH HAHA AHA AA HA OCOn AAA AON NAH OOSCSO 
OHOMANHHANMYHERANHADDWDOBAMWNMNHAMOSOH 
NNAHOCCCOCOCOHORGOMROCOCSCSOReooOeoCeS 
LA ASSSSSSANORSSARGSGS 2S AQISSS 
HHH HFA TH AHH HAHAH HAHOKR HON TAA nH HAHOn NAA eA 
NMNOMANCOMAHHAHADWDANRHONHADAMNMNAAeL WO 
RAAQAAAASSSOAASSASHASSSREARS 
ee ee ee ee ee ie ee 
HAH OCT ANNAMAMNDARHWOErNOAMDOONKEHMO Oe 
ARADASASSSSSSSSSRASGARGRONN AY 
HHH HAHAH AA nA dn AAA HOCH OH OHnHOOnrFN 
ANRDOONNAMEEHADHORBOYAHMONOMN HHH Oo 
AA AASLNSSSSRAARSSCSCORGSRARGGAS 
BHA nn An AAA AA OCOCOH On FA nA HAOCOnNnNOONnAeA 
AHADMDDADNHAHNAHAHANNANDADONOHDOMHWABWANA 
AAASSRASSARAARARRORARORARARS 
HA ne ACOH HH OHOCOCCCCOHnHONnOHOCOOCOSOCoOoSorT 
Ke HDMNADANHRADHMNNONMWOFHALOCWOrHADA 
ARASSAQSRAAARIRGSARSARRRGRASG 
ACOH HOCH OH HOH OHH OOCONHoCOCOOCOoONnOon 
FMODOOMNAHWARANCODNOHAMODOAWOWAWDWWM MO 
WARAALAARAACS HOO RARHAAHLADAARARR 
BMH OCOCHOOCOCOHHHOHOOCOSCOOCOSCOCSCSCoCooCoCoSoS® 
OE EFArFDHDOROMRBRODNYHTEFNAAGAHr- OMWMD 
FOF EF EHACHAHARDSCrHOANDSOAMAAWDAAWBDAMD 
B). See ed? Che es Be © © © © #0F @ (8 w (6) Ce chee 68) 8) eee eee 
mMHOoCCOOnH HH OCOOCOHOHOOHOOCOSCSCCCoCSoCSoCSS 
WDOSTCOMOONHOANHNANNWNEANANTHNOHHOWOKOhe Or 
Hr HARABMRDHOMDMAANMWMAOCAMAAAIWNDAHADAWAANH 
LR AQAARACAARAR ,AAGLARAARAAARAGHAAR SG 
cocooooonnoocooonocoocoocoooooeo“ecocy([ 
ADDN OFNDHOMONANAMNANDONMWNMOAH WH 
Hod HD © OH WD ARRRAODARADA AD WD QW A  O 


G 


SoscsccSosccoSO coco So coH SSS COSCO SOSSCS 


° 
SOHHANNDHDHHMODOOKEDODBROSCHHANAND 
Aron nintnte 


(B) One neighbor (179 runs used) ,>»¢ 


To side, diameter 


Ahead 


diameter 


3.0 


2.5 


0.5 275 10 1,25 1.50 2.0 


0.0 


MNUANSOrwH oO ONDHA HS 
HooonHoOSCoOSCHHAHAAde 


enononononaononoak 


SOHHANTATWKHKDOOrE 
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TABLE 1.—(CONT’D) 


(B) One neighbor (179 runs used),}»° 


Ahead To side, diameter 
diameter 
0.0 0.5 15 1.0 1.25 1,50 2.0 2.5 3.0 
8.0 0.96 
8.5 0.96 
9.0 0.96 
9.5 0.96 
OKO 1.04 
10.5 0.92 
11.0 0.98 
17.5 0.93 
E20 O97 
12-5 1.01 
13.0 0.98 


(C) Two neighbors symmetrically spaced to left and right of test cylinder path 
(64 runs used),°,° 


Ahead To each side, diameter 
diameter 
KE 2 4 6 8 
0 1.33 
1 0.98 1.09 pes 
Z 1G 117 1,02 fer) 1,16 
4 1,21 1,49 1.06 1,26 1,14 
6 PLT 1,27 0.97 1.08 1.14 
8 1,16 0.96 0.99 1.05 
10 0.98 1.02 1.04 1.00 
12 0.97 0.99 0.98 1,02 
14 0.90 1,14 
| 16 0.99 
18 1.04 


De 
(D) One neighbor (55 runs used) ° 


, Ahead To side, diameter 


diameter 
0 1 2 4 6 8 
OS ee eee 
7 0 1.45 1,17 ; 
x 1.16 1,10 
& 2 0.48 1.48 0.99 1,14 0.99 
4 0.80 1.18 1.01 1,10 1,10 
‘ 6 0.90 1,32 0.95 1,05 1.09 
: 8 0.89 0.92 0.99 0.98 1.08 
e110 1.04 1.06 0.99 
B 12 1,01 1,01 1,02 


_ 4 positions of neighbors are given relative to the axis of the test cylinder, in cylin- 
r diameters, Ratios were taken at the bottom of the first swing after release from rest. 
Cylinder diameter = 2 in, © Amplitude of oscillation = 18 in, ¢Cylinder diameter = 1.2 
€ Amplitude of oscillation = 6.4 in. 
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It was the object of this experiment to demonstrate that wave forces ona 
circular cylinder could be influenced significantly by parallel neighboring cyl- 
inders, and particularly by eddies shed by these neighbors. 


EXPERIMENTAL EQUIPMENT 


The horizontal components of the wave-induced water particle velocities 
past vertical cylinders were modeled by the motion of the cylinders performing 
simple harmonic motion in otherwise still water. The heavily weighted four- 
wire pendulum which was used to produce the motion is shown in Fig. 1 for the 
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FIG, 1.—APPARATUS DIMENSIONS 


1/2-in. and 2-in, cylinders. The test cylinder was mounted on a strain bar at 
the center of the pendulum frame. As listed in Table 1, the neighboring cyl- 
inders of the same diameter as the test cylinder were spaced, with respect to 
the test cylinder, from touching up to 8 diameters to the side and 18 diameters 
ahead of the test cylinder in the direction of the first forward swing after 7 
lease, 

For the 1/2-in. cylinder, the strain bar deflections were sensed by a ca. 3 
pacity type transducer, and recorded by an oscillograph. Since this strain bar 
was one-directional, drag and transverse, or lift, forces had to be record 
during different runs by turning the bar through 90°. The forces on the 2-in 
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ylinder were detected by strain gauge bridges and recorded by the oscillo- 
raph. Asimultaneous trace from a 2 gaccelerometer was recorded with each 
in of the one-directional bars for both the 1/2-in. and 2-in. test cylinders. 
bout half the runs were taken with a one-directional bar, Introduction of a 
vo-directional strain bar allowed simultaneous recording of drag and trans- 
erse forces, but the accelerometer trace was omitted. Calibrations of the 
rain bars were made before and after each set of runs by applying known 
ees at the waterline of the test cylinder by means of balanced strings, pul- 
‘ys, and weights. 


TABLE 2,—MAXIMUM FORCES ON 2-IN, TEST CYLINDER, IN GRAMS, REFERRED 
TO THE WATER LINE, 
—<————————————————— 


vo neighbors spacing in diam. 
I 1/2 ahead 2 ahead 3 1/72 ahead 


11/2 aside 2 aside 2 aside 
Drag Lift Drag Lift Drag Lift 
151 143 167 - 42 204 82 
159 - 68 219 -125 197 - 95 
166 122 175 91 159 89 
181 75 204 31 189 89 
182 157 212 - 84 204 94 
182 79 212 99 197 84 
189 - 88 219 -125 189 - 53 
182 99 219 120 212 84 
204 -131 233 - 68 241 - 79 
175 57 219 wert) LOT, -100 
204 -125 219 - 79 204 63 
182 - 73 acc ere 204 120 
204 31 aa erate 197 -125 
225 -136 AAS a5 204 47 
175 136 ot 5 A Py 60 a0 
189 94 oes tae Erene ° 


Average magnitudes 


185 101 208 86 200 86 
pe 
2 The two parallel 2-in. neighboring cylinders were spaced symmetrically to left and 
sht of the test cylinder path. Amplitude of oscillation, 18 in, Positive and negative 
; forces to opposite sides. Forces are maxima with test cylinder following neighbors, 
; not necessarily on first swing after release. 


To avoid the vibrations caused by suddenly applied accelerations, motion 
1g started from maximum deflection by relaxing the tension in a hand-held 
rd attached to the center of gravity of the pendulum. 

During many runs, a 16 mmcamera recorded the motion of the eddies shed 
the cylinders, the traces of the magnetic oscillographs, the pendulum posi- 
n, and a timer. Thus, visual impressions of the action could be checked by 
wing the film in slow motion, or frame by frame. 


ANALYSIS OF DATA 


Since the chief aim of this research was to discover effects of neighboring 
inders on a test cylinder, the forces onthe test pile by itself were made the 
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basis for comparison, The drag force histories for runs without neighbor 
were averaged from the oscillograph records andthe extremes were recorde 
as shown in Fig. 2(a). The curves were not corrected for velocity variations 
The drag difference curves shown in Fig. 3, resulting from subtracting th 
average force curve from the force curve for the run under consideration 
emphasized the effects of the neighboring cylinders, Fig. 3 shows typical tim 
position histories of test cylinders, test cylinder eddies, and eddies from tw 
neighboring cylinders showing comparison with the corresponding time fore 
difference curves (1/2-in. cylinders). 

The neighboring cylinders further decreased the velocity and amplitude o 
the pendulum motion over that caused by the test cylinder alone, For eac! 
number of cylinders, the decrement ofthe amplitude was measured so that th 
maximum velocity at the bottom ofthe first swing from release could be com: 
puted from the theory of damped simple harmonic motion. A multiplicativ 
correction factor was computed for each of the cases by dividing the square o 
the computed maximum velocity of the no-neighbor runs, by the square of th 
computed maximum velocity of the runs with different numbers of neighbors 
For one neighboring cylinder this factor was 1.03; for two it was 1.08. Each o 
the values in the tables has been multiplied by the appropriate correctio! 
factor, 

The lift forces proved to be non-replicative so they were not averaged 
Typical curves are included in Figs. 2(b), (c) and ), and maximum value: 
are listed in Table 2. From successive frames of the motion picture films 
measurements of the positions of the test cylinder, and of the eddies shed br 
the test cylinder and its neighbors, resulted in curves of cylinder and edd: 
distances from the release point of the test cylinder as functions of time, suc] 
as shown in Fig, 3 


RESULTS AND DISCUSSION 


It was found that the neighboring cylinders significantly influenced thi 
forces on the test cylinder as shownin Figs. 2, 3, 4, and 5 and Tables 1 and 2 
The curves shown in Fig. 4 were not corrected for velocity variations. Whe: 
the test cylinder was in the wake from its neighbors there was a general de. 
crease in the average drag force to minima of about 20% of single cylinde: 
values at mid-swing. When the test cylinder was not in a wake, there was a 
increase in mid-swing average drag to maxima of over 150% of single cyl: 
inder values. Superimposed on this general increase or decrease of averagi 
drag were fluctuating drag and lift forces caused by eddies passing the tes 
cylinder. The largest fluctuations of mid-swing lift and drag forces were ove 
90% and 150%, respectively, of the maximum drag force on the test cylinde 
without neighbors, Since the test cylinder experiences increased average dra 
forces when outside the wake of a neighboring cylinder and maximum lift an 
fluctuations of drag forces only when it lies within the wake, there is littl 
possibility of complete reinforcement of these effects. Nevertheless, reinforce. 
ment did occur, as shown by the highest and lowest recorded drag forces-a 
mid-swing of 1.75 and minus 0.17 times the average mid-swing value for sing] 
cylinders, q 

The average resistance forces on a single 2-in. cylinder are shown in Fig 
2(a) and the envelope enclosing all the corresponding resistance forces me 
ured is shown in Figs. 2(a) and A. The drag coefficients corrected to inf 
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cylinder values were the usual 1.3 to 1.4 for vertical cylinders extendin; 
through the water surface. The corresponding lift or side force envelope i 
shown in Fig. 5. Data for the 1/2-in. cylinders were the same but for scale 
Typical resistance and lift forces on the test cylinder in the company of ; 
pair of neighbors placed as in Fig. 6 make up Figs. 2(b), (c) and (ad). In these 
significant departures of the resistance forces from those for the cylinde 
alone are evident. The closer the spacing, the greater and the earlier are th 
departures. 


FIG, 6.—CYLINDER AND EDDY MOTION 


In Figs. 4 and 5 the principal effects of a single 2-in. neighbor direct 
ahead of the 2-in test cylinder can be seen as departures from the ranges ¢ 
resistance and lift forces found without neighbors. Again, the deviations fro} 
single cylinder behavior decrease as the distance between the cylinders ir 
creases, and become unimportant at 4 diameters. 

Typical differences between the resistance forces on the 1/2-in test a 
inder without neighbors and with two 1/2-in. neighbors are shown in Fig, 
with typical test cylinder and eddy loci corresponding to these force trace; 
The eddies from the neighbors caused the fluctuations in the force trace; 
Because the times of shedding of the eddies from the test cylinder do not cor 
respond to the force trace deviations, these deviations must be attributed | 
the approach of the eddies from the neighboring cylinders, rather than to th 
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shedding of eddies from the test cylinder. It is evident that the closest ap- 
roach of the first shed outboard eddy pair to the test cylinder coincides with 
he same point on the first dip of each of the resistance difference curves. 
These eddies have a direction of rotation that causes a decrease in the veloc- 
ty of the test cylinder relative to the surrounding water and must therefore be 
he reason for the negative force contribution. Subsequent outboard eddies 
‘ause Similar negative contributions. For the run at 8 diameter spacing, two 
airs of inboard eddies could be seen onthe film. Their loci precede the cor- 
‘esponding outboard pairs and they must, because of their direction of rotation, 
vause the increase in drag contributions shown before each decrease. Un- 
ortunately, on the films the inboard eddies were very difficult to follow be- 
vause of the surface disturbances behind the cylinders. As a result, none were 
slearly defined for the 4 diameter and 12 diameter spacing runs of Fig. 3, One 
nay, however, infer their existence because inallof the runs in which fluctua- 
ions of drag force were noticeable the first deviation was an increase, All 
lata support the facts that the shedding of eddies from the test cylinders when 
lone or with neighbors, caused no significant change in the resistance force, 
ind that large changes were caused when the eddies from the neighboring cyl- 
nders passed the test cylinders. 

The maximum recorded lift force of 157 gm is shown in Table 2, This force 
s some 80% of the average single cylinder drag force and about triple the 
erage single cylinder lift force. Lift forces of 50% of the single cylinder 
lrag force were common. Drag forces were often increased by 50%. Thus, the 
ift and the drag force acting simultaneously may augment the single cylinder 
orce by at least 60% and cause oscillating forces on piling. 

Often the motion pictures showed that an eddy was shed almost simultane- 
usly from each side of the test cylinder without causing a lift force. In many 
tases, the eddy that formed first was not shed first. It sometimes shrank while 
he opposite eddy grew and was shed. These starts were often followed by ulti- 
nate shedding as in the single cylinder von Karman street. When neighbors 
vere close, however, the test cylinder sometimes shed only one eddy, and 
ometimes none, as shown in Fig. 3. The frequency of shedding shedding of 
‘ddies from the neighboring cylinders can be judged, from Fig. 3, to agree 
vith that given by Landweber’. The velocities of these shed eddies were 
neasured in this and other figures for 1/2-in. and 2-in. cylinders and were 
ound to agree with Landweber’s values for eddies shed by cylinders traveling 
it the mid-swing velocity. 

The lift force fluctuation periods agreed more closely with the single cyl- 
nder eddy pair shedding period given by the von Karman vortex street anal- 
isis than with the periods for eddy pair shedding from the neighbors as given 
lsewhere for medium to wide spacing of the cylinders. With these spacings, 
Here was little measured difference between the lift force records and the 
ingle cylinder lift force records. It is therefore indicated that the lift forces 
ere caused mostly by the shedding of eddies from the test cylinder with but 
tle effect from neighboring cylinders, unless the neighbors were within 4 


6 “Flow about aPair of Adjacent, Parallel Cylinders, Normal to a Stream: Theoreti- 
1 Analysis, U. S. Navy, The David W. Taylor Model Basin,” by L. Landweber, Report 
, July, 1942, 
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diameters ahead and 3 diameters to the side. The directions in which the lif 
forces were exerted were not replicative, as shown in Table 2. 


CONCLUSIONS 


1. The presence of neighboring cylinders may increase the drag force of a 
cylinder by 75% or decrease it to less than minus 15%. 

2. Lift forces of 75% of the single cylinder drag force may result from the 
presence of neighboring cylinders. 

3. The average drag force ona cylinder may be increased by as such at 
45% if outside, or decreased by 80% if inside, the wake of close neighbors 

4, Fluctuating drag forces are caused mostly by eddies from neighborin; 
cylinders. 

5. One neighbor has approximately the same effect as two. 

6. Neighbors more than 10 diameters from a test cylinder cause onk 
minor drag and lift disturbances. 


. 
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FLUME STUDIES OF FLOW IN STEEP, ROUGH CHANNELS* 


By Dean F. Peterson,! F. ASCE and P. K. Mohanty” 


SYNOPSIS 


Based on observations of both natural streams and flume experiments, 
wee classifications of flow in steep, rough, open channels are proposed: 
‘anquil, tumbling, and rapid. The tranquil and rapid regimes are dominated 
y subcritical and supercritical flow, respectively. In the tumbling regime 
1e roughness elements induce supercritical flow over their crests followed 
y individual hydraulic jumps to the tranquil state. The characteristics of 
ow in each of the proposed regimes are described, and criteria for deline- 
ting the regimes are proposed. 


INTRODUCTION 


Uniform flow is commonly considered to occur when the water surface is 
irallel to the bed slope. For such flow the depth, water area, velocity, and 
ischarge are the same at each section of a channel reach. Whereas these 
mditions may be essentially correct for open channels in which roughness 
lements are small, relative to depth, they are hardly correct for channels 
ving large roughness elements such as occur in mountain streams or in al- 
vial channels in which large dunes are formed. Existing formulas relating 
e gross variables of discharge, water area, depth, bed slope, and roughness 
at are based onthe concept of uniform flow, could hardly be expected to apply 
the latter case, especially if regions of supercritical flow form over the 


a 
Note,—Discussion open until April 1, 1961. To extend the closing date one month, a 
titten request must be filed with the Executive Secretary, ASCE, This paper is part 
the copyrighted Journal of the Hydraulics Division, Proceedings of the American So- 
ety of Civil Engineers, Vol. 86, No, HY 9, November, 1960. 
4 Presented to the Hydraulics Division, ASCE, July 1, 1959, Fort Collins, Colo. 
Prof, of Civ. Engrg. and Dean, Utah State Univ., Logan, Utha. 
Research Assoc., College of Engrg., Utah State Univ., Logan, Utah. 
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tops of the roughnesses. For this type of flow in nature, the basic variables 
can probably no longer be described by simple numbers, but become statistical 

The purpose of this study was to gain some insight into the basic dynamic 
relationships of flow in open channels, especially in steep channels, in whicl 
roughness is an appreciable portion of depth. Since conditions in nature art 
extremely complex, laboratory flume studies under idealized conditions 0 
geometry were made as a first step in trying to understand the problem mort 
fully. 

The design of the experimental program and the apparatus was based, amon 
other things, on consideration of many informal observations in the field, es: 
pecially of Logan River located in northern Utah. This river is a typical, steey 
(50 to 75 ft per mile slope in its lower reaches), boulder-strewn mountail 
stream with discharge normally ranging annually from a minima of aroun 
100 cfs to a maxima approaching 2000 cfs. The writers were impressed by thi 
fact that the roughnesses commonly appeared to induce a “tumbling” phenom 
ena whereby localized areas of supercritical flow followed by hydraulic jump: 
occurred. Fig. 1 shows this condition on the Logan River, that is typical o 
rough, steep streams observed at many other locations. 

It seemed desirable that the flume be capable of rather large slopes, an 
that the roughness elements should be such as to induce large variations i 
depth and velocity. Since these latter effects were of primary importance, an 
because of a desire to separate them from other considerations and make ther 
as simple as possible, two-dimensional roughness elements seemed appropri 
ate. A second step could be to introduce three-dimensional roughnesses, an 
a third, to investigate similitude by using similar roughness elements of dif 
ferent size. As a fourth step, consideration was given to the introduction ¢ 
roughness elements having a distribution of sizes, however, limitations 0 
time precluded this. 

A number of preliminary theoreticalinvestigations were made utilizing di 
mensional and mechanical considerations. These indicated that some form ¢ 
Froude number, together with slope, some description of roughness heigh 
spacing, and intensity should be the important variables. These analyses coul 
be developed to some degree for specialized ideal cases in the flume, ht 
yielded little information regarding the general case. 

This paper summarizes the principal results and conclusions reached t 
the study. More detailed information, including a complete tabulation of dat 
are included in Mr. Mohanty’s dissertation.? 4 

Notation.—The letter symbols adopted for use in this paper are define 
where they first appear, in the illustrations or in the text, and are arrange 
alphabetically, for convenience of reference, in Appendix II. | 


FLUME EXPERIMENTS : 

43 

‘ 

Flume.—Opportunities at hand dictated the design of the water suppl 
flume, and measuring devices. The installation was made in the out-of-doo! 
at the site of the Utah Water Research Laboratory (Fig. 2). Water was-s 


phoned from the reservoir above the 40-ft high State Power Dam to a 7-ft-b’ 
3 


3 “The Dynamics of Turbulent Flow in Steep, Rough, Open Channels,” by P. K, u 


hanty, thesis presented to Utah State Univ., at Logan, in 1959, in partial fulfilment 
the requirements for the degree of Doctor of Philosophy. = 
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5-ft-by-7-ft headbox through a 10 in. pipe controlled by a gate valve at th 
headbox. A maximum discharge of 4 cfs was available. 

The flume was 2 ft wide, 16 in. deep, and 64 ft long. It was constructed « 
plywood on longitudinal timber girders supported by a hinge at the center ar 
by two differential chain hoists 9 ft from each end. It could be set on any slog 
between 0 and 8.5%. At the headbox connection, a sliding arrangement consist 
ing of a non-stretchable piece of nylon cloth 16 ft-by-14 ft with a hinged ar 
sliding door proved leak-proof and otherwise satisfactory. A sliding tail gat 
was provided at the discharge end of the flume. 

The discharge from the flume spilled into a stilling pool excavated in tk 
earth but lined with concrete. At the end of the pool, a V-notch weir was ir 
stalled. Depth over the weir was measured at a Stilling well with a micrc 
hook-gage. A floating slab of plywood served as a wave suppressor. As 
check, a commercial flow-meter was installed in the siphon pipe. The wei 
measurements, however, were uSed in the analysis. 

Two-dimensional roughness elements consisted of 0.3 ft-by-0.3 ft or 0.1% 
ft-by-0.135 ft square timber bars (or a combination of them). Three 
dimensional elements were 0.3 ft-by-0.3 ft-by-0.3 ft timber cubes. Thes 
were machine planed and fastened to the flume bottom by toenailing. A ten 
plate was used during installation to assure proper spacing. Roughness ele 
ments having finite width were chosen because of ease of attachment and be 
cause some closer resemblance to natural roughness, having finite widt 
could possibly be implied. 

Measurements.— The slope of the flume was carefully checked at 2 ft st 
tions using an engineers level. Depth measurements were made by means 
a portable micro-point-gage. Measurement of the depth Y, over the leadil 
edge of the roughness was repeated at from 5 to 10 roughness stations in 0: 
der to check uniformity of the cyclic pattern. The average value was used 
analysis. Depth was also measured along the longitudinal center line in ord 
to obtain the flow profile. A large number of closely spaced piezometer hol 
were provided along the longitudinal center line. These were connected to 
manometer board and recorded localized piezometric head at the bounda: 
that usually differed appreciably from the depth of flow measured at the san 
station. 

Discharge was measured using the 90° V-notch weir as described previous] 
and the temperature of the water was noted for each run. 

Procedure.—After the roughness elements had been installed and the slo 
fixed, the gate valve was opened enough to allow the desired flow. After t 
flow was established, discharge measurements were made. The rubber tub 
leading tothe manometer board were drained to remove entrapped air. Man 
meter readings and observed values of depth were recorded after allom 
sufficient time for attainment of equilibrium. 

Using a reflex camera, at least two photographs of the flow pattern we 
taken for each run. A journal describing the physical characteristics of fk 
and noting peculiarities, if any, was kept for each run. The following sho1 


the range of the variables: 4 
Variable Symbol Range 3 
Discharge Q 0 to 4 cfs 3 
Slope of channel s 0.01 
0.0246 
0.05 


0.083 
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Roughness height K 0.3 ft-by-0.3 ft bars 
0.135 ft-by-0.135 ft bars 
0.3 ft-by-0.3 ft-by-0.3 ft cubes 


Longitudinal spacing if 0K 
of roughness 10K 
15K 


total of 163 runs were made. 

Dimensional Analysis.—A number of possibilities for measuring and char- 
terizing velocity and depth for the flume studies are apparent. Several of 
ese were tried during the investigation. However, if a cyclically uniform 
eady state is established, a particular velocity and depth should serve as an 
dex to which all other elements of the velocity-depth cycle may be referred. 
yr the purpose at hand, the depth Y, and velocity at the front edge of the 
ughness element proved most satisfactory (Fig. 3). Letting the correspond- 
g velocity V; serve as the dependent variable, 


Vp =9,8pKbBLtS5 8, n, A} esis slappiaas (1) 


which K, b, L, and t are, respectively, the height, top width, longitudinal 
d transverse spacings center to center of the roughness elements, S is the 
annel slope, and g and v are the acceleration of gravity and kinematic vis- 
sity of water, respectively. The factor \ represents a roughness intensity 
ctor. Since roughness elements with sharp edges were used, fixed lines of 
paration were anticipated, and thus the form drag, and consequently, Vy 
ould be nearly independent of Reynolds number. Choosing g and Yq as re- 
ating variables and neglecting y and b, Eq. 1 may be reduced to 


V,/Ve ¥ = Fy = 69(K/¥p L/¥y, t/¥, 8,2)... (2) 
1 roughness intensity factor \ may be expressed as 
Sho KA) pt) (Eee LTR Ree 4 (3) 


which p is the number of transverse roughness spacings per unit width of 
e channel, and ais a shape factor. For two-dimensional roughness, p t = 0 
d =a K/L. For the two-dimensional roughnesses, 


Fi = $3(K/Yz L/Yz S, PRR ERD. RO (4) 
Since the roughness elements are of standard shape, a was not a variable 
r these experiments. 


é PRESENTATION OF RESULTS 


Classification of Flow.—Early observations of natural streams focused the 
tention of the writers on a condition of flow dominated by scattered regions 
alternate acceleration and deceleration through critical flow over the large 
ughness elements. This condition was termed tumbling flow because the 
iter appeared to tumble over the roughnesses. At lower discharges, in 
aches of lower slope or smaller roughness, however, the standing waves 
a with the hydraulic jumps were no longer in evidence, and the flow 


. 
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appeared to be in atranquil state even over the tops of the roughness element 
At higher stages, the subcritical depth pools below the large roughnesses we 
swept out by the increased momentum, and the flow appeared to skim over t 
roughnesses at supercritical velocity. These observations were repeated 
the flume experiments, and the flow was, accordingly, classified into thr 
major regimes; (1) tranquil flow, which, essentially was maintained throug 
out, (2) tumbling flow—supercritical flow occurred over the roughness el 
ments with hydraulic jumps and subcritical flow in the regions between, a 
(3) rapid flow—supercritical flow occurred generally throughout the stream 

In the flume, the flow could be readily classified into the principal regir 
by visual inspection. The depth Y; also helped distinguish the rapid and tra 
quil flow regimes because Y, was smaller than Y¢ in the rapid flow regir 
and appreciably greater than Yq in the tranquil flow regime. No hydrau 
jump between roughness elements formed for either the rapid or the tranq} 
flow regimes. 

In planning the program, the writers were quite preoccupied with the tul 
bling flow phenomenon and, accordingly, most of the experimental runs we 
consciously chosen to assure that tumbling flow was induced. Nevertheles 
many runs occurred in the tranquil and rapid regimes. In presenting the r 
sults, the latter are presented first, even though they are not as extensive 
for the tumbling problem. 

Tranquil and Rapid Flow Regimes.—For flows in either of these regime 
successful correlations in the form of Eq. 4 were possible. Except for t 
basic difference of subcritical and supercritical velocities, the correlatio 
obtained were similar for both regimes. Fig. 4 shows a plot of K/ Y; vers 
F, for tranquil flow over two-dimensional roughness elements. For this p 
all roughness elements were 0.3 ft-by-0.3 ft except for those points repr 
sented by the triangles with the right half blacked, (S = 0.01, L = 15 K) tl 
were 0.135 ft-by-0.135 ft. The curves in Fig. 4 are typical in shape to tho 
for various slopes in both the rapid and tranquil flow regimes and for both 1 
two-and three-dimensional roughness elements. 

Study of Fig. 4 indicates that for intermediate values of roughness spall 
the Froude number is smallest, and that for low values of relative roughne 
the Froude number tends to increase. These observations are consistent w 
those made by Morris¢ in connection with pipe flow. Morris proposed 1 
following classifications depending on roughness spacing: 


1. Isolated roughness flow.—The roughness elements are sufficiently aif 
so that the wake and vortex generating zone at each element is comply 
dissipated before the next zone is reached. 

2. Wake interference flow.—The roughness elements are close enougt 
interfere. 

3. Quasi-smooth or Skimming Flow.—The roughness elements are aa 
ciently close that the flow essentially skims the crests of the elements. 


¥ 


The decrease in Froude number for the intermediate values of roughné 
spacing could be associated with wake interference, whereas, the increz 
associated with decreased height and spacing of roughness could indice 
skimming. Similitude appears to existas is indicated by the curves for S = 0 


4 «Flow in Rough Conduits,” by Henry M, Morris, Jr., Transactions, ASCE, Vol. : : 
1955, pp. 373-398. 


y 9 FLUME STUDIES 61 


0.3'x03' ROUGHNESS 
@ L= 5K, S=0.0246 
© L= 5K 
® L=I0KS=0.010 
© L=I5K 


Del ia SES PEs 
Polbeedieeanicearerdim 
ca ff tt 
pasbakababels tal [Ir 


FIG, 4,—= VERSUS F, IN TRANQUIL FLOW 
il 


0.4 
L=5K S=0.050 


@ 0.30 ftx0.30 ft bar roughness 
© 0.135 ftx 0.135 ft bar roughness 


0.3 


Value of Y, 
° 
iy 


0.1 


yee of Y, 


FIG. 5.—-SIMILITUDE IN RAPID FLOW 


«4? fie 


62 November, 1960 HY 


and L/K = 15 in Figs. 4 and 5 that are drawn using similar roughness config 
urations of different size. 

Fig. 6 shows occurrence of tranquil, quasi-smooth flow in a fairly roug! 
steep, natural stream. The roughness elements are of quite high intensity ar 
their heights are estimated to be of the order of 0.15 Yy. The velocity i 
probably approaching, but is less than F, =1. Fig. 7 shows natural rapid flov 
also probably in the skimming region. The flow of Fig. 7 differs from that « 
Fig. 1 inthat essentially it is completely supercritical; no hydraulic jumps ar 
able toform between the roughness elements, and the magnitude and directic 
of velocity remains essentially the same over the roughness elements and be 
tween them. Fig. 8 shows tranquil and rapid flow in the flume. 

Tumbling Flow.—The flow in this regime was characterized by alternat 
acceleration and deceleration from subcritical to supercritical velocities, an 
vice-versa, in a cyclic order. The spacing between two successive rough 
nesses marked the length of one cycle. 

Hydraulic jumps formed between successive roughness elements. Fig. 
shows typical flow patterns in this regime for the two- and the threes 
dimensional roughness elements installed in the flume. Fig. 1 illustrates 
corresponding flow pattern in nature. 

Based on considerations from dimensional analysis, a number of correla 
tions were tried. Whereas some appeared to be successful, the simple rela 


tionship between critical depth Xo Va2/e and Yj was most satisfactory} 
The flow characteristics were independent of K/L and S. 

As in Figs. 1 and 9, it can be seen that each roughness element acts as a 
overfall control. In the two dimensional studies, the control extends entire] 
across the channel. For the three-dimensional arrangements, the extent ¢ 
the control depends on the configuration of the roughness pattern and the lat 
eral spacing of the roughness elements. In the latter case, the stream tend 
to flow at subcritical velocities between the roughness elements and at super 
critical velocities over them. The flow, however, is carried by momentur 
both over and around the roughness elements and consequent acceleratio 
causes a region of supercritical flow much larger than that actually occupie 
by the roughness (Fig. 9). A region of supercritical flow induced by a solitar 
rock in Logan River is shown in Fig.10. The rock is only about one-third th 
width of the area disturbed. One would expect that a relatively low intensit 
of three-dimensional roughness elements would induce the equivalent of fw 
overfall control in a rough reach of stream under tumbling flow. Thus, th 
discharge under conditions of fully developed tumbling flow would be inde 
pendent of K, Land S. For this case, the depth Y,; over the tops of the rough 
ness elements would be less than critical depth computed on the basis of aver 


, 3 : 
age q, that is, Y, = V/ q2/, =/¥,+c> Forvthe tumbling regime, therefore, th 


formula for the critical flow weir, q = CVeg cs should form the basis for th 
discharge-depth relationship in which C and c are constants depending o 
roughness geometry and intensity; that is, ie 3 


In Fig. 11 is a plot showing Yo= Va2/e versus Y, for the 0.3 ft-by-0.3 
bar roughnesses. The plots for the 0.135 ft-by-0.135 ft bar roughness (n 
shown) fall on exactly the same curve. From this figure C = 0.83 and c =| 
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FIG, 6.—NATURAL TRANQUIL 
QUASI-SMOOTH FLOW 


FIG, 9.—TUMBLING FLOW IN FLUME 
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q= 0.83 Tito. .cgmnrs carrie ( 


Fig. 12 shows a similar plot for the 0.3 ft cube roughnesses. For this casi 
C =1.0 and c = 0.09, thus, 


therefore, 


The foregoing hypothesis permits one to substitute the concept of an ove1 
fallweir for a reach of tumbling flow. The single overfall for natural stream: 
instead of being located at a particular cross section, is now divided into mar 
smaller overfalls statistically distributed throughout the reach. Such a reac 
may be regarded as a statistical critical depth weir for which the discharg 
coefficients depend on the relative size, distribution, shape, and intensity « 
the roughness configuration. The depth-discharge relation depends on neithe 
the overfall height for a weir, nor on the length and slope of a tumbling flo 
reach of channel. 


FIG, 10,—EFFECT OF 
ISOLATED 
ROUGHNESS 


In nature, in which roughness elements are variable in size and intensity 
one still has the problem of deciding how to characterize the flow parameter 
of velocity and depth, that are now statistical. This foregoing concept of th 
virtual weir may furnish a means of approaching the problem. In such a weil 
the crest elevation might be represented by the proper statistical function de 
scribing roughness distribution for a tumbling reach of stream corrected fc 
separation and roughness intensity effects. Flow would be expected to occu 
at a minimum of total energy. In Fig. 13, 


WwW : 

a (z+y+v2/2e)y Viey dx. ee ( 

with the condition that % 
WwW = 

oe y WAV, dxdt eee, ee ( 


in which Z, y and Vy, are, respectively, the height of equivalent roughnes: 
depth of flow and velocity at a vertical section of the virtual weir, and y i 
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ie unit weight. P is the power and h is the elevation of the energy line. For 
minimum, dp/dh = 0. The writers were unsuccessful in using this approach 
mee V, is an unknown function of x. However, it now appears that statisti- 
ully within the reach qy = Q/w = Vx y dx may be a valid assumption. Using 
achine computations this would make evaluation of the integral possible. 

The limited data available need to be supplemented by flume studies of flow 
yer roughnesses having statistical distribution of size and spacing and by ex- 
msive study of natural streams in order to develop a generalized concept of 
imbling flow and to infer the best means for describing it quantitatively. 

Instability of Flow.—Roll waves occurred in the flume under certain con- 
itions; Fig. 14 shows their occurrence in the flume. The leading crest of the 
ave is at A. The trough is at B, building up to a crest at C. 

Roll waves have been reported in natural and artificial channels by a num- 
2r of writers, and several criteria have been proposed to predict the condi- 
ons under which they willappear. Depending on whether Chezy’s or Manning’s 
mula was used, the most frequently presented criterion is a lower limit of 


FIG. 14.—ROLL WAVES IN FLUME 


roude’s number of 1.5 or 2.0, respectively, implying that flow above this 
lue of Froude’s number would be inherently unstable, and that roll waves 
uld develop. Koloseus, attempting to connect the Darcy-Weisbach resis- 
nce coefficient with free-surface instability, concluded from experimental 
idies that for F < 1.6, theflow was stable and the resistance coefficient in- 
pendent of Froude’s number. When F > 1.6, the resistance coefficient was 
pendent on Froude’s number, the flow was unstable, and roll waves would 
rm if the channel were sufficiently long. ; 
In contrast with the results reported in the literature, the writers did not 
erve roll waves in rapid flow even though the Froude number frequently 
ceeded 1.6 and the flume should have been of adequate length. Rather, roll 


“The Effect of Free Surface Instability on Channel Resistance,” by H,. J. Koloseus, 
is presented to the State Univ., of Iowa, at Iowa City, in 1958, in partial fulfilment 
the requirements for the degree of Doctor of Philosophy. 
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waves formed in the transition ranges between tranquil and tumbling flow o 
between tumbling and rapid flow. As reported by others, they formed onl 
when the depth of flow was relatively shallow. No roll waves were observe 
using the cube roughness elements. Ordinarily, the large roughness element 
damped out instabilities, and a stable condition was quickly established. How 
ever, in the transition region a delicate condition of stability evidently exists 

During the present investigation, data on roll waves were collected onl 
incidentally and are, thus, too limited for systematic study. Two interestin 
observations were made: (1) For L/K = 5 roll waves occurred during th 
transition from tumbling to rapid flow only when F,/VS_ = 6, in which $i 
the channel slope, (2) For L/K = 10 or 15 roll waves occurred during tk 
transition from tranquil to tumbling flow when F; = 0.9. 

Air Entrainment.—The writers had anticipated that air entrainment mig] 
create a problem, especially with regard to measurements. Except for tk 
case of three-dimensional roughness under extreme conditions of slope ar 
discharge, however, air entrainment did not affect the measurements signifi 
cantly. For both tranquiland tumbling flow, the approachvelocity to the rougt 
ness crests was, of course, subcritical. 


ANALYSIS 


Flow in Steep, Rough, Channels as Uniform Flow.— Uniform flow, as corm 
monly defined, implies that depth does not change from section to sectio: 
This is a relative matter—truly uniform flow would be approached only in ver 
smoothchannels at low velocities, because any irregularities would cause som 
variation in depth, even though small. Where one draws the line is thus 
matter of judgment. For example, the flow shown in Fig. 6 would ordinari 
be classed as uniform. However, what about the flow in Figs. 1, 7, and 9 

In the flume studies under a regular roughness pattern, the water surfac 
varied in a regular repetitive cyclic way throughout the length of the flum 
once equilibrium was established. One might, therefore, regard such flow ¢ 
uniform, at least ina macroscopic sense. Over a reachof natural channel o1 
would suspect also that there would exist a statistical uniformity of flow, aga 
in a macroscopic sense. Perhaps a more general concept of uniform flo 
would be more appropriate based on equilibrium between the weight compone 
and the boundary drag (shear plus form drag of roughness elements) over 
sufficiently long reach of channel so that statistical description would be cot 
stant. Under this definition, except for the runs where roll waves forme 
macroscopically uniform flow was quickly established for all of the runs pei 
formed by the writers. 

Delineation of Flow Regimes.—The relation between discharge, slope, a1 
height and spacing of roughness elements determined the flow regime in tl 
flume. This is, undoubtedly, also true in nature, although the description 
the roughness may be quite complex. Since tumbling flow separates the re 
gimes of tranquil and rapid flow, an understanding of the conditions und 
which this can occur should lead to criteria delineating the principal regim 

The distinguishing characteristic of the tumbling regime is the occurren 
of stationary hydraulic jumps. Thus, an understanding of the conditions leadil 
to the formation of hydraulic jumps in sloping channels should be helpful. T! 
common equation derived from consideration of conservation of momentu 
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nd continuity between two successive cross sections, 
y d?/2 - y d2/2 + W Sin 6 = pra (gro Vy eves apes stte (10) 


3 difficult to solve because of the difficulty of determining the weight, W, of 
1¢ water prism. 

Hydraulic jumps in sloping channels have been studied by a number of re- 
earchers, and the results have been summarized by Bradley and Peterka.§ 
he writers found that these data may be summarized by the generalized em- 
irical formula, 

Mp oR gutta 40,253) /0.1 38 oon sn ae ceria binn (11) 


his is plotted in Fig. 15. 
For hydraulic jumps to be maintained, slope S must be sufficient to accel- 
rate all or part of the flow to supercritical velocities, and the height of the 
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FIG. 15,—-HYDRAULIC JUMP IN SLOPING CHANNELS 


yughnesses must be sufficient to induce the required sequent depth. The slope 
ust also be such that the form drag of the roughness elements and the shear 
1 the channel bottom may be overcome by the component of the weight paral- 
lto the bed. Form drag depends on the flow distribution around the rough- 
28s and, thus, on the flowregime. From rational considerations together with 
npirical evidence, tentative criteria delineating the various flow regimes in 
e flume are proposed (Table 1). Analysis leading to these criteria are in- 


uded in Appendix I. 


pn : : 
“Hydraulic Design of Stilling Basins: Stilling Basin with Sloping Apron (Basin V),” 
t nN. Bradley and A, J. Peterka, Proceedings ASCE, October, 1957, Vol. 83, No, 10, 


. 1-32, 
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Flow Formula.—Open channel flow formulas in common use may be e) 
pressed in the general form, 


V2.0 g ON" Biss tae eee oe (1 


in which m and n are constants and C’ is a function of the channel roughnes 
Although these formulas were originally proposed for tranquil flow and con 


TABLE 1.—CRITERIA FOR FLOW REGIMES 


Criterion Condition Type of Flow 


Two-dimensional bars 


(Fy ¥,/Y2)" ris 
tranquil or 
tumbling 

dg rapid 
tranquil or 
tumbling 

L Skimming 
wake interference 
or isolated roughnes 

L/(¥,-LS) tranquil 

K <i/2(¥0/Yg)-8/2(¥o/¥q) 

2 1| tranquil 


Three-dimensional cubes 


(1 Yy/Y2)” $(X1 + ¥,)/ (1.16 J ¥, + 2¥g) | rapid 


< (Yy+ Yq) / (1.16 j Yj + 2¥q) | tumbling or tranquil 


General 


tumbling will 
not occur 


tranquil 
tumbling or rapid — 


paratively low relative roughness, their use in the past has been extended 
include supercritical flow as well. 

In the tumbling regime, the depth of flow over the roughness crest tend 
pass through the critical value regardless of S, L, and K. Thus, in Eq. 12 
values of C', m, and a would tend to become 1, 1/ 2, and 0, respectively. ’ 
the other regimes Ce m, and n would approach those for a smooth c 
Within the transition zones, the values of C', m, andn apparently will vary d 
pending on proximity of the flow balance to tumbling flow. q 


i 
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Due to the extremely rough channel bottom as compared to its sides, the 
hannel should behave essentially as avery wide channel, and m should be very 
lose to 1/2. This strongly suggests the possibility that a general open channel 
ow formula may be expressed as 


1 which both C' and n are functions of relative roughness geometry only. 

The interdependence of Fi, K/Yj, and L/Yy, discussed previously, also 
uggests that the most important parameter representing roughness geometry 
ay be the roughness intensity A expressed possibly in the form 


MeL = Diy el oe te Ce eee (15) 


1 this form, A is the ratio of the projected effective area of the roughness 
srpendicular to the direction of flow to the total area of the channel bed. 

It is believed that both C' and n are functions of the relative roughness 
csometry, that is 


Cl= (05 B/%) K/%) 2. 2 (16) 


id 
n= Fylde LN g, BLYigia wasikes's macueee: (17) 


owever, lack of adequate data in allsub-regimes of flow prevents quantitative 
valuation of Eq. 16 and 17 for even the simplified conditions of the flume 
udies. 

Roughness in Hydraulic Design.—The writers believe that use of roughness 
lements have not been fully exploited in hydraulic design. With better under- 
anding, their rational use could very possibly result in improved perform- 
ice and economy of certain structures. Some specific possibilities might in- 
lude the following: 


1. Energy Dissipator.—A chute with a single large stilling basin is fre- 
ently very expensive. By proper use of roughness elements it might be 
ractical to establish tumbling flow throughout the entire length of the chute. 
he energy would then be dissipated by a number of overfall sinks throughout 
s length rather than by a Single large one. Velocities would never greatly 
<ceed critical. 

2. Use of Steeper Slopes in Channel Design.—In earthen channels steep 
opes resulting in scouring velocities are usually undesirable. In lined chan- 
Is of steep slope there is the danger of instability, impact waves, and forma- 
on of roll waves. Use of roughness elements to control the flow at a suitable 
locity might permit the use of steeper slopes where this is desirable. 

3. Increased Sediment Transport.—By using roughness to induce greater 
rbulence, the sediment carrying capacity of a designed channel may be in- 
eased. Benedict, Albertson, and Matejka’ used roughness elements in a 
ume todevelop sufficient turbulence to suspend the bed load ofa stream, thus 


7 «Total Sediment Load Measured in Turbulence Flume,” by P, C, Benedict, M. L. 
srtson, and D, Q. Matejka, Transactions, ASCE, Vol. 120, 1955, pp. 457-484, 
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making possible the measurement of total transport. This same idea could be 
utilized to insure transport of sediment loads through channels. 

4. Stream Measurement.—Through a thorough understanding of the role of 
extreme roughness elements in a naturalstream, an inexpensive and yet satis- 
factory method of estimating discharge of mountain streams may be developed 
eventually. Also, because uniform flow is rapidly induced by large roughness 
elements, a roughness flume might be developed for water measurement. This 
would seem to be a promising device for intermittent streams transporting 
large quantities of sediment. A properly designed tumbling-flow flume might 
permit use of a flume having approximately the same natural slope and rough- 
ness as the stream, thus avoiding the problems inherent in providing a stilling 
pooland sediment drop out, but still having the desirable characteristics of an 
overfall weir. 

5. Roughness Standard.—An artificial standard of roughness for flow in 
open channels comparable to that of Nikaradse for pipes is desirable, but has 
not been fully developed. 

6. Flow Analysis.—Understanding of the pressure distribution around 
roughnesses in open channels under various flow conditions would permit dy- 
namical analysis of many flow problems. 


CONCLUSIONS 


1. Flow in extremely rough steep open channels may be classified intc 
three regimes: tranquil, tumbling, and rapid. 

2. In the tranquil and rapid regimes, flow may possibly be described by ar 
equation of the type Fy = C'S", in which C' and n depend on the roughness in- 
tensity and configuration. 

3. In the tumbling regime, the hydraulic jumps form over the roughness 
elements to the degree that a reach of stream acts as a control. Sucha reac] 
may be regarded as analogous to a critical flow weir. 

4. For shallow depths under certain conditions, roll waves may form a 
the transition zone between tranquil and tumbling flow and rapid flow. 

5. Based on better understanding of their action, the use of large rough: 
ness elements could result in improved hydraulic design for energy dissipators 
improved sediment transport, and stream measurement. 
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APPENDIX I. CRITERIA FOR FLOW REGIMES 


Dynamic.—In Fig. 16(a) assume 
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From the observed data, this was true within 2%. Considering 1 ft of widtl 
and applying the principle of momentum, 


v(xi/2 4 v2/2) +WSin 6 - T+F,=pq(V2-Vj).----- (19 
in which F, is the pressure on the downstream face of the roughness. 
From Fig. 16(b) 
W Sin'@) = TeFp vastness cae (20 


and 
y Y2/2+pqVy+Fpt Fy = y ¥3/2+pa Vo iia (21 


For supercritical values of Vj, if the right side of Eq. 21 is greater than the 
left side, a hydraulic jump can form. Otherwise, the flow will tend to skim 
F, is part of the form drag on the roughness; thus, one can summarize 


Fe tp Cpl Kp. Vi/2 7: 2 eee (22 


in which Cp is a form of drag coefficient and j is the transverse width 0 
roughness per unit width of channel. Substituting and rearranging yields 


ae - 
F,(¥1/¥o)" = (¥y + ¥g)/ (Cpi¥y+2Y¥g).......-- (23 
For two-dimensional roughness, j = 1, using Ch = 2 arbitrarily, 


(Fy Uta): SE facccn we the eee (24 


For cube roughness using j = 0.3 Cp = 1.16, approximately, 
2 ; 
Fy (¥1/¥q) = (¥y + Yq) / (0.35 ¥y +2 YQ) -........ (25 


If (F, Y1/Yg)" exceeds the value of 1/2 or of (¥1 + ¥g)/ (0.35 ¥y +2 Yo) 
rapid flow would be expected; if not, either tumbling or tranquil flow woul 
occur. 

Sequent Depth.—For large L, Fig. 16(c), d would approach that compute: 
from Chezy’s formula for the appropriate bed, excluding the large roughnes: 
elements. If the required sequent depth dg exceeds K + Yq, rapid flow pre 
vails. That is if 3 


do = Fy dy (S.+ 0.153) /0,133 > KY, oss sca (26 


the flow will be rapid. 

Roughness Spacing.—In Fig. 16(d), 1 is the distance to the point at which th 
centerline of the jet intersects the line joining the tops of the roughness ele 
ments. If L <1, the flow will skim. Let x and z be the horizontal and vertica 
co-ordinates of the jet of velocity Vj at time t, then x = Vy t and Z = g t2/2 


For z = Y;/2 
x=1=5 Vi V¥i/e ee ee ee (2% 


Submergence.—Keutner® found that for Hy/H, > 0.85, Fig. 16(e), the nagl 
did not plunge, and the flow remained tranquil. Assuming Hy = Y, and Hg = Y 


8 “Neues Berchnungsverfahren fur den Akfluss anWehren aus der Geschwindegklit 
verteilung des Wassers uber der Wehrkone,” by C. Keutner, Die Bautechnik, p, 575, 192 
; 
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LS, Keutner’s criterion reduces to ¥1/(¥y - LS) > 0.85 would imply tranquil 
ow. 

Length of Jump.—If d is the uniform flow depth established fora channel 
«cluding the gross roughness elements, and do is replaced by Y, + K, Eq. 10 
2comes 


y d2/2 + pq V, +WSin 6 = y (¥, + K)*/2+pq Vo ere (28) 


q- 28 may be used to compute W. Assuming an average depth of jump to be 
4 + Y, + K)/2 (or any other appropriate average depth), compute the length 
jump 1j necessary to provide the weight. If this length exceeds the rough- 
SS Spacing, the flow cannot tumble. 

Residual Slope.—Part of the channel slope may be thought ofas providing 
‘avity force to overcome the drag offered by the roughness element. Let 
is be designated Sp. If S - Sp < Sc, in which S,; denotes critical slope, the 
slocity cannot reach critical, and tranquil flow will prevail. 

Roughness Height.— Let the velocity of approach, Fig. 16(c) be V,. Assum- 
g no energy loss, the equations of continuity and energy may be combined 
id rearranged to form 


i a [1/2 (¥,/¥,)° - 3/2 (Yo/Yq) + 1 pies ee (29) 


alues of K less thanthat given by Eq. 29 will be too small to induce tumbling 
OW. 


APPENDIX II. NOTATION 


The following letter symbols, adopted for use in the paper and for the guid- 
ce of discussers, conform essentially with American Standard Letter Sym- 
ls for Hydraulics (ASA-Z10.2), prepared by a committee of the American 
andards Association, with ASCE representation, and approved by the Associ- 
ion in 1942: 


b = top width of roughness; 

Sc = constants; 

Cy = drag coefficient; 

dy, do = sequent depth, hydraulic jump in a smooth channel; 
F, = V1/Vg Yz = Froude’s number at the upstream crest of roughness; 
Fy = drag force; 

g = acceleration due to gravity; 

Hy = upstream head over the weir; 

Ho = downstream head over the weir; 


= transverse width of roughness per unit width of channel; 
= height of roughness; 
= longitudinal roughness spacing; 
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= length of jump; 
= length of backwater curve; 
= exponent of depth; 
= exponent of slope; 


= number of transverse roughness spacings per unit widtl 
of channel; 


= discharge per unit width of channel; 


= channel slope; 


i 


critical slope; 


= slope used to provide gravity energy to overcome forn 
drag; 


= transverse roughness Spacing; 

= velocity at upstream roughness crest; 

= velocity of approach; 

= critical velocity; 

= weight of small prism of water in hydraulic jump; 
= width of flume or channel; 

= horizontal distance; 

= depth of flow above upstream crest of roughness; 
Schaal 

= depth of flow approaching the roughness; 

= vertical distance; 

= shape factor; 

= specific weight; 

= density; 

= symbols denoting functions; 

= shear; 
= angle of slope; 


= kinematic viscosity; 
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TESTS ON PRESTRESSED CONCRETE EMBEDDED CYLINDER PIPE 


By Hugh F. Kennison,! M. ASCE 


SYNOPSIS 


This paper presents the results of the first full-scale combined internal 
pressure and external load tests ever made on 60-in. prestressed concrete 
embedded cylinder pipe, designed according to AWWA Specification C301. The 
structural behavior of this pipe under combined loads is reported, based on the 
development of stress cracks in the concrete and on strain recordings from 
SR-4 strain gages. 


INTRODUCTION 


Prestressed concrete cylinder pipe has had world-wide acceptance since its 
introduction in 1942. In this original pipe, the prestressed wires are wrapped 
directly on the steel cylinder encasing the concrete core (Fig. 1). The pipe is 
manufactured in diameters from 16 in. to 48 in. In 1952 the first prestressed 
concrete embedded cylinder pipe was commercially produced. In this pipe, as 
distinguished from prestressed concrete cylinder pipe, the prestressed wires 
are wrapped about the concrete core in which the steel cylinder is embedded 
(Fig.2). This pipe has been manufactured in diameters from 24 in. to 120 in. 
and is particularly applicable for large diameter and high pressure pipelines. 
Both types of pipe are described in the American Water Works Association 
(AWWA) Specification C301-58. 

- Prestressed concrete embedded cylinder pipe consists of a steel cylinder 
with bell and spigot joint rings welded to its ends. This steel membrane is 
~ Note.—Discussion open until April 1, 1961. To extend the closing date one month, a 
written request must be filed with the Executive Secretary, ASCE, This paper is part 
of the copyrighted Journal of the Hydraulics Division, Proceedings of the American So- 


ciety of Civil Engineers, Vol. 86, No. HY 9, November, 1960. 
I vice Pres., Engrg. and Research, Lock Joint Pipe Co., East Orange, N. J. 
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FIG, 1.—PRESTRESSED CONCRETE CYLINDER PIPE 
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FIG, 2.—PRESTRESSED CONCRETE EMBEDDED CYLINDER PIPE 
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drostatically tested for water tightness and strength. The cylinder is then 
mcased in acore of quality concrete consolidated by high frequency vibration. 
\fter the core has been properly cured, it is helically wrapped with prestressing 
yire under ahigh uniform tension. The wire is then protected with a concrete 
yr mortar coating. 

To obtain realistic information onthe physical behavior of prestressed con- 
‘rete embedded cylinder pipe, the American Concrete Pressure Pipe Associa- 
ion (ACPPA) sponsored a series of tests at a laboratory in Wharton, N. J. 
This three-day test program was witnessed by consulting engineers and major 
sonsumers from New York, N. Y., Philadelphia, Pa., Baltimore, Md., Kansas 
sity, Mo., Dallas, Tex., Fort Worth, Tex., Ames, Ia., Denver, Col., Los Angeles, 
Talif., and San Francisco, Calif. These tests were designed to determine the 
structural behavior of 60-in. pipe when subjected to combinations of external 
oading and internal pressure, since under normal use these simultaneous con- 
litions would be in effect. 


PREVIOUS TESTS ON PRESTRESSED CONCRETE CYLINDER PIPE 


Earlier combined load tests onsmall diameter prestressed concrete (lined) 
xylinder pipe have indicated a cubic parabola relationship between three-edge 
yearing load and internal hydrostatic pressure, as shown in Fig. 3. 

The basic curves resulting from these earlier investigations are the “design” 
urve, and the “ultimate” curve. The “design” curve has been the principal 
ontrol for a combinedanalysis of prestressed concrete (lined) cylinder pipe. 
[he combination of normal operating pressure and dead load has been confined 
within a curve which is 90% of the “design” curve. 

The “ultimate” curve represents the combinations of pressure and load which 
vill cause failure. This curve is by far the least realistic of the curves and 
‘onsequently seldom used. Under the three-edge type of loading this curve does 
1ot give an accurate representation of the pipe’s performance under field con- 
litions. As the pipe nearsits ultimate capacity it behaves as flexible pipe, and 
inder field conditions this capacity would be far greater than indicated, due to 
he extra side support given by passive earth pressures. In actuality, a much 
\igher bedding factor should be applied with the ultimate curve, to obtain real- 
stic results. 

The terminal ordinate and abscissa values of these basic curves can be cal- 
ulated and/or verified by relatively simple tests. These terminal values are: 
NQ.001 in, = three-edge bearing load producing the first 0.001-in. crack (in- 
-ipient cracking) ; Wut = three-edge bearing load producing ultimate failure 
f the pipe; Pp = pressure causing zero concrete stress in the core; and Ph 


: bursting pressure of the pipe. 


PIPE DESIGN FOR 60-IN. COMBINED LOAD TESTS 


~ Class 160 psi, 60-in. prestressed concrete embedded cylinder pipe was 
selected for this test program, since it falls approximately in the middle of 
he range of sizes and pressures for this type of construction. Seven identical 
est specimens were constructed and tested at various combinations of inter- 
al pressure and external loads. These test pipe were 4 ft long, 60-in. inter- 
al diameter and designed for a combination of a working pressure of 160 psi 
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id a three-edge bearing load of 6,175 lb per lin ft. A common design was 
iosen for all specimens to minimize the many variables which occur in such 
series of tests. The design, consistent with AWWA Specification C301, is as 
liows: 

Concrete core thickness = 44 in. (steel cylinder 14 in. from inside of core); 
merete coating thickness = 14 in.; cylinder gage = #16 U.S.S. hot rolled 
\STM A245); prestress wire = #6, MB untempered (ASTM A227); steel cyl- 
der area = 0.718 sq in. per lin ft; steel cylinder diameter = 63 in. (outside 
ameter); prestress wire spacing = 0.392 in., c. toc.; prestress wire area 
0.883 sq in. per lin ft; Py = normal operating pressure = 160 psi; Po = zero 
mpression pressure = 216 psi; P, = bursting pressure = 535 psi; W, = al- 
wable three-edge bearing load at operating pressure = 6,175 lb per lin ft; 
ult = load to cause three-edge bearing failure at zero pressure = 27,200 Ib 
r lin ft; gross wire wrapping stress = 140,000 psi; resultant wire stress 
105,000 psi; and wire stress at Py, = 111,570 psi. 

Fig. 4 shows the design curve, the 0.9 design curve, the ultimate curve, the 
yrmal operating pressure, and the normal three-edge bearing load for these 
st specimens. These curves will be reproduced on all subsequent plots of 
e test results for reference purposes. 


DESCRIPTION OF TESTING EQUIPMENT 


For these tests the ASTM standard three-edge bearing load was chosen as 
e method of applying the external load because it was felt that this loading 
uld be reproduced more consistently than a simulated trench condition. Ex- 
nsive research work previously performedat Iowa State College (Ames, Ia.) 
, Anson Marston, W. J. Schlick, and M. G. Spangler provides the necessary 


TABLE 1.—THREE-EDGE BEARING LOADS CONVERTED TO EQUIVALENT 
DEPTHS OF COVER 


Three-Edge Cover with Ordinary Bedding, in Feet 
searing Load, in mae = +2f | Bys=D,+ 3 ft 
d=D + 3ft | By=Do d =o 
Pounds per = 100 pef w = 120 pef w = 120 pef 
Linear Foot “4 as ; 


(1) 


x al, and 
This table of equivalent depths of cover is based on the studies of Marston, et al, 
sumed eae bedding (a bedding factor of 1.5 when compared with three-edge 
aring loads) and soil constants of KH and KH" of 0,130, the ordinary maximum for 
ay soils. 
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information for converting the three-edge bearing loads to various types of 
trench conditions. The use of this excellent material eliminates the additional 
variable of evaluating a simulated trench condition and allows a reasonable 
interpretation of relatively simple and routine crushing tests. Table 1 has 


Three edge 
bearing load 


Internal 
Hydrostatic 
Pressure 


FIG, 5.-SYSTEM OF FIG, 6,—60-IN. PRESTRESSED ase 
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LOADING DED CYLINDER PIPE 
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SEMBLED FOR COM- 
BINED LOAD TEST 


4'-0" 
Three-Edge Bearing Load 


60" Internal Diameter 


q 
FIG, 7.—SCHEMATIC VIEW OF PIPE SPECIMEN ASSEMBLED 4 
ON INTERNAL BULKHEAD . 


been prepared for converting subsequently recorded three-edge bearing te: 
results to the more realistic trench conditions, in feet of cover over the pipe; 
in Table 1 Bg is the trench width, D, denotes the outside diameter of the pipe 
and w is the specific weight of the soil. : 
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The combined conditions of three-edge bearing external load andinternal 
ydrostatic pressure under which these tests were carried out is illustrated 
n Fig. 5. The 60-in. test pipes, with 18-in. long thinsteel cylinder extensions, 
ere placed over a tubular internal bulkhead. The steel cylinder extensions 
educed the supporting effect of the bulkhead on the three-edge bearing load. 
‘he test specimen completely assembled in the combined load testing facility 
s shown in Fig.6. The arrangement of the test specimen mounted in the tubu- 
ar bulkhead is illustrated in Fig. 7. The support which was provided by the 
ydrostatic bulkhead was carefully determined andall three-edge bearing loads 


FIG, 8.—INTERIOR VIEW OF TUBU- FIG, 9.—EXAMINATION OF THE 
LAR BULKHEAD WITH CONCRETE COATING 
THICK GLASS WINDOW, 

SUBMARINE LIGHT HOUS- 
ING AND LONG FOCAL 
LENGTH MICROSCOPE 


2re corrected by this amount. The tests and results for evaluating this sup- 
ting effect are presented subsequently. 

Throughout the tests it was necessary toobserve the inside concrete lining 
the pipe to detect and measure the development of lining cracks. To ac- 
mplish this the tubular bulkhead was constructed with a thick glass window, 
rough which the inside lining of the pipe could be surveyed with a 40X micro- 
ope. This window was placed opposite the bottom of the pipe, which would 
the location of the maximum tensile stress in the concrete lining. A sub- 
arine lighting device was constructed within the bulkhead to illuminate the 
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inside bottom lining of the pipe over the area covered by the microscope. The 
arrangement of the glass window, the microscope, and the submarine light is 
shown in Fig. 8. 


PROCEDURE USED TO OBSERVE TEST RESULTS 


One of the major indications of adequacy of design in prestressed concrete, 
as in reinforced concrete, is the degree to which stress cracking of the con- 
crete is controlled under load. Reinforced concrete pipe must necessarily 
have tensile stresses and resulting cracking of the concrete whenever the ten- 
sile strength of the concrete is exceeded. Prestressed concrete pipe was de- 
veloped to prevent or reduce such stress cracking with resulting improvec 


Load 


Top Section Quarterpoint 


10 (a) 10(b) 


FIG, 10,.—LOC ATION OF SR-4 STRAIN GAGES 


performance. A considerable margin exists between the load or pressures 
causing first visible stress crack, and that causing ultimate failure. Betweel 
these two points, lies the critical stress crack condition which many author- 
ities consider to be in excess of 0.01 in. for conventionally reinforced pipe 
In prestressed concrete pipe, the critical stress crack occurs in the exterio’ 
coating and is conservatively defined as baying, a sustained maximum width i 
excess of 0.005 in. 

These tests report on the cracking of the lining and coating caused by th 
combined load application. The cracking behavior of the pipe was obtained b 
a careful examination of the concrete lining and coating with a 40X microscop 
equipped with a calibrated grid to accurately measure crack widths. The con: 
dition of “incipient” cracking or first tensile failure can be recorded by traine 
personnel through the use of sucha microscope. This microscopic crack i 
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uot readily visible to the naked eye and is defined as a crack, 0.001 in. in width 
and 12 in. long. The first visible crack is defined as acrack 0.002 in. wide and 
12 in. long and can usually be observed with normal eyesight if a meticulous 
sxamination of the surface.is made. The first visible crack supplies a consis- 
ent criterion of the stress cracking of the concrete and is so used in the re- 
cording and interpretation of the results of these tests. Fig. 9 shows the sur- 
vey of the pipe coating by microscope to detect and measure stress cracks 
9ecurring in the coating while the pipe is under a combined load. 

To further study the structural behavior of the pipe under combined loads, 
strain recordings were obtained through the use of SR-4 gages. These gages 
were installed in three test specimens at various locations as shown in Fig. 
10(a). During the tests the pipe were placed with the strain gages at specified 
ocations (Fig. 10(b)) relative to the points of application of the external load. 
strain readings were then taken during combined load applications. 


PRELIMINARY INFORMATION AND TESTS 


The seven test specimens were numbered 1834 through 1840 for purpose of 
-ecord. Table 2 gives the history of these seven test pipe. 


TABLE 2.—MANUFACTURING AND TEST DATES (IN 1955) AND CONCRETE 
CYLINDER STRENGTHS 


Date Tested Concrete Strengths, in 


Pounds per Square Inch 


Com- 


bination Core Coating 


7 Day |28 Day| 7 Day 28 Day 


(9) 


11/7 
9/9 & 9/13 


10/6 


‘Pipe #1836 was used in the test conducted to evaluate the effect of the hydro- 
tatic bulkhead in supporting external loads, as mentioned previously. In this 
“st the pipe alone was subjected to increasing three-edge bearing loads in 
500 lb per lin ft; increments up to 19,000 lb per lin ft. Vertical deflections 
ere measured to the nearest 0.001 in. at each increment of load. This test 
as repeated to record any effect of concrete tensile failure. The lower curve 
10wn in Fig. 11 resulted from these two tests. The 18-in. long steel extension 


linders were then welded to the ends of pipe #1836 and the unit assembled in 
hydrostatic bulkhead. The complete assembly was then subjected to the 
e-edge bearing test as before, the results of which appear as the upper 
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urve of Fig. 11. The difference in three-edge bearing load between the two 
urves, for a given deflection, is thentakenas the supporting effect of the hydro- 
tatic bulkheadat that load. Fig. 12 shows the resulting correction factor which 
vas used to obtain the true three-edge bearing load. 


TEST PROCEDURE AND RESULTS BASED ON CRACKING BEHAVIOR 


The following discussion contains a separate section on the tests and re- 
ulting data, based on cracking behavior, for each of the individual test speci- 
mens. Results obtained through the use of the SR-4 strain gages will be pre- 
ented separately following this section. A variety of test procedures was 
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FIG, 13.—COMBINED LOADS ON TEST PIPE NO, 1834 


arried out on each test pipe. These include: (1) subjecting the specimen to 
ombined loadings that trace the 0.9 design curve and the design curve; (2) 
ubjecting the specimen to a variety of combined loadings sufficient to produce 
incipient” cracking (0.001-in. x 12-in. crack); (3) subjecting the specimen 
» combined loadings sufficient to produce the first visible crack (0.002-in. x 
2-in. crack); (4) subjecting the specimento combined loads in which both the 
iternal pressure and three-edge load were increased by a predetermined 
mount; (5) subjecting the specimen to an increasing three-edge bearing load 
10 internal pressure) resulting in ultimate failure of the pipe; and (6) sub- 
‘cting the specimen to an increasing internal pressure (no three-edge bearing 
ad) resulting in ultimate bursting of the pipe. Each section will contain a 
sscription of the tests performedon the specimen along with a plot of the re- 
iting data on a combined load curve for this design such as was shown in 
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Figs. 13, 14, 15 and 16 show clearly the conservatism of the design methoc 
because of the very adequate margin between the theoretical design, or incip- 
ient crack, curve for the pipe tested and the actual pressure - load combina- 
tions which caused these cracks. 

Pipe #1834.—Pipe #1834 was first subjected to a series of combined loads 
that traced the 0.9 design curve. These combined loads are plotted on Fig. 13. 
It was next subjected toa series of combined loads that traced the design curve 
No microscopic cracking of the lining or coating was detected during these twe 
series of combined loads. Following this the internal pressure was held a 
160 psi (normal operating pressure) and the three-edge bearing load was in- 
creased until the first “incipient” crack (0.001-in. x 12-in.) formed. This firs 
“incipient” crack appeared at the top quarter point of the spigot end. The las 
test consisted of recording the cracking behavior of the pipe as the interna 
hydrostatic pressure (no three-edge bearing load) was increased to the burst- 
ing pressure. All loads applied to pipe #1834 are plotted in Fig. 13. 

Pipe #1835.—This test specimen was subjected to an increasing three-edg 
bearing load (no internal pressure) until it reached ultimate failure. The ver- 
tical and horizontal diameter changes were recorded along with the behavio: 
of the crack formation. The relationship between three-edge bearing load ant 
the diameter changes are shown in Fig. 17. The various cracking loads an 
cracking information is plotted on the ordinate of the combined load curve i 
Fig. 14. The results of the hydrostatic tests on test specimen #1839 and #184 
are also plotted in Fig. 14 as will be discussed subsequently. 

Pipe #1836.—This test specimen was used to develop the correction facto 
to be applied to the three-edge bearing machine load caused by the supportin 
effect of the hydrostatic bulkheads. This test was discussed and the result 
presented in Figs. 11 and 12. ’ 

Pipe #1837.—Thistest specimen was first subjected to various combination 
of internal pressure and external three-edge bearing loads which produced th 
first visible crack (0.002-in. x 12-in.) condition, The test procedure was t 
hold pre-established three-edge bearing loads and increase the internal pres 
sure until the first visible crack occurred or reoccurred under subsequent loads 
The combined loads resulting in the first visible crack are plotted on the com 
bined load curve in Fig. 15. In this series the first visible crack occurred i 
the pipe coating under all load combinations. The controlling crack under in 
ternal pressure and no external load occurred at the top of the pipe while th 
controlling crack under external load and no internal pressure was located : 
the springline (90° from the top of the pipe). : 

Next this pipe was subjected to three combinations of internal pressure ar 
three-edge bearing loads in which these loads were increased simultaneous! 
to 2.0, 2.5, and 2.75 times the normal operating loads of 160 psi internal pres 
sure and 6175 lb per lin ft three-edge bearing loads. The maximum crac 
width at each of these three combination loads was recorded. This informs 
tion is recorded on the combined load curves of Fig. 15. , 

Pipe #1838.—As in pipe #1837 this specimen was first subjected to variot 
combinations of internal pressure and three-edge bearing loads which produce 
the first visible crack (0.002-in. x 12-in.) condition. The procedure on th 
specimen was reversed so as to hold a predetermined internal pressure 11 
increase the three-edge bearing load until the first visible crack formed. TI 
combined loads resulting in the first visible crack are plotted on the combi 
load curve in Fig. 16. The controlling crack, in the case of pipe #1838, occu: 
red in the coating and started at the upper quarter point (45° from the tomy 
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FIG, 15.—COMBINED LOADS ON TEST PIPE NO, 1837 
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FIG, 16,—COMBINED LOADS ON TEST PIPE NO, 1838 
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the pipe) and descended diagonally downward toward the springline (90° fron 
the top of the pipe) of the pipe. 

Next, as in pipe #1837, this test specimen was subjected to three combina. 
tions of internal pressure and three-edge bearing loads in which these load: 
were increased simultaneously to 2.0, 2.5, and 2.75 times the normal operatin; 
loads of 160 psi internal pressure and 6175 lb per linft three-edge bearing load 
The maximum crack width at each of these loads was recorded. This informa 
tion is also recorded on the combined load curves of Fig. 16. This specime! 
was then subjected to a series of increasing internal hydrostatic pressures wJ 
to 300 psi for the purpose of strain recordings, the results of which will bepre. 
sented later. 


FIG, 19.—FAILURE OF PIPE NO, 1840 UNDER HYDROSTATIC 
TEST PRESSURE OF 605 PSI 


The final test on pipe #1838 consisted of increasing the three-edge bearin 
load (no internal pressure) to ultimate failure and recording the cracking be 
havior throughout. This information is likewise recorded in Fig. 16. Previou 
to this last test to ultimate failure this pipe had been subjected to numerou 
high combination loads including one above the theoretical ultimate strengt 
of the pipe. This resultedin a substantial increase in the flexibility of the pip 
as is shown by the vertical diameter deflection in Fig. 18. In spite of this th 
three-edge bearing strength under each crack condition was greater than th 
conservative design strengths used for standard combined load calculations 

Pipes #1839 and #1840.—These twotest specimens were subjected to an in 
creasing internal hydrostatic pressure (no three-edge bearing load) upto burst 
ing pressure of 620 and 605 psi respectively. The cracking behavior of th 
specimen was observed and is shown on the combined load curves of Fig. 1 


ly 9 PRESTRESSED PIPE 93 


long with the results of the three-edge bearing load tests on pipe #1835. Fig. 
9 shows test specimen #1840 at the instant it failed under hydrostatic test. 


TEST PROCEDURE AND RESULTS BASED ON STRAIN RECORDING 


Pipe #1834.—Pipe #1834 was constructed with SR-4 strain gages distributed 
adially through the walls of the pipe at three different circumferential sections. 
fhe radial distribution of the SR-4 gages is shown in Fig. 10(a) and the spacing 
f the circumferential sections is shown in Fig. 10(b). The circumferential 
sections were spaced at 45° so that the strains recorded during the combined 
oads werethose occurring at; (1) the top of the pipe directly under the three- 
,dge bearing loads; (2) the quarter point, 45° from the top of the pipe; and (3) 
he springline or 90° from the top of the pipe. 

The strains occurring at these locations in pipe #1834 were recorded during 
he combined load applications which traced the 0.9 design curve and the de- 
sign curve as shown previously in Fig. 13. Since the strain recordings in the 
wo series of combined loads were of essentially the same pattern but slightly 
lisplaced in magnitude, only the records obtained during the tracing of the de- 
sign curve will be presented. Some of the gages were damaged during con- 
struction of the test specimenor for some other reason failed to record properly 
juring the test. Forthis reasonthere is not arecord for eachof the SR-4 gages 
shown in Fig. 10(a). Fig. 20 shows a plot of the strains recorded during the 
ypplication of the combined loads on pipe #1834 which traced the design curve. 
The letters on the abscissa of the strain curves denote the particular load ap- 
jlication, as shown on the combined load curve. 

The solid lines in Fig. 20 represent the strains as recorded by the SR-4 
rages. The dotted lines are computed values of the strain obtained froma theo- 
retical analysis of that particular section under a combined internal pressure 
und external three-edge bearing load. For the computed curves an elastic 
nodulus of 28.5 x 108 psi was used for the steel and 3.7 x 106 psi for concrete. 
As seen on the curves of Fig. 20 the computed and recorded strains in general 
‘ollow the same pattern. The displacement of the curves representing gage #7 
at the top of the pipe is difficult to explain unless it is a zero shift of the re- 
order that occurs only when the pipe is loaded. The larger variation between 
he computed and recorded strains for the concrete is to be expected because of 
he non-homogenity of the material. 

Pipe #1837.—Pipe #1837 had SR-4 strain gages mounted on the prestress wire 
und the concrete coating at the springline on both sides of the pipe. There were 
wo gages on the prestress wire and two on the concrete coating at each location. 
Type A-7 SR-4 gages were used for mounting on the coating and Type A-9 for 
mounting on the prestress wire. 

Strain recordings were taken on pipe #1837 during the combined loads to 
produce the first visible crack as was previously shown in Fig. 15. The strains 
resulting from this series of combined loads is shown in Fig. 21. The strains 
shown in the upper diagram of Fig. 21 are plotted directly above the combined 
oads, plotted in the lower diagram, which produced those strains. Since there 
was reasonable agreement between the four SR-4 gages on the prestress wire. 
he average strain recorded by these gages has been plotted. The average of 
the two coating gages on each side of the pipe have been plotted separately. 
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fe) 100 200 300 


95 


96 


(microinches / inch) 


STRAIN 


THREE-EDGE BEARING LOAD 


November, 1960 HY 


600 


-® 
Average of Four Gages on Prestress Wire - 
500 “AY ne - 
\- @ 
Si 
eo.|0 
400 dee 
ae Oa Rees 
en Peat a Meter o 
oOo _-@ 
hee 
@ 
300 
Oi jE tae Oo ee 
Cae ee Poe 


200 Average of Two Coating Gages on One Side 


Average of Two Coating Gages on the Other Side 


Final Zero Readings Initial Zero ee 


100 


F E D CG B A 
LOADING 
20 See 5 
en rs 
Ds 
15 meee 
eC ex 
— 
\ 
g ai, ‘0, 
= 6 First “Visible” Crack Curve L~ 
N\ ag 
ao a 
2 Letters Denote Loading NPA 
S 5 S 
3° Load "F" Produced "Incipient" Cracking Only \ 
= \ 
\ 
(e) 


) 100 200 300 
INTERNAL PRESSURE  (Ibs./ in2) 


FIG, 22,-COMBINED LOAD, STRAIN DATA FOR PIPE NO, 1838 


¥9 PRESTRESSED PIPE 97 


onsiderable shift in the zero readings of the coating gages is shown by the 
easured strain following removal of the last load (load K). 

This strain data was recorded when the first visible (0.002-in. x 12-in.) 
rack formed, which, in all cases occurred in the coating. Since any cracking 
" the concrete coating will necessarily relieve some of the tensile stress in 
ie coating, the strains as recorded by the coating gages would be affected by 
ie location and general characteristics of the coating cracks. 

Pipe #1838.—This test specimen was instrumented with SR-4 strain gages 
1) exactly the same manner as pipe #1837 which has been discussed in the pre- 
ious section. Again, as with pipe #1837, strains were recorded during the ap- 
lication of the combined loads to produce the first visible (0.002-in. x 12-in.) 
rackas was discussed previously and the results shown in Fig. 16. Thestrains 
scorded during this series of tests are shown in Fig. 22. Again the recorded 
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FIG. 23.—-HYDROSTATIC LOAD STRAIN DATA FOR PIPE NO. 1838 


‘rains fromthe four SR-4 gages on the prestress wire have been averaged and 
.e two recorded strains on the concrete coating on each side of the pipe have 
2enaveraged. The strain pattern exhibited by the prestress wire in pipe #1838 
;nearly the same as that exhibited by pipe #1837, However, the indicated 
rains of the concrete coating on pipe #1838 are considerably greater than on 
ipe #1837. As discussed in the previous section the strains measured by the 
yating gages could be affected by the cracking characteristics of the coating. 
his is apossible explanation for the large variation in strains indicated by the 


yating gages on the two pipe. 
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In both cases (pipe #1837 and #1838) the strain in the coating remains sub. 
stantially constant while the prestress wire strain decreases with a decreasé 
in internal pressure. 

Following the increasing combination loads on pipe #1838 (shown in Fig, 16) 
strain recordings were taken onthis specimenas the internal hydrostatic pres- 
sure was increased from 0 to 300 psi in 30 psi increments with no three-edge 
bearing load. The resulting strains are plotted in Fig. 23. Since there was 
good agreement between the various gages, the average, of the four prestress 
wire gages and the average of the four coating ggaes are plotted. Shown in Fig 
23 also are; (1) the theoretical strain - pressure relationship for the prestress 
wire with no coating effect with and without all the inelastic losses occurring 
after the core is prestressed; (2) the theoretical strain pressure relationshiy 
for steel only. The displacement of the coating strain curve to the left of the 
prestress wire strain curve can possibly be explained by the stress relief ir 
the concrete coating caused by the cracking of the coating. It must be remem- 
bered that previous to thistest pipe #1838 had been subjected to numerous load- 
ings which caused cracks in the concrete coating. 


CONCLUSIONS 


The prestressed concrete embedded cylinder pipe tested in this program 
were designed for operation at 160 psi internal pressure and 6175 lb per lin fi 
three-edge bearing load. This represents an equivalent trench condition of 1! 
ft of cover with ordinary bedding in a trench whose width is 3 ft greater thar 
the outside diameter of the pipe. These normal operating conditions fall on the 
0.9 design curve which is 10% below the design curve. 

The test results illustrate the conservative design of prestressed concrete 
embedded cylinder pipe, and clearly demonstrate the adequate stress crack 
control provided by the present design. No stress cracking was observed or 
the lining or coating for any combined load falling inside the design curve. It 
fact as pointed out under the heading “Test Procedure and Results Based or 
Cracking Behavior” there is a considerable margin of safety in the incipien 
crack condition. A combined load of 400 psi internal pressure and 15,400 Ik 
per lin ft three-edge bearing load was required to produce the first critica 
cracks (0.005-in.) in the concrete. This is a combined load equal to 2.5 times 
the design operating conditions for both internal pressure and external load, 
Upon removal of the excessive loads, these cracks closed to the point where« 
they were only visible through the 40X microscope. This illustrates that mino1 
cracking as a result of stresses caused by live loads or water hammer whicl 
may temporarily exceed normal design conditions is not serious, since the 
stress crack will close when the excess load is removed. 
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TRANSLATIONS OF FOREIGN LITERATURE ON HYDRAULICS 


Third Progress Report of the Task Force on List of Translations 
of the Committee on Hydromechanics of the Hydraulics Division 


INTRODUCTION 


In an attempt to bring tohydraulic engineers an up-to-date listing of trans- 
ations of foreign literature on hydraulics, itis the intent of the task force pre- 
aring this report to issue, from time to time, progress reports containing 
uitable items. Interested readers are urged to submit discussions to (a) add 
) this list, (b) offer suggestions for improvement, and (c) otherwise assist the 
ask force in fulfiling its aims. The present list is the third issued since the 
ublication, in 1957, of ASCE Manual 35, entitled “A List of Translations of 
‘oreign Literature on Hydraulics.” The first list appearedin the August 1959 
ournal of the Hydraulics Division and the second in the June 1960 journal. 
‘his list is to be considered as an addendum to Manual 35 and will be included 
1a revision of that Manual when the time is deemed appropriate. 


NONYMOUS. “Model Experiments in Hydraulic Engineering.” Statens Vatten- 
fallsverk, pp. 75-76, 1951. (On file at Army Map Service.) Ref. No. TC 89, 


F 97, Eng. Tr. 


NONYMOUS. “Model Experiments in Hydraulic Engineering.” Statens Vatten- 
fallsverk, pp. 73-74, 1952. (On file at Army Map Service.) Ref. No. TC 89, 


| eee tae fi 

$97, Eng. Tr. 
NONYMOUS. “Model Experiments in Hydraulic Engineering.” Statens Vatten- 
_fallsverk, pp. 80-81, 1953. (On file at Army Map Service.) Ref. No, TC 89, 


§ S97, Eng. Tr: 


en Tne ab itted | 
-Note.—Discussion open until April 1, 1961, Separate discussions should be submitted 
or the individual papers in this symposium, To extend the closing date one month, a 
mitten request must be filed with the Executive Secretary, ASCE. This paper is part of 
copyrighted Journal of the Hydraulics Division, Proceedings of the American Society 


f Civil Engineers, Vol. 86, No. HY 9, November, 1960. 
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ARMANET. L. “Genissiat (Power Station) Turbine Butterfly Valves” (Vannes- 
papillon des Turbines), From book entitled “Genissiat” published by La 
Houille Blanche, pp. 199-219. Tr. by P. J. Colleville, (On file at the Hy- 
draulics Division, Fritz Engineering Laboratory, Department of Civil En- 
gineering of Lehigh University.) Ref. Tr. No. 1, 1953. 


BELOVA, M. N. “The Theory of Pneumatic Wave Damping” (K voprosu o 
teorii pnevmaticheskogo volnogasheniya). Gidrotekh. Stroit, Vol. 1956, No. 
5, pp. 46-48. Tr. by Technical Information and Library Services, Minstry 
of Supply, London, England. (ASTIA Order No. AD 215819). 


CHECHEL, N. S. “The Effect of the Number of Runner Blades on the Power 
Output of High Speed Variable Pitch Blade Turbines.” Transactions of the 
Leningrad Polytechnical Institute of M. I. Kalinin, Ministry of Higher Educa- 
tion of the USSR, No. 193, pp. 38-50, 1958. (On file at DSIR, LLU.) Ref. 
RTS 1327. 

CONSTANTINESCO, GOGU, AND HADDON, WALTER. “Means for Checking 
the Effect of Inertia in Fluid Columns.” Tr. of Dutch patent 4,270, 16 pp, 
22 July 1959. (On file at LC and SLA) Ref. 59-18792. 


DECHEV, V. I. “An Investigation of Profile Cascades of a Reversible-flow 
Pump Turbine in a Cavitation Tunnel.” Energetika, Vol. 1958, No. 10, pp. 
130-143. (On file at DSIR, LLU) Ref. RTS 1346. 


DREYER, A. A., AND SELITSKAYA, Ye. S. “USSR Oceanographic Research; 
Measurement of Wind-raised Waves During the Second Voyage of the Com- 
plex Antarctic Expedition (and) Study of Variations in Sea Level.” Mezh- 
dunarodnyy Geofizicheskiy God. Informatsionnyy Byulleten (USSR), Vol, 
1959, No. 7, pp. 57-64, (On file at OTS) Ref. JPRS L-1166-N, 60-11164. 


DUMOV, V. I. “Design of Centrifugal Pump Stages with Preconnected Axial 
Wheels Having Improved Anticavitational Characteristics.” Teploenerge- 
tika, Vol. 6, No. 6, pp. 35-39, June 1959. (On file at U. S. Army Engineer 
Research and Development Laboratories, Fort Belvoir, Va.) Ref. T-703. 


FILATOV, I. N, “Effect of Hub and Throat Ring Form on the Cavitation Pro- 
perties of Variable Pitch Blade Hydraulic Turbines.” Transactions of the 
Leningrad Polytechnical Institute of M. I. Kalinin, Ministryof Higher Edu- 
cation of the USSR, No, 193, pp. 66-74, 1958. (On file at DSIR, LLU) Ref. 
RTS 1328. 


FRANKE, P, “Theoretical Considerations Concerning Jet Contraction in the 
Case of Flow Under Gates” (Theoretische Betrachtungen zur Strahlkon- 
traktion beim Ausfluss unter Schitzen). Die Bautechnik, Vol. 33, No. 3, pp. 
73-77, March 1956. Tr. by J.C. Van Tienhoven, WES. (On file at WES) Ref. 
Tr. No, 60-2. 


GIERS, I. V., AND SRETENSKI, L. M. “Influence of Varying the Principal 
Dimensions of a Ship on its Wave Resistance” (Vliianie izmeneniia glavnykh 
razmerov korablia na ego volnovoe soprotivlenie). Prikladnaia Matematika 
i Mekhanika, Vol. 10, pp. 21-32, 1946. Tr. by M. H. Brode, D, Fisher, and 
A. C. Thorpe, DTMB. (On file at DTMB) Ref. Tr. No. 282. 


HASKIND, M. D. “The Oscillations of a Floating Form on the Surface of a 
Heavy Fluid” (Kolebaniia plavaiushchego kontura na poverkhnosti tiazheloi 
zhidkosti). Prikladnaia Matematika i MeKhanika, Vol. 17, pp. 167-178, 1953. 
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Tr. by M. H. Brode, D. Fisher, andA. E, Thrope, DTMB. (n file at DTMB) 
Ref. Tr. No. 283. 


LLARIONOV, V. I., AND BENDITSKY, N. D. “A Measurer of the Velocity and 
Direction of Flow (and) Measurement of the Temperature of the Surface 
Layer of Water in the OpenSea” (Struyemer - Izmeeritel' Skorosti iNapra- 
vleniya Techeniya (and) Ob Izmerenii Temperatury Poverkhnostnogo Sloya 
Vody v Otkrytom More). Meteorologiya i Gidrologiya (USSR), No. 6, pp. 
49-54, 1958. On file at LC and SLA) Ref. M168, M169; 59-19065. 


NDRI, E. “Experimental Research on Models of Constrictions for Surge 
Tanks” (Ricerche sperimentali su modelli di strozzature per pozzi piezome- 
trici). L' Energia Elettrica, Vol. XXXIV, No. 6, pp. 554-569, June 1957. Tr. by 
J. C. Van Tienhoven, U.S. Army Engineer WES. (On file at WES) Ref. Tr. 
No. 60-3. 


VANOV, A. A. “Some Conclusions from an Analysis of Observational Data on 
Windexcited Disturbances” (Nekotoryye Vyvody iz Analiza Dannykh Nablyu- 
denii nad Vetrovym Volneyem). Akademiya Nauk SSSR. Morskoy Gidrofizi- 
cheskiy Institut. Trudy, Vol. 8, pp. 44-62, 1956. Tr. by Robinson. (On file at 
IC and SLA) Ref. 59-22613. 


AFFRY, P. AND HOURS, R. “The Study of Littoral Drift Using Radioactive 
Tracers.” Cahiers Oceanographiques, pp. 475-498, 1959. (On file at Army 
Map Service) Ref. No. 339696. 

ALIKHMAN, L. E. “Heat Exchange in the Vicinity of the Forward Stagnation 
Point in an Axially Symmetrical and Two-dimensional Gas Flow” (Teploob- 
men v okrestnosti perednei kriticheskoi tochki v osesimmetrichnom i plos- 
kom potoke gaza). Akademiia Nauk SSSR, Izvestiia, Otdelenie Tekhniches- 
kikh Nauk, Vol. 8, pp. 11-16, 1955. (On file at the Ministry of Aviation, Gt. 


Brit.) Ref. N-45780X. 


OLMOGOROV, A. N. “The Breaking Up of Droplets in a Turbulent Stream.” 
Akademiya Nauk SSSR., Doklady, Vol. 66, No. 5, pp. 835-838, 1949. On file 
at ATS) Ref. ATS-07L37R. 


ONONKOVA, G. Ye. “Measurement of Waves on the Surface of the Sea by 
Underwater Pressure Transmitters.” Nauchaye Doklady Vysshey Shkoly. Fi 
iko-Matematicheskiye Nauki, No. 2, pp. 228-234, 1958. Onfile at Ministry 
of Defense, Joint Intelligence Bureau, Division of Scientific Intelligence, 
London, England) Ref. Tr. No. 520. 

OSTERIN, S. I., SHEYNIN, B. 1, AND KATABRZHIS, A. K. “Experimental 

Characteristics of Stratified Flow in a Steam-water Mixture in a Straight, 
Horizontal Pipe” (Opytnyye kharakteristiki razdelennogo techeniya parovo- 
dyanoy smesi v pryamoy gorizontalnoy trube). Teloenergetika (USSR), Vol. 

3, No. 1, pp. 22-26, 1956. (On file at LC or SLA) Ref. 59-16543, CTS No. 

345. 

OVALEV, N. N., et al. “Investigation of the Pressure and Stresses Acting on 

the Runner Blade of Hydraulic Turbines at the Narva Hydroelectric Power 

Station,” Energomash, No. 1, pp. 29-32, 1959. (On file at DSIR, LLU) Ref. 

RTS 1113. 
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KRASSOV, I. M., AND TURBIN, B.G. “A Possibility of the Determination of 
the Axial Hydrodynamic Force in a Valve” (Obodnoy vozmoshnosti oprede- 
leniya osevoy gidrodinamicheskoy sily na zolotnike). Avtomatika i Tele- 
mekhanika (USSR), Vol. 19, No. 3, pp. 217-220, 1958. (On file at OTS) Ref. 
59-11958, No. H-2056A. 


KRISAM, F. “Influence of Volutes on Characteristic Curves of Centrifugal 
Pumps” (Der Einfluss des Leitvorrichtung auf die Kennlinien von Kreisel- 
pumpen). VDI Zeitschrift, Vol. 94, No. 11-12, pp. 319-322, April 1952. Tr. 
by A. Ostapenko and J. B. Herbich. (On file at the Hydraulics Division, 
Fritz Engineering Laboratory, Department of Civil Engineering of Lehigh 


University.) Ref. Tr. No. T-5. 


KRYANIN, I. R. “Cavitation Erosion in the Blades of High-power Hydraulic 
Turbines.” Energomashin, No. 3, pp. 14-18, 1955. (On file at DSIR, LLU) 
Ref. M. 943. 


LAMBOSSY, P. “Forced Oscillations of an Incompressible and Viscous Liquid 
in a Rigid and Horizontal Tube. Calculation of the Force of Friction” (Oscil- 
lations forcées d'un liquide incompressible et visqueux dans un tube rigide 
et horizontal. Calcul de la force de frottement). Helvetica Physica Acta 
(Switzerland), Vol. 25, No. 4, pp. 371-386, 1952. (On file at SLA) Ref. 59- 
17964, Tr. No. 192. 


LAUCK, A. “The Nappe Over a Weir” (Der Ueberfall uber ein Wehr). Zeit- 
schrift fur angevandte Mathematik und Mechanik, Vol. 5, No. 1, pp. 1-16, 


Feb. 1925. (On file at WES) Ref. Tr. 59-6. Tr. by J. C. Van Tienhoven. 


LAVRENTIEFF, V. M. “Theory of the Ideal Cavitating Propeller” (Teoriia 
ideal'nogo kavitiruiushchego dvizhitelia). Doklady Akademii Nauk SSSR, 
Vol. 50, pp. 89-93, 1945. Tr. by B, V. Nakonechny, DTMB. (On file at 
DTMB.) Ref. Tr. No. 279. 


LETERRIER, G. “Behavior of Rock in Steel-lined Tunnels; Results of Meas- 
urements Made at Randens, Montpezat, Brévierés, and Serre-Pongon” 
(Le comportement du rocher dans les galeries blindées; résultats des 
measures effectuées a Randens, Montpezat, Bréviéres, Serre-Poncon), 
La Houille Blanche, Vol. 11, Special No. A, pp. 144-170, March-April 
1956. Tr. by J. C. Van Tienhoven, U. S. Army Engineer WES. (On file at 
WES) Ref. Tr. 60-4. 


MAMEDOV, G. A. “Crossflow of Fluids Between Two Adjacent Beds of Differ- 
ent Permeability.” Izwestiya Vysshikh Uchebnykh Zavedeniy. Nefti Gaz 
(USSR), Vol. 2, No.4, pp. 49-55, 1959. (On file at ATS) Ref. ATS-08L35R, 


MEIXNER, J.“Flows of Fluid Media with Internal Transformations and Com- 
pression Viscosity” (Str6mungen von fluiden Medien mit inneren Unwand- 
lungen und Druckviskositat). Zeitschrift fur Physik (West Germany), Vol 
131, No. 4, pp. 456-469, 1952. (On file at SLA) Ref. 59-17940, Tr. No. 188 


MIKHAYLOV, I. G., AND SHUTILOV, V. A. “A Simple Method of Detectiny 
Cavitation in Liquids.” Akusticheski Zhurnal, Vol. 5, No. 3, pp. 376-378 
1959. (On file at DSIR, LLU) Ref. M. 1504. ; 


OBUKHOV, A. M. “Energy Distribution in the Spectrum of Turbulent Flow.’ 
Izv. AN SSSR Ser. geogr. i geofiz, Vol. 5, No. 4-5, pp. 453-463, 1941. (Or 
file at DSIR, LLU) Ref. M. 1375. 
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ISTOSLAVSKII, I. V., AND GRUMONDZ, T. A. “On the Relation Between the 
Generation of a Lift Force on a Wing and the Character of the Flow in the 
Boundary Layer.” Izvestiia Vysshikh Uchebnykh Zavedenii, Seriya Avia- 
tsionnaiia Tekhnika, No. 1, pp. 27-36, 1958. On file at NASA) Ref. NASA 
TT F-26. 


ADIONOV, S. I. “Wave Dissipation by Compressed Air (Pneumatic Break- 
water).” 47 pp, “Ocean Transport”, Moscow, 1958. Tr. by Hildegard 
Arnesen, University of California, Dept. of Engineering Technical Report 
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at U. S. Army Research and Development Laboratories, Dept. of the Army 
Corps of Engineers, Fort Belvoir, Va.) Ref. Tr. No. T-646. 
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RESISTANCE EXPERIMENTS IN A TRIANGULAR CHANNEL# 


Closure by Ralph W. Powell and Chesley J. Posey 


RALPH W. POWELL,?! F, ASCEand CHESLEY J. POSEY,22 F. ASCE.—The 
thors wish to express their appreciation for the contributions of the discus- 
rs. They also wish totakethis opportunity to report on additional tests which 
ve been made since the original paper was prepared. Table 3 gives the re- 
lts of 19 more runs on the smooth channel and 27 runs on the channel with 
2 walls coated with paraffin. This was'a mineral wax withsmall polyethylene 
ntent, suchas is used for waterproofing paper milk cartons, which was heated 
d rolled onto the inside surface of the flume with painting rollers. As the 
uterial solidified upon contact with the relatively cool surface of the flume, 
became slightly tacky so that the rollers pulied it up into a pebbled surface 
ich like that of fairly smooth concrete. Plaster casts of this surface were 
<en, and a method found whereby the casts could be hardened without distor- 
n. Equipment, designed to measure this roughness and analyze the meas- 
ements, was built, but so far it has not been perfected. 

The numeration of runs in Table 3 is not continuous with that of Table 1. 
1 runs 158-174, a false bottom was installed in the channel, making it a tra- 
zoidal channel. The results for the trapezoidal shape will be reported else- 
ere. 

Only one of the 19 runs on the smooth triangular channel had a Froude num- 
rover 0.85 (No. 202 with 0.88), and it showed nogreater resistance than the 
1ers, so that these runs are all classed as tranquil flow. The average n by 
inmning’s formula was 0.0092, the same as for the 23 runs already reported. 
e@ average discrepancy was slightly more, however, bringing the average for 
» entire 42 runs with tranquil flow up to 3.8%. 

The Colebrook formula for circular pipes is 


45, jira Sra Sed 
Ubtituting 2 R for ro this reduces to 


1 € 2.52 
—=-2 08 (seh RF | 
aie 14.83R Rv 


. Ackers says that for the triangular channela better fit is obtained if the 2 
changed to 2.2. (He also uses 14.8 instead of 14.83 and 2.51 instead of 2.52; 
y 

a May, 1959, by Ralph W. Powell and Chesley J. Posey. 


21 Prof, Emeritus of Engr. Mechanics, Ohio State Univ., Columbus, Ohio. 
* Dir., Rocky Mountain Hydr, Lab, and Head, Dept. of Civ. Engr., State Univ. of Iowa, 
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TABLE 3 - SUMMARY OF TEsT RESULTS (CONTINUATION OF TABLE !) 
RUN DATE 
NO. 


+ 2-58 
92-58 
9-2-58 
9-2-58 
9-2-8 
9-2-4 
9-2-8 
9-2-8 
9-2-8 
9-2-3 
9-8-58 
9-8-58 
9-8-58 
9-8-58 
9-5-8 
9-8-58 
9-5-58 
7-5-3 
9-35-58 
73059 
730-59 
73059 
72759 
727-59 
7AT-59 
7-27-59 
720-59 
7-18-59 


& 720-59 


203 
204 
205 
206 
207 


8-4-9 
8-3-59 
8-3-59 
8-359 
8-1-59 


& 8-3-5) 


208 
209 
210 


2il 


2l2 
213 
214 
215 
2l6 
2I7 
218 
219 


81-59 
810-59 
g:10-59 
8-4-5) 
65-59 
6-559 
88-59 
12-59 
81959 
8-8-59 
71-59 
871-59 
a71-59 
9- 1-59 
82659 
1-99 
827-59 


bc 828-59 


220 


9-1-59 


ae Wee Y, Q Frovoe R , CHEZy Mans 
TIME ROUGHNESS cS = Hes SLOPE AAT C £ 7 
11:00-11:73 PARAFFIN Y 04 a s 0 -0153 0-34 0-42] 1-62 1336 88-7 -0327 -Oll 
10:31-10:46 « V 0686 18 0153 0-439 0.69 1-89 233 103.3 0241 -OIL 
956-1014 V 0.889 11-4 -0153 0-560 1.8652 1-97 45 107-5 -022 -O10 
9:19- 9:35 V 1-084 10.6 0153 0.6% 3.04 189 446 102.8 -0243 -Oll 
8:46- 9:06 V wachos 0153 0.762 4.22 2-07 553 112-9 -O202 -Oll 
4:AT- 4:59 1 0-311 14-5 0152 0.776 4.55 2:14 63 117-2 -0187 -OlE 
4:23-4:37 14 0.459 14-7 0152 1-006 8-18 2.01 926 109-8 -02I3 -Oll 
3:57-4:13 «~~ 0.600 14.9 0152 1-165 12-27 2-09 1207 114.4 -019%6 -Oll 
3:30-3:44 «= 140-740 15-0 0152 1.2% 1689 2.20 1497 120-5 -0I77 -Olt 
2:59-3:16 " 50-898 14-9 0152 1.434 22-72 2-30 1812 125-9 -Ole2 -Olt 
3:40-3:58 “  V 0-349 14.6 00874 0-313 0.1818 083 57 59-6 -0724 OI 
2:57-3:16 © Vi 0-500 14-4 -00874 0.405 0.442 1-06 1733 763 -042 OF 
1:25-2:17" V-O-7II 14-0 -00874 0-534 1-059 1-27 222 91-6 (306 Oli 
1:02-1:24 «©  V 0-973 12-3 -00874 0.683 232 1-50 36 108-5 0218 Ol 
3:35-3:55 4 0-301 13-5 00875 0-886 4.33 1-46 40 1055 021 -Ol 
10:08-10:30.  “  -V _1-275 11-3 -00875 0-891 4.56 1-52 532 109-3 -0215 -Ol 
2:49-3:13." ~~ L 0-439 13-7 -00876 1.088 7-65 1-55 781 III-2 -0208 -O1 
[1:55-2:32  » — . 0-623 13-7 -00876 1-335 13-00 1-58 1078 113-5 -0200 -O1 
10:5] -11:35 . 4 0-836 13-0 00876 1-565 20.36 1-66 1415 119-4 -0180. -0) 
{1:20-1:26 SMOOTH V 0-352 135 -00144 0.343 0-1856 0.67 59-8 119-4 -0180 -00 
10:41-11:00 = V 0-549 12-8 -0044 0-516 0-556 0-73 J/7-0 12-9 -0155 -00 
9:39-10:05 «= V 0-708 11-8 -00144 0.649 1-049 0:77 170-1 137-0 -0137 -00 
4:22-4:43  « VV 0.884 13-0 “00144 0-786 1-826 0-83 252 147-5 -O118 -O0 
2:52-3:54 «  V 1-228 13-2 -00144 1-106 4.15 0-81 410 143-0 -01% -O 
(0:20-3:28 =" ~——sL10-349 12-1 00147 1-247 5-41 0.78 462 «136-5 0138-00 
10:01-10:55 § 1s 0.596 13-0 -00I47 1-€8412-15 0-82 784 1446 -0123 -0 
3:04-4:24 «0-814 12-2 -00150 2-0) 19:55 0-85 1036 146-5 -O1I7 -00 
Hs 16-11:51 "010-897 11-9 -00150 2-11 22-69 0-68 1139 152-2 -DINT -00 
10:51- 2:29 
9:09- 9:44 ~~  -V_0-268 11-0 -00068 0-330 0-113] 0-45 35-4 116-6 -0189 -00 
3:09- 3:51 » VY 0-447 12-7 -00068 0-509 0-335 0-45 71-2 116-9 -0188 -00 
2:18-2:47 «© VV 0-724 12:3 -00070 0-786 1-108 0-50 150-5 125 -0156 -00 
1:25-1:34 SV E-OLL 12-3 -00071 1-070 2-55 0-4 255 1362 -0139 -OL 
1:21-1:46 «SVs 139 T-8 00071 1-208 3-44 «0-54 «= 301: 1356-0140 -Of 
9:15-10:50 
9:55-11:05 «© V1-301 H-3 09071 1-362 4:80 0-55 367 140-1 -0131 -Of 
11:39-(2:1 «4 0-480 13-5 -00074 1-704 8-75 0-58 S67 142-7 -0126 -Ol 
9:09-1:18 »  .90-564 11-8 -00074 1-878 11-17 0-58 6B 1428 -0126 -Of 
4:19- 4:57 
8:53-1:55 
3:22- 3:36 » Ls 0-669 12-0 -00074 2-09 14-49 0.57 737 141-9 -012B -O 
9:30-10:00 
2:49- 3:02 
9:49-1:12 
10:38- 12:04 «U3 0-817 12-5 -00074 2-29 19-66 0:62 923 153-1 -O10 -O 
W:57- 3215 PARAFFIN V 0-590 13-6 -00074 0-699 0-665 0-41 104-9 100-6 -024 -0 
10:40- 11:30 V 0-80! 12-5 -00074 0-904 1-426 0-46 169-6 113-4 -0200 -0 
9:14- 9:49 nV 1-047 10-7 -00074 1-160 2:79 0-48 245 = 118-9 -0182 -0 
4:59- 5:37 » 4-313 13-2 -0007%4 1-396 4-60 0:50 361 123-4 -0169 -0 
2:01-4:53 SV s«A-268 14-2 -00074 1-421 4:50 0-47 355 115-8 -0192 -0 
4:03-4:45 «~~ L 0.469 13-7 00077 1-752 8-45 0-52 534 126-1 -0162 -0 
xe an “ U5 0-664 13-9 00077 2:15 14-32 0:53 743 128-1 -0157 -0 
3:26-3:39 «0-665 13-9 00077 2-16 14.36 0-52 742 127-3 -0159 -0 


| 


pe 


4 
, 


' 


_ clam., Denver, Colo., 1953, 
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but this is not a significant difference.) Therefore we have called 


tis o € 2,52 
ES asin See “i 5 \ 


Ackers’ formula, and have recomputed much of our data by this formula. It 
should be noted that when 2 is changed to 2.2, the value of € is automatically 
changed. 

The 19 new runs in the smooth channel were computed by this formula with 
€ = 0.00030 ft, which was the value Mr. Ackers found for the 23 runs we had 
already reported. The average discrepancy of the computed f from the ob- 
served f was 5.4%, which corresponds to an error of 2.7% inn or V. Combined 
with the 23 runs with an average error of 3.5% as found by Mr. Ackers, this 
makes a weighted average discrepancy of 3.1%. This is slightly less than the 
average discrepancy found by Manning’s formula or the formula suggested by 
the authors under the heading “Discussion of Results,” but the difference is 
hardly significant. 

For the tranquil flow in the triangular channel with the rolled-on paraffin 
surface, the average value of n was 0.0111 withan average discrepancy of 0.0003 
or 2.7%. Using Ackers’ formula an average value of « of 0.00269 ft is found. 
The f’s computed with this value agree with the observed f’s with an average 
discrepancy of 5.5% which is equivalent to a 2.7% discrepancy innor V. This 
is just the same as found for Manning’s formula, so it has no advantage. These 
data are for the 8 runs (Nos. 213-220) for which the Froude number did not 
exceed 0.53. No runs were made for this roughness with the Froude number 
between 0.53 and 0.83, but from the results with rectangular roughness it is 
assumed that tranquil flow would exist up to a Froude number of at least 0.70. 

When fitted to Nikuradse’s formula the € for these 8 runs is 0.0014 ft and 
the average discrepancy of the computed from the observed value of f is 6.3% 
which corresponds to a 3.1% error in V. But it is not surprising that this for- 
mula does not give good results here, because the flow is in the transition 


rather than the rough-channel range. The expression eR VE 

8V2R 
higher than 15.1 instead of 100 as it needs to have for “fully developed turbu- 
lence.” 

The 16 runs (Nos. 175-184 and 188-193) for which the Froude number was 
more than 1.27 gave an average n of 0.0111, just the same as for the tranquil 
runs. The average discrepancy was 0.0004 or 3.7%. Ackers’ formula gave an 
average ¢ of 0.00341 ft and the f?s computed from this value agree withthe ob- 
served f’s with an average discrepancy of 6.9% which is equivalent to a 3.4% 
discrepancy in n or V, not significantly better than Manning’s. The ratio of the 
observed f tothe f as computed by Ackers’ formula with ¢€ = 0.00269 ft was also 
investigated. The average of these ratios was 1.049 with an average discrep- 
ancy of 7.3%which is equivalent to 3.6% in n or V, not significantly better than 
Manning’s formula. It is surprising however that the f’s averaged 4.9% larger 
for rapid than for tranquil flow, while Manning’s n remained the same. 

Three of the runs (Nos. 185-187) had Froude numbers from 0.83 to 1.27. 
For these runs n varied from 0.0123 to 0.0173, « by Ackers’ formula varied 
from 0.00685 to 0.0332 ft, and the ratio of f to the f computed by Ackers’ for- 
mula with € = 0.00269 ft varied from 1.32 to 2.61. If we canbelieve the results 


TS Ee ocala Sel i tress eas aE een Tar RPeINL eT, 7a ea yal 
23 “Hydraulic Laboratory Practice,” Engineering Monograph No. 18, U.S. Bur. of Re- 


had no value 
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of these three runs, the resistance to critical flow is much greater than for 
tranquil or rapid flow. Runs 21 and 22 also fall in this group. Considering the 
waves which form, this high resistance is perhaps not surprising. It was sur- 
prising however that the run at a Froude number of 0.83 gave more resistance 
than the one with a Froude number of 1.06. (These are the Froude numbers 
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FIG, 17.—n AS A FUNCTION OF FROUDE NUMBER — 
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without the « term, but including it would not increase these values by more 

than about 3%.) | 
However, Mr. Ackers’ Fig. 1 indicates that even dividing the runs into three 

catagories, tranquil, critical, and rapid, is not sufficient. Fig. 17 is a plotting 
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of Manning’s n against the Froude number 
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While there is a great deal of scatter and the lines drawn are only tentative, 
it is evident that n is low at low Froude numbers, rises to a maximum as crit- 
ical flow is approached, and then falls off gradually rather than suddenly as the 
Froude number becomes more than one, Similar curves would result from 
plotting « against F as in Mr. Ackers’ Fig. 1. We thank him for calling this to 
our attention. 

In regard to his last question, the authors planto study other types of rough- 
ness, but probably will make no further tests on rectangular batten roughness. 


TABLE 4,—VALUES FOR a, b, AND c, 


Author 
(1) 


Zegjda Rectangular - 
Keulegan Wide 0.812 
Keulegan Trapezoidal 0.993 
Powell Wide 2.63 
Bretting = 
Thijsse Wide 0.981 
Owen Triangular 0.1515 
Sayre Wide 0.366 
Waterways Ex- 

periment Sta- Actual 

tion concrete 3,056 


Triangular 
Triangular 


Powell & Posey 


Ackers 


To answer Mr. Ackers’ other question, let the Colebrook formula be written 
in the more general form 
1 


i 2 


ea, € vi 52 1H R 
V0) 88 14.83R = RvE 


) and when R becomes very large Ene a log(2) 


(3 


+b and when e= 0, 3 =alog R Vf -c. In the literature we find the values 
f 
given for a, b, andc in Table 4. It seems, therefore, that while the value of a 
“may depend upon the shape of the cross section, it depends more on the judg- 
ment of the experimenter in fitting values of a, b, and c to his data, and, of 
‘course, also on the accuracy of his data. Also, it must be remembered that 
‘Colebrook’s formula is at least partially empirical and that it still remains an 
open question as to whether it is the best formula to represent the facts. ; 
_ From the purely empirical standpoint, it does not seem to us that Ackers 
formula fits our data as wellas some of the others. This is shown by the table 
of values of the average percentage discrepancy in the computed value of V 
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from the observed value for tranquil flow given in Table 5. For the smooth 
channel, Ackers’ gives the best fit, and for the paraffin lining it ties with Man- 
ning’s. But forthe 12 in., 6 in., and 3 in. spacing of battens it gives the poorest 
fit. We have not tested the fit for rapid flow except for the paraffin lining. 
Here the average percentage discrepancy inthe velocity by the three formulas 
was 3.7 by Manning, 3.6 by Nikuradse, and 3.4 by Ackers. 

We will turn now to Mr. van Malde’s discussion. He is very kind to refer 
to our work as accurate. We of course tried to be as accurate as we could, 
but the sources of error which we havé already described must have had a siz- 
able effect, as the scatter of points in Fig. 17 indicates. 

Varying the slope while keeping the depth constant is a good way of varying 
the Froude number. Since we have shown above that the resistance to flow de- 
pends on the Froude number, we agree with Mr. van Malde that w is influenced 
by S. Unfortunately, we had not foreseen this effect so that we did not make a 
wide variation in our slopes for some of the roughnesses. We thank Mr. van 
Malde for calling attention to Kirschmer’s papers. Wehad seenonly the earlier 
one in its French form,” and that not until we had torn out our 3-in. batten 
spacing. Since none of Kirschmer’s papers are widely available in the U.S. we 
are reproducing Fig. 14 from his French paper as our Fig. 18. The observa- 
tions are by Varwick, a student of Kirschmer’s. The ordinate | is our f/4and 


TABLE 5,—AVERAGE PERCENTAGE DISCREPANCY IN V 


Battens spaced, in inches 


Formula Smooth Paraffin 


Manning 


Prandtl- 
Nikuradse 


Ackers 


the abscissa is our R/4. With our apparatus it would have been very difficult 
to observe a large number of slopes with the same depth, but we are glad that 
Varwick was ableto do it. Similar results had already been obtained in Russia 
by Zegjda.24 

We agree with Mr. van Malde that these tests give no certain answer as to 
what formula for channel resistance is “best.” About as far as we can go is 
to say that for the range covered by these experiments, Manning’s formula 
with n constant for any given absolute roughness and any given Foude number, 
gives discrepancies no larger than might be expected due to experimental er- 
rors, Other formulas investigated do not fit the data well enough to warrant 
their use. This is especially striking since the authors had begun their study 
with the belief that a formula definitely better than Manning’s could be found. 
The wide variation of n for rough channels with changes in the Froude number 
is perhaps the most important finding of this study. a 


9 “Pertes de charge dans les conduits forcees et les canaux decouverts,” Revue Gen- 
erale de l’Hydraulique, Paris, France, Vol, 15, No. 51 (May-June, 1949), pp. 115-138, 
“Teoriia Podobiia i metodika rascheta gidroteknicheskikh modelei,” (Theory of 
Similarity and Methods of Designof Models for Hydraulic Engineering,) by A, P. Zegida, 
Leningrad, Union of Soviet Socialist Republics, 1938, 
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We believe one of Mr. van Malde’s remarks in the first paragraph of his 
discussion might be misinterpreted. For tranquil flow in the smooth channel 
we found 7’ = 0.93 which means that f averaged 7% less than by Prandtl’s for- 
mula. This means that our observed V’s averaged 3.5% greaterthan Prandtl’s 
formula would give, or that if our data are correct, the average error in V by 
Prandtl’s formula would be -3.5%. But for the rough channel, using the € which 
gave the best fit, the average discrepancy in V was 2.18% or say 2.2%, with no 
minus sign, as this is the average of the absolute values of the discrepancies 
without regard to sign. 

In a private communication to the authors, H. J. Koloseus, M. ASCE, has 
made a very thoughtful study of why the two equations given under the heading 
“Results from Rough Channel,” seem to be about equally satisfactory while the 
first makesf dependonly on the relative roughness ¢ /R (but with € a function 
of the batten spacing A), and the second makes it depend on A and fp which is 
a function of the Reynolds number. He introduces another dimensionless par- 


: V i..0 

ameter which he calls S, = v2¢8 € ™* (which by the way isthe square root of 
= } 

twice the expression X = gS (634) 2) used by Bretting.25 Koloseus shows that 


2 
= Ae LENS so that f could be expressed in terms of R andS, instead of 
R 

4R/e. Then if ina given set of data Sx is approximately constant,f seems to 
be a function of R alone, or of e/Ralone. He points out that with the 12 in. 
batten spacing there was very little variation in S,. The same was true with 
the 24 in. spacing. However with the 3 in. and 6 in. spacing there was quite a 
little variation in S, and he submits Fig. 19 which shows an interesting rela- 
tionship. He also points out that Nikuradse’s formula for rough pipe withe = 
0.00040 ft fits the smooth channel tranquil data with an average discrepancy in 
f of 6.6% which is somewhat better than the 3.5% for V (7.1% for f) which Mr. 
Ackers reported for his modified Colebrook formula. It may be added that Mr. 
Koloseus26 has made a very careful study of resistance to flow at high Froude 
numbers in a rectangular channel roughened by closely spaced cubes. 

The authors wondered why Buckley called his formula Manning’s formula. 
This matter has now been cleared up by one of the authors. 

The authors are sorry that there was no discussion of the “capillary rise” 
mentioned under the heading “Results from Rough Channel.” It seems that this 
little known phenomenon may account in whole or in part for the apparent in- 
crease in roughness in critical and rapid flow. 

During the summer of 1959 a very large number of pitot measurements of 
the velocity were made for uniform flow at a relatively large depth. In order 
to do this the flow had to be held steady for many hours at a time. This took 
less man-power than establishing new steady flow rates at frequent intervals. 
Advantage was taken of the long duration of the uniform flow condition and of 
the availability of personnel to determine water-surface and bottom profiles 
‘more thoroughly than before. Values were read every 2 ft along the length of 


f 


25 “A Set of Practical Formulas Based on Recent Experimental Research,” by A. E. 
Bretting, Report, 2nd Meeting of the Internatl, Assn, for Hydr. Research, 1948, pp, 399- 
415, 
26 pissertation, by H. J. Koloseus, Presented to the State University of Iowa, Iowa 
‘City, Iowa, in partial fulfilment of the requirements for the degree of Doctor of Philos- 
Py August, 1958, Presented at the October, 1958, ASCE Convention. 

_ 27 “History of Manning’s Formula,” by Ralph W. Powell, Journal of Geophysical Re- 


‘search, Vol, 65, April, 1960, pp. 1310-1311, 
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¢ flume with the exception of only a few points that could not be seen from 
é level stations. 

The results are shown in Fig. 20 compared with the values for the mean 
ater-surface and bottom profiles as ordinarily determined, that is, by the 
sual method with bottom elevations determined every 8 ft and water-surface 
levations every 24 ft. It can be seen that for this flow condition the ordinary 
ethod determined the slopes fairly well, but underestimated the depths by 
srhaps 0.06 ft or about 2.85%. The Froude number for this run was 0.88 which 
eant that surface waves hada relatively large effect. Together with the large 
sflection of the bottom due tothe heavy water load, the underestimation of the 
»pth in this test undoubtedly represents nearly an extreme value for the pres- 
it experiments. 

Attempts were also made to get transverse water-surface profiles, but set- 
ng gages to the fluctuating water-surface proved too difficult a task. With 
e slight irregularities of the flume obviously causing much larger water-sur- 
ce waves, they would be of doubtful value. 

The fact that the channel showed less resistance to flow than the correspond- 
g smooth pipe, has troubled both the authors and the discussers. However, 
ansford and Robertson28 reported that tests made at the University of Illinois, 
7 W. M. Owen, showed f for a smooth triangular flume less than for the cor- 
2sponding pipe. For values of R up toa little over 100,000, the tests gave 
= 0.295/R9.25, It is well knownthat inthis range of Reynolds numbers, Pran- 
Vs formula is practically equivalent to the Blasius formula f = 0.316/R9-25, 
herefore Owen’s ~ would be 0.295/0.316 = 0.934 which is almost exactly the 
936 given by the average ofour 42 runs of tranquil flow in the smooth flume. 

On the other hand Kirschmer? gives a graph showing tests on a smooth tra- 
2z0idal channel by Varwick and on a rectangular one by Raju. These indicate 
about 1.04 for the trapezoidal channeland 1.12 for the rectangular one. Kir- 
‘hmer feels that this is to be expected because the effect of secondary flow 
juld be greater in the rectangular section than in the trapezoidal. It would 
em that a triangular section should give » at least as great as a trapezoidal 
1e. 

The elevation of our channel bottom may have averaged as much as 0.01 ft 
‘low the values measured at the 8-ft panel points, so that the tabulated values 

Yn may have been as much as 0.01 ft too small. The 42 runs for the smooth 
ume with tranquil flow were recomputed witheach yp increased 0.010 ft from 
e value given in Tables 1 and 3. The resulting y was 1.007 instead of 0.936. 
was also found that the average discrepancy from the Prandtl-Nikuradse 
rmula for smooth pipes given by these assumed values was only 5.4% for f, 
nich is equivalent to 2.7% for V or n, compared with 3.3% using bottom ele- 
tions as recorded. But the close agreement with Owen’s results leaves some 
iestion as to whether this adjustment should be made. 


28 Discussion of “Open-Channel Flow at Small Reynolds Numbers,” by Wallace M, 
nsford and James M, Robertson, Transactions, ASCE, Vol, 123, 1958, p. 707. 
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THE SETTLING PROPERTIES OF SUSPENSIONS? 
Discussion by N. Claes H. Fischerstroem 


N. CLAES H. FISCHERSTROEM, ! M. ASCE.—Mr. McLaughlin’s paper is a 
uluable contribution to our knowledge of the sedimentation process, especially 
nce it includes results from determinations of the settling of suspensions, as 
ell as practical hints regarding the execution of the tests. 

Although, according to the introduction, the paper is limited to the settling 
-operties of suspensions and tothe use of the results in computing removals, 
author, in section 4, also discusses “practical engineering problems,” He 
akes some statements regarding the influence of “other factors” onthe design 
a settling tank. The writer has some questions to ask regarding the author’s 
Minions, and presents some further conclusions from the results. 

Mn section 2 the author starts with an analysis of the settling, based on the 
mount of suspended solids which passes a certain plane, calculated by the 
ean concentration and mean velocity of the solids under non-disturbing flow 
mditions. The equation arrived at covers the settling in both the free and the 
ndered zone. The main value of the equation, according to the author, is to 
10w what research and development are necessary for complete calculations. 
it are there really any new items of research mentioned in the paper which 
‘e not already known? However, the analysis is a good basis for the exten- 
ve discussion of the tube test procedure. 

The flow of water ordinarily takes place mainly inthe free settling zone and 
e equation is not very useful for the purpose of drawing conclusions regarding 
is part of the settling basin. Such an equation is still lacking. 

In the end of section 2 the author discusses scale models. Settling studies 

models cannot ignore the flow conditions. The writer believes that under 
ost actual conditions it is more important to study the flow conditions in the 
‘ale model and the settling in the tubes. If the flow is satisfactory, and only 
en (not necessarily ideal) will the sedimentation in the prototype approxi- 
ately confirm to that in the tube. 

Section 3 gives valuable information regarding results and discussion of 
sts by the pipette and the multiple-depth analysis. The latter will be discus- 
id here. 

In Fig. 5 the settling in a tube is illustrated by a concentration diagram. If 
2 dashed lines represent concentration at constant z/t, shouldn’t they all 
art at ¢ (z,t) = 100%? Figs. 7 and 9 represents the settling in a tube of ac- 

1 suspensions. Shouldn’t the full lines (time curves) start at ¢ (z,t) = 0% 
i dashed lines (z/t curves) start at ¢ (z,t) = 100%? As is seen from 
2 figures, the shape of the upper part of the curves will be different if the 
ait values are taken into consideration. This first part of the curves can 


2 December, 1959, by Ronald T, McLaughlin, 
+ Chf, Cons, and Research Engr., Vattenbyggnadsbyran (VBB), Stockholm, Sweden, 
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never be vertical but will always be nearly horizontal or more or less sloped 
which is important when drawing conclusions. 

From Figs. 7 and 9 we can see that the settling properties of the suspen- 
sions are bad. If there has been no flocculation period there will be a chemica 
reaction and we are not studying a settling process alone. It is known tha 
Fitch’s sedimentation test is carried out after dosing of the (whiting and) fer- 
risul by a rapid mix, but without flocculation of any kind. Is the author’s tes 
with betonite clay and alum also carried out without preceding flocculation ‘ 
Or if it is flocculated, what are the details (stirring, time etc) ? Is the ma. 
terial tested really a suspension? Isn’t it a suspension under formation by % 
process which is not only time-consuming but also influenced by turbulence 
concentration, pH a.s.o. The writer does not meanto say that such suspension: 
would not occur, but nevertheless these tests might be too specific to alloy 
conclusions regarding the settling laws. 

In section 4, part b, the author states “This approximation will often be bet. 
ter than making an elaborate study of the hydraulics of the tank while ignorin 
the properties of suspension.” Why ignore any of the important factors ? Thi 
underestimation of the hydraulics seems tobe not uncommon for the laborator; 
fellow. In contrast to the field engineer, he gets no experience from the distur 
bances from density currents, side currents, wind, turbulence etc. The settlin; 
properties are easily determined, but the hydraulic properties are not. Tha 
is perhaps why so many existing basins, and also so many basins under con 
struction, are unsatisfactory in hydraulic respects. 

In the same section, part d, the author states, “It will now be demonstrate 
that whether overflow rate or detention time should be used depends upon th 
nature of the suspension. Indeed, for a single tank the removal ratio may b 
closely related to overflow rate for one suspension, to detention time for 
second, and to neither for a third.” This sounds simple, but Mr. Fitch wa 
perhaps nearer the truth when, in a discussion regarding influence on settlin 
of Class II suspensions by the overflow rate, namely, the detention time, h 
wrote, “both are working.” Will it ever be possible that a sedimentation onl 
depends upon time ? And what does the author mean by “neither for a third? 

To illustrate the preceding, the author starts with a theoretical discussio 
of the concentration diagram. Overflow rate, he says, is determining if 2z, 
curves are straight and parallel tothe z (depth) axis and detention time is de 
termining if the concentration (time) curves are parallel to the same axis. | 
we deal with suspensions which settle, the latter will never be possible becaus 
all time curves but the one at t = 0 must start at the concentration @ (z,t) = ( 
In the thin surface layer, sedimentation results immediately in a complete re 
moval. In the surface layer there is no flocculation. Also in the next, shallo 
layers, clarification is more rapid than in the deeper layers. Take, for in 
stance, the author’s curve for T = 700 sec (Fig. 7). Ata depth of 10cm w 
have a removal of nearly 50%, at 20 cm of about 40%, and at a depth of 100 rs 
only 23% (at the same time or volume load). At the same detention time 
shallow basin will always give a higher removal than a deeper basin. It i 
consequently, even with a suspension with poor settling properties, more € 
ficacious to use, for example, a 0.2 in. basin with 5 times as large surface t 
touse a 1.0 mbasin. The shallow basin furthermore will have lower turbulenc 
and higher stability thanthe deeper basin; when we alternatively insert 4 tr 
in the deeper basin, we also exclude disturbance from wind action in 4/5 of t 
basin. 3 
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To illustrate the relation between overflow rate and time for certain kind of 
uspensions we can study the writer’s Fig. 1, based upon the results published 
y Fitch. This diagram, which is of the same type as the author’s though re- 
rranged for technical purposes, shows the residual after sedimentation: 

a. at constant settling time (volume load) and variable depth (overflow rate, 
ull lines) 


30 


N 
oO 


PER CENT RESIDUAL 


0 1 2 3 4 5 
DEPTH OF SETTLING CHANNEL, ft 


FIG, 1.—PERCENTAGE RESIDUAL VERSUS DEPTH OF SETTLING CHANNEL 


: -b. atconstant overflow rate and variable depth (settling time, dashed lines). 
: "The time curves show us that at constant time, 5 ehn; for example, the re- 
idual increases with the basin depth, or with a deeper basin with a smaller 
rea. From these curves we may draw the conclusion that the best way -- from 
sedimentation point of view -- to use a given volume is to give the basin as 
arge an area as possible. 

a 
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The overflow rate curves indicate that the residual decreases with increas: 
ing depth, which means with increasing time or a decreased volume load. Thi 
increased efficiency must be bought bya greater basin volume. The effect ob- 
served, which depends upon an accelerated settling because of flocculation 
shows us, however, that there is a certain relation between the overflow rat 
and settling time for every required settling efficiency. 

We can, for example, take a recovery process with 90% settling efficienc 
(10% residual) and a purification process with 2% allowed residual. Corre: 
sponding overflow rates and detention times will be. 


Residual Overflow rate per detention time 
10% 1.1/2.5 hr 0.4/2 hr. 0.2/1.5 hr 
2% 0.5/5 hr 0.2/3 hr 0.05/2 hr 


From the examples we can confirmthe often expressed relation: “the lowe: 
the overflow rate the shorter the time,” in order to reach the same removal 
The only possibility to increase the volume load and to reduce the basin siz 
is in fact by decreasing the overflow rate, environmental factors being un 
changed. Very low residuals, near 0%, which we now require in certain pro. 
cesses, seem only to be attained at very low overflow rates and very smal 
depths with these extreme suspensions (or by flotation, of course). This ca: 
be achieved by shallow single story basins or multiple story basins of the typ 
suggested by Camp. Which type of basin ought to be chosen depends upon th 
cost per basin unit volume and other factors, such as available space, stabilit 
requirements, wind conditions, and so on. 

We can also consider the overflow rate and time influence in this way. | 
the free particles of a suspension at the time t = 0, for example, before an 
flocculation can take place, have certain settling velocities, which we represer 
in a settling diagram, we can always determine the overflow rate for arequire 
removal at depth = 0 or time = 0 for this suspension. Because we have a cer 
tain depth or time, the actual removal will be increased, and to keep the re 
quired removal we can increase the volume load. Because the incluence ¢ 
flocculation will increase the settling properties only up to a limited depth 
above which the particle speed will disperse the floc again - this improvemer 
by increased time or depth will take place only to a certain limit. Increasin 
the time above this point will not increase the removal ratio, and any attemg 
to increase the overflow rate load will be a failure. From the diagram we ca 
see that the constant time curves will always lead to zero at decreased over 
flow rate; curves for constant overflow rate would finally be parallel to th 
depth axis. 

To summarize, it will always be possible to get a required removal, eve 
by reducing the time ~0, by lowering the overflow rate, but this will not k 
possible by too high an eMection rate, even if the time is increased to — rs 
Although both overflow rate and time can be said to be limiting factors, it fol 
lows that the overflow rate might be said tobe of primary importance when a 
increased volume load is intended. ] 

It seems most probable that the very shallow basin will have a future, pre 
vided it is designed on technically correct principles. Our knowledge regardin 
short-time settling is unsatisfactory. To determine the settling of a suspensi 
in a very short water column ina short period will require new iaboratdl 
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sthods. It does not seem to be possible to still the water and avoid turbu- 
ee rapidly enough in the “tube test.” The writer believes that we perhaps 
ist change to a continous flow method. 

The investigations made by the author, Fitch, and others during recent years, 
ve given us an idea of the relation between the overflow rate and time for 
me special types of suspensions. Further investigations regarding other 
eS are desirable. What seems to be especially valuable is what has pre- 
usly only been assumedor felt intuitively or possibly doubted or denied, that 
 flocculated suspensions of Class II A type, there is no direct proportionality 
capacity tothe settling surface area; the capacity increases, but at a slower 
te, the speed depending upon the type of suspension. The investigations also 
nfirm, however, that the conclusions regarding the extreme form of a basin, 
awn from the overflow rate theory, are correct; the best basin for all types 
Suspensions is still the one with depth and time 0 and surface o. The in- 
stigations seem, in this respect, to have shown quite clearly, as mentioned, 
it the only constructional way to increase the volume load is to decrease the 
oth, not to increase it. The conclusion that a settling process can be depend- 
- only upon the time, seems to be drawn from concentration curves lacking 
» important part just below the water surface. 

The prior discussion is not meant as a criticism, only to give a more com- 
‘te view on the conclusions drawn from the author’s settling theories and 
serimental results. We must be grateful to the author for his excellent work 
1 the writer agrees with his conclusions, with the exception possibly for 
. 1. The writer believes, namely, that when a choice must be made, it may 
more important to study the flow in the basin than the properties of the sus- 
asion. Finally such a choice depends, however, upon the knowledge and ex- 
rience of the designer. 
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EARLY HISTORY OF HYDROMETRY IN THE UNITED STATES 
Discussion by Ivan J. Moskvitinoff 


IVAN J. MOSKVITINOFF, 84 & ASCE.—This paper is a valuable contribution 
the history of engineering science in the United States in providing the col- 
ction and systematization of the information scattered in many publications 
id magazines. The author paid a tribute to the American pioneers in hydrom- 
ry and to the engineers and professors who built the scientific and technical 
ise of modern hydrometry. The paper contains important information on the 
irallel development of the hydrometric methods and devices in some European 
untries. 

Describing the investigation on velocity pulsation at Burlington, the author 
mcludes: “The influence of pulsation cannot be eliminated by short observa- 
ons at a great number of points as it is generally believed.” The results of 
e writer’s measurements of the velocity distributionin many irrigation chan- 
Is of different sizes and inthe Syr-Daria River, in Turkeston, Russia, in 
114-1915, led him to the following conclusion: the magnitudes of the rate of 
yw based on velocity measurements at many points in a cross section, with 
short period of observation, were sufficiently accurate in comparison with 
e rates of flow obtained by two-point method with a long period of observa- 
mm necessary to eliminate the influence of pulsation. 

The writer vigorously supports the author’s proposal to use a more logical 
rminology as indicated in Appendix I. 

The writer hopes also that the author will be able to fulfil his intention to 
iblish the history of the modern hydrometry after the First World War, for 
is would be a valuable contribution to engineering science. 


a January 1960, by Steponas Kolupaila. 
84 Giy. Engrg, Dept., The City College of New York, New York, N. Y. 
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ENERALIZED DISTRIBUTION NETWORK HEAD LOSS CHARACTERISTICS® 
Discussion by F. P. Linaweaver, Jr., John C. Geyer, and Jerome B. Wolff 


F. P. LINAWEAVER, JR., 21 A. M. ASCE, JOHN C. GEYER, 2“F. ASCE, and 
tROME B. WOLFF,22 M. ASCE.—There are many problems that continually 
‘ise in the operation and planning of a water distribution system that require 
hydraulic analysis. If the basic analyses of average daily, maximum daily, 
1d peak hourly conditions are availablefor each pressure zone of a water dis- 
ibution system, some of these problems can be answered by the use of propor- 
ynal loads as presented by Mr. McPherson. For systems without storage or 
ith storage at or near the pumping station, proportional load hydraulics ac- 
rately produces the desired answers. The author has used test data to show 
at fora single pumping station, and a reservoir elsewhere on the system, the 
ad loss between the pumps and any given point in the system is a function of 
e ratioof pumping rate to demand rate. The discoverythat such relationships 
<ist should prove very useful, particularly when analyses are made by long hand 
ethods or by use of digital computer. 

The method appears to have usefulness for solving problems other than those 
volving pumping versus storage relationships. One recurring problem is to 
termine the effect of proposed housing developments on the distribution sys- 
m. This would normally be a concentration of loading at a point and, there- 
re, would be a non-proportional load similar to a fire flow. The flows into 
- out of equalizing storage facilities are in effect non-proportional loads, and, 
erefore, a procedure might be determined whereby Eq. 2 can be used for ap- 
ying fire flows, etc. 

In making analyses involving a storage reservoir, the fact of discontinuity 
the situation when pumps are changed or when the reservoir either fills or 
npties, must be dealt with. This requires, in particular, that the time of fill- 
 andof emptying the storage reservoir must be determined, because at these 
mes there occur drastic changes in the hydraulic situation. Following of the 
elical variations in storage levels would appear to require the reverse ap- 
ication of Eq. 2, with appropriate accounting for pump operation and pumping 
ation discharge curves. 

In the caseof a more complex system, as shown in Fig. 7, the determination 
‘ the appropriate Qp/@a ratio to account for various pump curves of a par- 
cular pumping station, and the changes inlevels of tanks and standpipes would 
spear to be very difficult. The author states that “The computations for the 
1ecessive hourly balancing of pump-network-storage functions can now be 
srformed much more simply using Eq. 2,...” In Eq. 2, Qg, n and @ would be 
1own; in order to determine Dh, Qp must be determined. If there is more 


January, 1960, by M. B. McPherson. 

21 Asst. Civ, Engr., Bur. of Water Supply, Baltimore, Md. 

22 Prof. and Chmn. of Sanitary Engrg., Johns Hopkins Univ., Baltimore, Md. 
23 Asst, Dir, of Public Works, Baltimore County, Md. 
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than one pumping station or more than one tank or standpipe, the distributi 
of the quantities of water being supplied by each source would not be known. 
would appear that in order to determine the ratios of Q)/ Qg as shown in co 
umns (3) and (10) of Table 4, an analysis of successive balancing would be r 
quired similar to the one described in another publication, 5 but accounting fi 
the changes inlevels of tanks and standpipes. Thus, a thorough analysis wou 
still be required. It would be interesting to know the author’s procedure f 
determining columns (3) and (10) of Table 4 especially in the Number 4 seri 
of runs. 

In analyzing complex systems on the Baltimore Analyzers, a procedure fi 
balancing is used, whereby the operating levels of the storage tanks or rese 
voirs are determined on the initial test, and then the system is operated throu 
the day just as the pump operator or an automatic station would operate. I 
stead of using an hour-to-hour balancing, the period between tests is vari 
as needed to obtain results within the desired accuracy. 

The operating level is defined as the high water level during the early a.1 
after several days of the same pattern and magnitude of demand. This level 
obtained by setting the demand at the average for the day being analyzed, a 
by reading the difference inhead betweenthe sources and the storage facilitie 

The question, “Is it not feasible to consider writing a design program for 
digital computer which would permit direct determination of the best combin: 
tion of pipes to satisfy prescribed conditions defined by means of Eq. 2?”, is 
very Challenging one. It could be expanded to adetermination of the best con 
bination of pumps, pipes, and storage facilities by modifying Mr. McPherson 
Kq. 2 to include the criteria presented by George G. Schmid.24 
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SCOUR AT BRIDGE CROSSINGS® 


Discussion by D. V. Joglekar, W. J. Bauer, L. J. Tison, S. V. Chitale, 
A. Rylands Thomas, Mushtaq Ahmad, and Pier Luigi Romita 


D. V. JOGLEKAR. 16_tt is agreed that each river and each reach must be 
idied to understand its individual, almost personalized, characteristics and 
at scour at bridge piers is closely related to (a) the river conditions up- 
ream and downstream (such as indicated by flood hydrographs from year to 
ar); (b) the meandering tendency of the river that, in turn, depends on the 
tritus load carried by the river during floods, (c) whether the river is flow- 
> in its alluvial plain or has its sides and bed resistant to scour, (d) whether 
2 river is flashy or has sustained floods, and (e) the constriction of the riv- 
section caused by the bridge. The reach of the river also has an important 
aring on scour at bridge piers, for example, bridges on hill torrents, bridges 
alluvial plains with mild slopes, and bridges on the tidal reaches of the 
yer, In a braided river near a hilly region the river changes are so rapid 
at it is very difficult to estimate the waterway required for bridges across 
rious streams of the braided river. An illuminating case of violent move- 
ant of Manas River!? in Assam can be cited. In 1909, three bridges were 
ovided on the three tributaries of this river: 


1. on the Manas River which had 10 spans of 100 ft, 
2. on the Bholookadoba Branch - 2 spans of 75 ft, 
3. on the Beki Branch - 4 spans of 12 ft. 


1e to changes in river courses, the bridge had to be rebuilt several times on 
e Bholookadoba Branch between the period 1909 to 1945: 


ae. 2 of 75 ft, 
2. 3 of 75 ft, 
3. 1 of 250 ft and 2 of 150 ft. 


ianges in the Beki Branch were more violent and larger waterways had to be 
Ovided during this period as follows: 


a. 4x 12 ft, 

. 4x 20 ft, 

- 9x 20 ft, 

. 1x 50 ft, 1 x 30 ff 6 x 19 ft, 
1x 50 ft, 2x 40 ft, 6x 19 ft, 
. 2x 40 ft, 3 x 100 ft, 

mt x. L00ft. 


1e river is still uncontrollable in spite of various training measures. 


AQ Fwp 


@ February, 1960, by E. M. Laursen, 
16 Adviser, Central Board of Irrigation and Power, Poona, India, 
, «River Training,” Railway Board Technical Paper No, 334. 
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Experiments!8, 19 Were carried out at the Central Water and Power Researt 
Station, Poona, India, in 1938 and 1939 for finding scour round a single pie 
placed in the center of aparallel sided flume, embeddedin sand of mean dian 
eter 0.29 mm. The following relation was worked out: 


2 
d > 
we 3 0.78 
> = 1,70 (i dio lat ie; wile <ee s0)kbp B. ia) Plane weer (14 


in which b is the width of the pier, q, is the discharge per foot in the center « 
the flume upstream of the pier, and d, is the maximum depth of scour at tl 
pier below water level. As it was difficult to correlate this with the depth ¢ 
scour at prototype piers, due to q (intensity of discharge per foot) dependir 
upon the curvature of the river upstream - which varies from river to rive: 
it was considered desirable to study cases of actual scour in prototypes ar 
work out a general, empirical relationship. Besides, it has to be remembere 
that the angle of repose of bed material in the model and the prototype is tt 
same, hence, the extent of scour in plan in the vertically distorted model : 
found always relatively greater than in the prototype. This in effect reduce 
the discharge intensity at the pier due to greater dispersion of flow and henc 
the depths of scour obtained in the model would be relatively less. Data! 
were, therefore, collected for scour round Prats piers of various bridges cor 
structed in India and a general relationship!9, 0 was worked out as follows: 
1 


d, = (2) 0.473(Q/t) > = 2 Ditiaeey) 6) os, 8 voles ae 


in which Qis the maximum flood discharge in cu sec; d, is the maximum dept 
of scour below highest flood level; f is = 1.76 Vm, and m is the mean grade « 
bed materialin millimeters. In Eq. 15, a representative f value has to be use 
From bore data, values of f for each strata is to be worked out to ascertai 
that the anticipated depth is not based on the f value which is higher than thi 
appropriate at that depth. Recent advances infoundation engineering have mat 
it possible to take the pier foundation sufficiently below the maximum probab! 
depthof scour to provide adequate grip length. Where this cannot be done, hig 
level stone protection, though costly in the long run, has to be employed. Ar 
other advantage of deep foundation is that, because of increased side frictic 
on the pier,sides embedded in sand, the load bearing capacity of the pier in 
creases considerably as compared to that of a pier with shallow foundatio1 
Generally this additional load bearing capacity is not taken into consideratic 
in the design but is kept as a margin of safety. The previous railway practic 
was to work out the depth of bridge foundations according to the observatio 
made by Spring22 and Gales.23 In the case of shallow foundation, protect@ 


4 a cree Annual Report, Central Water and Power Research Station, Poona, In 
ia, -39, % 
19 “Behaviour and control of rivers and canals with the aid of models - part II,”Re 
eon Publication No, 13, Chapter VIII, Central Water and Power Research Sta., Poo: 
ia, a 
20 Technical Annual Report, Central Water and Power Research Station, Poona, Ind! 
“Stable Channels in Alluvium,” by G, Lacey, Journal of the Institution of Civil E 
gineers, Paper No. 4736, 1929, 7 
22 “River Training and Control of the Guide Bank System,” by F. J, E. Spring, Rail 
way Board’s Technical Publication No, 153, “= 
23 « Principles of River Training for Railway Bridges and their Application to the Cai 
of the Hardinge Bridge over the Lower Ganges at Sara,” y R. Gales, Journal of the In 


titution of Civil Engineers, December, 1938, Paper No. 5167, 
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s to be provided by stone pitching and if this is at high level, a lot of the 
one (stone used generally weighs 80 to 120 lb and is called one man stone) 

washed away downstream due to turbulence and has to be replenished, even 
ring floods, to ensure the safety of the bridge. Due to this turbulence very 
ep scour occurs downstream of the piers as in the case of Hardinge Bridge 
the river Ganges, depth of scour being of the order of 200 ft. 

The current railway practice is to provide a grip length for the pier equal 

half the depth of scoured bed below H.H.F.L. so that the total length of piers 
low H.H.F.L. is 3D( Lacey). 

Experiments25, 26 were carried out at the Central Water and Power Research 
ation, Poona, India, for testing the design of training works, waterway, and 
ngth of piers of arailway bridge at Mokameh on the river Ganges near Patna 
3ihar State). The pier foundation of this bridge has been taken to a depth of 
out 200 ft below H.H.F.L. which is equal to 3D(Lacey). As the foundations 
“e deep enough, protection by way of stone pitching round piers is not pro- 
ded. 

Various cases of bridges, for which rivers had to be trained to ensure the 
ifety of piers, have been experimented on. These experiments are described 
‘the Technical Annual Reports of the Central Water and Power Research Sta- 
on, Poona, India, for the year 1937-38, 1938-39, 1939-40, 1940-41, 1944 to 47, 
349 and 1952 to 1958. 

In the case of flash flood type rivers, fixing the waterway is much more dif- 
cult. In such cases the flood rises and falls so rapidly that the river has no 
me to scour its bed with the result that the afflux (difference in water sur- 
ce upstream and downstream of bridge) increases enormously; the bridge is 
kely tofail by outflanking. A railway bridge on Luni River in Rajasthan State2 
iled in this manner. 

In the case of inerodible bed material, it is difficult to estimate the maxi- 
um depth of scour. Hydraulic model experiments are unable to reproduce 
‘is scour due to obvious limitations. Recourse has, therefore, to be taken to 
eld data. In some cases, the maximum flood level.is underestimated and the 
aterway provided is insufficient. If the bed is inerodible afflux increases be- 
ond the safe limit, with the result that standing wave conditions prevail down- 
‘ream of the bridge, which lead to undermining of bridge piers. The railway 
ridge on Yeshwantpur River28 in Andhra Pradesh failed and collapsed for 
imilar reasons. 

The writer agrees with the author that afflux caused by a bridge depends so 
wach on the erodibility of the bed material. Thus in the case of the railway 
ridge on the River Ganges25 at Mokameh, the afflux caused by the bridge was 
nly a couple of inches for a flood of 18,000,000 cu sec. This was because the 
ood was sustained and the river scoured its bed as the discharge increased. 


24 Investigations carried out by means of models at the Khadakwasla Hydrodynamics 
esearch Station, near Poona in connection with the protection of the Hardinge Bridge 
hich spans the river Ganges near Paksey, by C.C. Inglis and D. V. Joglekar, East Ben- 
al Railway, Public Works Department, Bombay, India, 1936, Technical Paper No. 55, 

25 Technical Annual Report, Central Water and Power Research Sta., Poona, India, 


950. 

26 «Manual on River behaviour, control, and training,” Ch. VI, Pub, No. 60, Central 
oard of Irrigation and Power. ; 
27 Technical Annual Report, Central Water and Power Research Station, Poona, India, 


4, 
28 Technical Annual Report, Central Water and Power Research Station, Poona, India, 
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In the case of rivers with flash floods or with inerodible bed, full afflux has t 
be allowed for in the design according to standard formulas. 

As emphasized by the author, more field experiments are necessary to im 
prove our method of estimating scour at bridge piers. 


W. J. BAUER,29 M. ASCE.—The writer was in attendance at the University 
of Iowa, when the research described by Mr. Laursen was getting under way, 
and is, therefore, aware of its pioneering nature. The purpose of this discus- 
sion is to consider the application of the results presented by the author to 
a particular approach to waterway design. 

The fundamental purpose of a bridge over a stream is to separate the 
human traffic and the flowing water. It is commonly accepted that provision 
for the passage of the maximum conceivable discharge through the waterway 
beneath the roadway is economically undesirable. It is good practice to pro- 
vide for the safe, although infrequent, overtopping of the roadway by a rare 
flood. If the design is adequate, after the passage of such a flood, the high- 
way is immediately ready for service with only minor repairs being required 
during subsequent routine maintenance. 

In order to achieve this objective, the backwater produced by the struc- 
ture must be small at the stage of incipient overtopping. The writer has used 
values of 0.5 ft to 1.0 ft as being reasonably small for the backwater at this 
stage. The design problem then becomes how to provide for the passage of 
the flood corresponding to the stage of incipient overtopping without exceed- 
ing an allowable backwater, 

Some of the flow area required will exist beneath the elevation of the low- 
flow bed of the stream in accordance with the reasoning set forth by Mr. 
Laursen under the heading ‘Local Scour at Piers and Abutments: Backwater 
at Bridge Crossings.” The extent to which such area is available to the flow 
during flood may be controlled by suitable scour protection. 

At bridge sites in West Virginia, at which scour protection of broken rock 
is readily available, the writer has used a design velocity beneath the bridge 
of between two and three times the typical velocity in the stream during 
flood, This was accomplished without exceeding the allowable backwater at 
the stage of incipient overtopping. Such a relationship, between the velocity 
in the contracted section and that in the approach flow, gives sufficient as- 
surance that the material deposited over the broken rock fill during low flow 
will be scoured out during flood, provided the depth of scour required is not 
excessive, Fig. 6 indicates the relationship between depth of scour at an en- 
croaching abutment and a parameter dependent on the degree of contraction 
of the natural waterway by the structure, Fig. 6 indicates the attaining of 
depths of scour equal to the depth of flow at relatively minor amounts of 
contraction. It is, therefore, not difficult to imagine circumstances in which 
it would be possible to provide as much waterway area beneath the low-flow 
bed of the stream as above it. ¥ 

In the instances that the writer has in mind, the waterway area provided 
beneath the low-flow bed of the stream was about 20% of the total waterway 
area at the stage of incipient overtopping. The depth of the waterway area 
beneath the low-flow bed was only about 40% of the nominal depth of flow at 
the stage of incipient overtopping. Nevertheless, significant cost saving: 
were achieved compared to a longer bridge with equal waterway area entire 


29 Cons, Engr., Chicago, Ill, 
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ibove the low-flow bed of the stream. Itis possible to imagine circumstances 
fm which the length of the bridge might be out in half, the waterway area re- 
naining the same. This has obvious economic significance, 

The natural filling of the scour hole during the recession of the flood tends 
© minimize the depth of the pool of water that might be left standing under 
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j FIG, 12.—PROPOSED DESIGN 
) j ly made 
he bridge during low flow. As the beds of such streams are common 

p of pecans pools and bars anyway, the addition of another small pool 
eneath the bridge is of little consequence. Fig. 12 shows the type of design 
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L. J. TISON.29—The relationship proposed by Mr. Laursen for the prediction 
of scour at piers and abutments for the case in which sediment is supplied to 
the scour hole, in extremely interesting. 

The author presents an analysis in which he compares the flow without 
contraction and the flow due to a long contraction. He then uses, for both 
flows, an approximate form of the total sediment load relationship recently 
proposed by himself and consideres a bridge crossing as a long contraction, 
foreshortened in such an extreme that it has only a beginning and an end. 

The writer has found that a scour at piers and abutments also took place 
when the contraction of the flow was without significance and introduced the 
idea that the scour at piers may be affected by the curvature of the stream- 
lines around the piers. 

When a pier with an arbitrary shape is placed in a stream, it produces 
(Fig. 13) a first curvature with a center 0; and a radius pj, followed by a 
second curvature 02 Pg. In the neighborhood of the bank of the river, the pier 
exerts no action on the direction of the streamlines, cd. But, considering the 
region of the first curvature in the neighborhood of the surface, a line such 
AB, tangent in each point to the principal normals of the streamlines, the 
following relationship may be written: 


in which z is the vertical height, p is the pressure, v is the velocity and yis 
the specific weight of the liquid. The first supposition was that the motion is 
parallel to the bottom. 

In the neighborhood of the bottom, another relationship of the same nature 
gives: 
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B and B' are on a same line perpendicular on the bottom. 
The variation of the pressure between Band B' follows the hydrosatic law, 
so that: 


4 
p po! 4 
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From these three relationships 


with v greater than v' (v'is the velocity in the neighborhood of the bottom 


30 Prof,, Univ. of Ghent, Ghent, Belgium. 
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A consequence of Eq. 19 is that the motion cannot be parallel to the 
bottom, the trajectories must dive. 

This diving motion will have another consequence : a local attack of the 
bottom under the influence of the first curvature 01 91. Evidently, the im- 
portance of scour will depend on the value of the vertical component of the 


Bank 


FIG, 13 


diving motion, and this vertical component will depend on the value of the 
second member of Eq. 19. 

For example, small values of p will increase the depth of scour. Experi- 
nents with different models of piers, with the same length and breadth and 
with the same discharges and heights, gave a confirmation of this result. 


136 November, 1960 HY 9 


Rectangular piers (reduced values of P) produced a scour of 113 mm with 
a discharge of 301 per sec in a flume with a breadth of 70 cm and a height of 
10.5 cm, The length of the pier was 24 cm and its breadth was 6 cm. 

The simple rounding of the edges of the piers reduced the scour to 91 mm, 
whereas a triangular shape reduced it 70 mm. 

An aerodynamic shape gave no further reduction, but the shape of a lens 
produced a reduced scour of 54 mm. (large values of p). 

The protection realized by a single pile before the lens-shaped pier re- 
duced the scour to 38 mm, This pile produced bigger values of p and worked 
as a reduction of the width-length ratio. 

A change in the repartition of the velocity will have an influence on the 
secondmember of Eq. 19 and,therefore,on the value of the scour. A higher 
roughness of the bottom on a distance before the pier (with pebbles on the 
bottom) is realized. In Eq. 19, v' was reduced whereas v was increased, 
with the conclusion that the erosion had to be increased. That is what was 
found with an erosion of 71 mm with the lens-shaped pier (54 mm with the 
fine sand). 

The consideration of Eq. 19 shows that it would be possible to consider- 
ably reduce the erosion by adapting the variation of the radius p with that 
Orvis. 

Therefore, the repartition of the velocity from the bottom to the surface 
was measured and the lens-shaped pier was given a variable radii of cur- 
vature corresponding with the values of v2. The result was a non-prismatic 
pier with a larger base. The erosion around this new model was reduced to 
5 mm maximum, 

It is easy to see that the second curvature 09 (9 will give a relation 
(Eq. 19) with a negative second member, A rising motion will, therefore, - 
follow the first diving motion. This was observed on all the models. For the 
rectangular pier, the small radii of curvature around the upstream edges 
produced a quasi-vertical rising motion. : 

Behind the pier, with the lens-shaped section, the curvatures of the type 
2 was followed by another curvature of type i, and the formation of a sec- 
ondary smaller scour was to be observed just behind the pier. 

Many other results could be deduced from three first considerations and 
also from the consideration of the spiral motion induced by the diving motion 
deviated by the reaction of the bottom. Further results can be found in some 
of the writers publications, 

The same theory can be used for the study of the action of groynes and of 3 
the motion of sediments in derivations. : 
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Shag CHITALE.?!_In connection with estimation of scour at piers and abut- 
ments of bridges, Mr. Laursen has stated that comparison of model and pro- 
totype data indicated that the depth of scour could betreated as simply another 
length and that equilibrium depth of scour obtained in the field. This concept 
of equilibrium depth is due to the author and more clearly enunciated by the 
following quotation:2 


“At least as a first approximation the equilibrium scour depth, with 
certain qualifications as to the flow conditions, appears to be a function 
only of geometry, i.e. the relative depth of flow, the shape of the pier and 
the angle of attack. . . velocity of flow and sediment size (and, therefore, 
rate of transport and intensity of boundary shear) do not influence the 
equilibrium depth of scour... .” 


- Some investigations in models have been made on the subject of scour depth 
at bridge piers inIndiaat the Central Water and Power Research Station, Poona 
which the writer thought would be of interest in context with the author’s find- 
ings. 

The first series of experiments were conductedin connection with the Hard- 
inge Bridge over Ganga River and were reported in the Annual Reports of the 
Station for the years 1938 to 1942. A geometrically similar replica of Hardinge 
Bridge pier was embedded in Ganga sand of mean diameter of 0.29 mm ina 
parallel sided channel. The results of these experiments gave the relation 


ees 1 70 pws rede oee tat (20) 


b 


in which ds = maximum depth of scour below H.F.L.; d= depth of flow in the 
flume; b = widthof pier; and q, = discharge per foot run upstream of the piers. 
This relation is not dimensionally correct and cannot therefore be adopted for 
general application. 

Basic experiments were subsequently conducted in 1941 with the object of 
testing the influence of upstream depth andsand diameter onscour round piers. 
‘The pier tested in these experiments was also 1/65 scale model of that of the 
Hardinge Bridge. It was rectangular in section, of length 1 ft, width 0.6 ft and 
: semicircular cut and ease waters. The bed of the flume in these experiments 
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31 Chf, Research Officer, Flood Control Div., Central Water and Power Research 


Station, Poona, India. 
_ 2éseour Around Bridge Piers and Abutments,” by E. M. Laursen and A, Toch, Iowa 


‘Highway Research Board Bulletin No. 4, May, 1956. 
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FIG, 14,—PIER EXPERIMENTS 
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was laid with sand of 0.32 mm while the following materials were used just 
yround the pier in succession. 


Sieve number Type of sand Mean size, in mm 
1 screened dam sand m = 0.16 
2 White ‘V’ sand m = 0.24 
3 Nala sand m = 0.68 
4 Nala sand m= 1.51 


The sand round the piers was laid flush with the upstream bed level. A 
constant discharge of q = 1 cfs per ft was run and water level was adjusted to 


TABLE 2,—RESULTS OF EXPERIMENT 


Exper- Mean 
iment diameter 
No. of sand 

around 


pier, inmm 


0.1005 
0.1261 
0.1428 
0.1789 
0.2254 
0.3799 


0.1005 
0.1261 
0.1631 
0.2254 
0.3366 
0.4002 


0.1490 
0.1841 
0.2140 
0.3012 
0.4762 


0.1261 
0.1631 
0.2254 
0.3366 
0.4575 
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Pier experiments: In the 8 ft chamnel with sand of m = 0.32 mm 
Scale of pier: 1/65 b = width of pier = 0.57 ft = 37 ft 
; = length of pier = 0,925 ft = 63 ft 


with constant discharge q=1 cfs per ft 
ds = maximum depth of scour round the pier 


d = flow depth in the flume 

: 

get a particular depth, the depths varying from 0.5 to 1.45 ft. Each experiment 
was continued until final maximum scour was obtained round the pier. 

In afew tests in which the upstream depth was less than stable depth the up- 


stream bed scoured and blanketed the scour pit around the pier. ‘In such cases, 
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the maximum depth of scour at the nose was measured just before deposition 
in the scour hole of sand from upstream occurred. 

In experiments when sand around pier was coarser than the bed material 
upstream, the bed around the pier was laid higher than upstream level, to get 
scour round the pier for the upstream depth laid. 
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FIG, 15,—1/16 SCALE PIER EXPERIMENTS, ds : d FOR q=1 cfs 


Table 2 shows the results obtained important conclusions, as follows: _ 


1. With axial flow, maximum depth of scour was always at the nose of the 
pier, scour at sides being less by 5 to 15%. 
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2. The ratio of scour at the nose and depth of flow in the channel bears a 
imple relation with the approach velocity in the channel (see Fig. 14). 


3. The depth of flow on the upstream has also an influence on scour at the 
ier nose (see Fig. 15). 
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z FIG. 16.-RESULTS OF BASIC EXPERIMENTS 
i | 
It will -be seen that correlation of depthof scour either with upstream depth 
Fig. 15) does not appear to be satisfactory. The writer, therefore, further 
nalyzed the experimental data and found that the Froud number provides a 
etter criterion. Fig. 16 shows plot of depth of scour against the Froud num- 
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ber, Although some scatter of points is evident in this plot, the general trend 
is remarkable, statistical equation obtained over the range of experimental 
data being 


ti OSbs 6.08 eh 40 xO ol a a ee (21) 


It is thus seen that experiments at the Research Station do not lend support 
to the author’s equilibrium scour theory, and it appears that further investiga- 
tions both in the field and in the laboratories are desirable before it can be 
accepted unreservedly and with confidence. 

Grateful acknowledgment is made of the kind permission given by M. G. 
Hiranandani, Director, Central Water and Power Research Station to refer to 
experiments previously conducted at the Station and also for useful suggestions 
offered by him. 


A. RYLANDS THOMAS,?2 F, ASCE.—In Fig. 17, the author’s design curve 
for piers (Fig. 9), with a factor of 0.9 for application to piers with semicir- 
cular nose form, is compared withresults of experiments with models of piers 
carried out by the writer in association with Sir Claude Inglis.33,34 These 
piers were models of the piers of the Hardinge Bridge over the Ganges River 
(now in East Pakistan) which were 37 ft wide and 63 ft long, including the semi- 
circular nose and tail. The modelpiers were fixed ina parallel-sided channel 
with abed of incoherent Ganges sand about 0.3 mm mean diameter and the pro- 
cedure was to run a constant discharge without sediment load until scour had 
ceased. 

It will be seen from Fig. 17 that the author’s design curve, and still more 
his curve from the Iowa data, Fig. 3, lie well above the Poona results. The 
difference appears to be too great to be explained by sediment load, 21 though 
this cannot be ruled out as the upstream depth may be reduced more by sedi- 
ment load than is the depth of scour at the pier. It would be of value if the 
author presented the data on which he based his design curve, showing partic- 
ularly the effect of load. 

It is also most desirable to compare small-scale results with observations 
made under full-scale conditions. Such data are difficult to obtain because of 
the need to take observations during high floods which are very often of short 
duration. It is easier to measure depth of scour at the pier than in the chan- 
nel upstream, where a number of observations must be taken to average oul 
the effect of bed waves. Even when this is done it is not certain that the bec 
was in equilibrium at the time of observation. The Poona experiments showec 
thata reduction in channel depth, due for example to resistence to scour, woul¢ 
increase the depth of scour at the nose of the pier. q 

The relationship between depth of scour at the nose of the pier and the depth 
of channel upstream is, therefore, perhaps not the most suitable one for com- 
parison of full-scale data. Nor is it generally the most suitable for practical 
application of a design formula, because the upstream depth during a maximum 
flood is not often known beforehandand would have to be calculated. An erro 
in assessing this depth would lead toa corresponding and greater error ines 
mating the level to which scour is liable to occur. ‘‘ 


32 Cons, Engr., London, England. 3 


33 Annual Reports (Technical), Central Irrigation and Hydrodynamic Research Stati 
Poona, India, 1938-39, p. 39; 1939-40, p. 33; 1940-41, p. 35, 

34 “The Behaviour and Control of Rivers and Canals,” by C. C. Inglis, Central Wate1 
power Irrigation and Navigation Research Station, Poona, India, 1949, p. 327, 
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A more basic relationship is that between depth of scour below water level, 
)s, Width of pier, b, and upstream discharge per unit width, q. The relation- 
hip indicated by the Poona results is 


D 2/3 0.78 
_ = 1.70 (3) Beene Vomeer mas tee (22) 


n ft-sec units, applicable to scour infine sand within the limits of q2/3/b or 
o/b at least from 1 to 6 (see Fig. 17) and not greatly affected by the grade of 
and. It is necessary to take into account the effect of obliquity of approach 
low, which may be done by using for b the projected width of pier ona plane 
ormal to the direction of approach. A factor of safety should also be used in 
esign. 


Design curve, Bulletin 4 
Semi circular cutwater, axial approach) 


Poona Models 
Semi-circular cutwaters 
axial approach 
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scale 
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FIG, 17,—-SCOUR AT BRIDGE PIERS 


The discharge per unit width to be used for design must, in most cases, be 
‘reater than the mean value because in rivers the deep channel meanders from 
ne bank to the other. It is therefore advisable to use, at least as a check, the 
lationship proposed by Inglis34 giving the maximum depth of scour around 
pier as approximately twice the Lacey regime depth, that is, 


Der Gt. > wd scie penis pe ONE CES: 


n which Dg is the depth of scour in feet below water surface, Q ig the dis- 
harge in cu sec, and f is the Lacey silt factor (approximately 1.8 Vm where 
n = mean diameter of bed sand in mm). Eq. 23, derived from data of several 
ridges in India and Pakistan, takes into account obliquity of approach and con- 
| ntrationof discharge but not the width of pier, which is clearly an important 
actor. The data related tonormal widths of pier and to depths of scour ranging 


rom 25 ft to ALT tt. 


144 November, 1960 HY | 


MUSHTAQ AHMAD.35_The problem of correct estimation of depth, shape 
and extent of localized scour is important from the point of view of establishing 
sound design practice and safety of hydraulic structures such as bridge piers 
abutments, spur dikes, groins, or pitched islands. This problem is much mort 
important for hydraulic structures in alluvial rivers of West Pakistan wher 
fine sand is found for hundreds of feet in depth, and scour depths of 40 ft to 8 
ft below water level are common. The author has, therefore, dealt with a sub. 
ject of special importance and great utility to this region. 

It is proposed to discuss the author’s approach to the problem with specia 
reference tofield andlaboratory experience in West Pakistan. The author ha: 
assumed the depth of scour as another length and has related it to the norma 
depth or the width of a pier, and maintains that the depth of scour does not de. 
pend on the degree of concentration until scour holes around neighboring ob: 
structions overlap. He maintains that there is equilibrium or limiting dept 
and believes that for a given mode of sediment movement, the depth of scou: 
depends only on the geometry of contraction and the approach depth. He hold: 
that the effect of velocity and sediment size can be neglected as that of secondar: 
order. The writer agrees that there is an equilibrium and limiting depth o 
scour and that the effect of sediment size can be neglected. Tests made oi 
depth of localized scour at spur dikes36 also showed that localized scour dept! 
does not vary with grain size in the range (0.1 m to 0.7 m) usually met witl 
in alluvial plains of West Pakistan. However, this concept may not be vali 
for the entire range of bed material sizes ranging from fine sand to gravel an 
boulders, 

The author’s view that scour depth is another length connected with norma 
depth or width of pier implies that it varies with the discharge, or more cor. 
rectly with the discharge per foot run. The author’s Eq. 8b and c holding fo: 
rivers and streams in West Pakistan are: 
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)which K may be a function of boundary geometry of the contraction at the 
ridge, abutment shape and thickness and shape of the nose of the piers, and 
) forth. Because the depth is a dependent variable on discharge intensity, the 
itter may beused in preference to depth. Here only lies the difference in ap- 
roach between the author and the writer in studying this problem. As shown 
reviously the two approaches are not very different. The writer prefers to 
tudy the variation of Dg/a2 3 = K as a function of boundary geometry, shape 
f pier nose, and abutment, characteristics of bed material, and distribution of 
elocity in the cross section at the piers representing the concentration of 
ow. The functional relation for study may be of the type: 


D 1/2 JV 
ie> Vv gys B 
a2 73 ( eae | ae aes ’ S, B, = =) (4 108 6 TRU SCAES Poses (28) 


1 which @ is the angle between the current and the spur dike, an abutment or 
pier, vi/ 2/V is the ratio of mean velocity in half the channel width on the 
ide of spur dike or abutment to the mean velocity in full section. It will bea 
1easure of concentration of flow; W is fall velocity of bed material; B is the 
iver width upstream of the bridge; B; is distance between the abutments; n 
3 the number of piers,; and t is thickness of each piers. 

The writer does not agreethat the depthof scour does not dependon the de- 
ree of concentration. It has been found that by changing the concentration or 
slocity distribution at bridge site bytraining works upstream, the scour depth 
an be considerably reduced. The variationof scour depthas aresult of change 
1 velocity and flow distribution due to the training works is explained in Figs, 
B, 19, and 20. Figs. 19 and 20 depict a model of the Ling Stream showing a 
ridge below a curve. The concentration of flow near the left flank is notice- 
ple in Fig. 19. A pitched island and a spur on the left bank fans out the flow 
- the bridge site to obtain more uniform velocity and lesser scour as shown 
1 Fig. 20. 

In the method of plot adopted by the author in which scour depth in relation 
) U/S depth has been studied, the non-dimensional scour perimeter may not 
ave significant relation with velocity, for, in open channels as the velocity or 
ischarge increases so does the depth, although not necessarily by erosion of 
1e bed. For a uniform approach velocity, a mean depth can be assumed for 
3e in the author’s relations. But, generally, abnormal scour depths at piers 
r abutments are associated with concentration of flow resulting from the var- 
ition of velocity and depth in the approach section and the selection of repre- 
sntative approach depth for use in the author’s relation and the estimation of 
ischarge distribution in the channel and on the berm becomes very difficult. 
1fact, the author’s term (Qc + Qo)/Q¢ in Eq. 9 is, also, a measure of con- 
sntration of flow. However, the estimation of Qe, Qo, or Qw is difficult for 
ractical use in the computation of maximum scour depth in case of a curved 
sproach. The writer has shown9 that change in flow concentration above a 
pur dike, can be depicted in terms of velocity distribution and common types 
re depicted in Fig. 21 by I, Il and I. The type Ill velocity distribution gives 
1e maximum concentration onthe outside of abend as due to negative velocity 
nthe inside, a part of the water way is blocked by reverse flow. For the three 
mpes of velocity distribution depending on different types of approaches, the 
lues of V,/V and K as determined from writer’s studies onspur dikes (which 


ay be applicable to abutments), are given in Table 3. Since maximum scour 
an occur on any one or a group of piers in a meandering alluvial channel, the 
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jepth of foundation of piers has to be computed from the values of K as fixed 
yy the river curvatures likely to occur, keeping in view the restraints imposed 
mm river by the upstream training works. The method commonly used in West 
Pakistan supported by Laboratory studies and field observations consists in 


working out q = oak 


TABLE 3.—VALUES OF vi /V AND K 


Approach condition 


(1) 


Scour below a severe bend on the concave side ac- 


companied by a swirl on the convex bend. 2 = 2525 
Moderate Bends. 1.5 -1,75 
Straight obstruction placed at any angle of 30° to 90° 1.2 -1.5 
Straight obstruction placed at an angle of 90° - 150° 

to the flow 1.5 - 1.75 


TABLE 4.—MAXIMUM SCOUR DEPTHS OBSERVED AT BRIDGE PIER AND 
ABUTMENTS ON DIFFERENT MODEL STUDIES AT IRRIGATION RESEARCH 
INSTITUTE, LAHORE.) 


River Site Scour depth K Thic’ y.la Dg 
Gare vite Ds ness| —° . = ds 
Oy SNe Nc of b + Y 
rect- q : ° 
ae pier 
scour (b) 
depth, 
(1) (2) (7) (8) | (9) | (10) | (11)}(22)) (13) 
Ravi Shahdara 64.0 1.68] 26.4] 10 |2,64) 2.06) 47 
Bridge 
Jhelum | Jhelum 39 1,60 | 23 10 |2,.3 |1.98 42.8 
bridge® 
Maagla 90 TGONs28 oF (AUCH Nec ch ieee || ote 
New bridge ment 
Deg S-harakpur 45 1.384} 21.5] 6.25 |3.45]2,22| 35.5 
Rohi Kasur 39 1.4 |14 |6,25 |2.25/1.97) 26.3 
Sohan Dhok Pathan 42 1.36} 18.6|9.0 |2.07)1.94 36,1 
Indus Thatha 1,100,000 68 1.6 |25 | 7.0 |3,.57)2,25| 40.8 
| Sujjawal 
a proposed 
bridge. 


a Computation from Fig. 9, Bulletin No. 4. 
ig Bed material - five sand. 


_ At the bridge site, the value of K canbe selected from the maximum curva- 
re likely tooccur. This will give scour depth below water level, and the depths 
elow bed level can then be worked out. Maximum depth of scour from model 
udies at bridge piers and abutments of different rivers and streams, after 
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correcting the results for the effect of scale distortion on geometry of scour, 37 
are summarized in Table 4, In this table, scour depth has been computed from 
the author’s design curves of Bulletin No, 4, in Fig. 9. The scour depths ob- 
tained after correcting for the model scale distortion are generally higher than 
those of the author. In computing Dg value by his method, the actual mean 
value of Y, has been taken from the cross section above the bridge. The value 
of K does not exceed 1.7 because severe bends are not possible in these cases 
due to the presence of guide banks. 


TABLE 5,—SHAHDARA RAILWAY BRIDGE ON RIVER RAVI* 


Q H.F.L.| Gauge 


¢ A spur constructed u/s to current approach. f 
aa 


Fig. 22 shows abnormal scour at a guide bank head. A heavy embayment is 

noticeable on the right and construction due to silting caused by roller on the 

left. The K value of 2.3 can be obtained in such a case. i 
The heavy embayment resulting in abnormal scours illustrated in Fig. 2 


with K value greater than 2 are possible only on spur dikes or guide bank head 


37 “Effect of Scale Distortion, Size of Model Bed Material and Time S$ 20- 

’ cale on Geo- 
metrical Similarity of Localized Scour,” by Mushtaq Ahmad, Proceedings, Internatl, 
Assoc, of Hydr, Research, The Hague, Netherlands, 1955, 4 


a 


—— 
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und are not to be allowed on bridge piers or abutments, In fact, for alluvial 
rivers or streams of West Pakistan, guide banks at least equal to the length 
af the bridge, with curved heads and expanding water way on the upstream side, 
are provided to shift the maximum embayment and abnormal localized scour 
0 the guide bank heads instead of allowing it to occur near the bridge with 
short abutments. Under these conditions, the values of K for estimating the 
maximum probable depth of scour can safely be taken between 1.7 and 2.0. In 
case of a shortabutment length before the conditions for abnormal scour where 
K greater than 2 are obtained, the approachroads or railway will be threatened 
and cut by the embayment formed by an alluvial meandering river. It is there- 
for necessary to design pier and abutment depths for scour calculated from 
K 1.7 to 2.0 and provide proper guide bank to keep the road or railway ap- 
proaches safe from embayment and to keep abnormal scour away from the 
main bridge crossing. Actual scour depths observed on the railway bridge on 
Ravi near Lahore and the computed values of K as recorded in Table 5 also 
supported the recorded values of K. 


PIER LUIGI ROMITA.°8—The author’s long lasting efforts to cast light upon 
the phenomenon of scour around bridge piers and abutments are to be greatly 
commended because of the great practical importance of the problem, The 
clear and condensed presentation in this paper of the results of these efforts 
is a substantial step towards a generalized solution of the problem, and will be 
of great use to designers and tothe administrations responsible for the con- 
struction and maintenance of roads and railways. To this end a coordinated 
effort should be made, in order to insure the necessary verification on proto- 
types of the proposed relationships. 

There is no doubt that one of the most important factors of scour, when 
piers have a sufficiently high length per width ratio, is the angle of attack be- 
tween the pier and the flow. Even small deviations of this angle from the zero 
value are responsible for fast increases of the scour; this usually represents 
a much greater danger to the stability of the pier than any underevaluation of 
the maximum possible flood, The collapse of many piers during floods is, in 
fact, to be ascribed more to the angles of attack, which may occur due to un- 
usual cross-currents and deformations of the river bed, than to the unexpected 
entity of the discharge. 

In view of all this, the writer carried out, some years ago, a systematic 
model investigation of the influence of the angle of attack upon the depth of 
scour. The experiments were carried out ina glass-flume of the Hydraulic 
Laboratory of the Polytechnic Institute of Milan, supplied with clear water, 
using uniform sand of 1 mm diam as bed material, and in such conditions that 
there was no general bed movement but only localized scour around the pier. 
The shape of the pier was not particularly studied, and it reproduced a rather 
widely used type of pier; its length per width ratio was about 5, Two series of 
tests were carried out with different values of the discharge, while the angle 
of attack was varied from 0° to 90°. 

The data obtained at a zero angle of attack are in rather good agreement 
with the curve representing Eq. 11 Fig. 9, but for the fact that they consistently 
lie a little below this curve. An increase by 60% in the discharge (which cor- 
responded to an increase by 20% in the average velocity of flow) brought around 
only a 20% increase in the maximum depth of scour, at the same angle of attack. 


38 Assist. Prof., Hydr. Structures, Politecnico di Milano, Milan, Italy. 
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The tests with varying angle of attack have shown the basic importance of 
this factor. For an angle of 15° the maximum depth of scour was 1.9 times 
that corresponding to a zero angle; for an angle of 30° the increase of scour 
depth was in the ratio 2.6:1 in respect to the zero angle, for an angle of 45° in 
the ratio 3.0:1,and for an angle of 60° in the ratio 3.3:1. Increases of the angle 
of attack beyond 60° and up to 90° did not bring around any further appreciable 
increase in the scour depth. These values are considerably in excess of those 
indicated by Fig. 10. 

Another interesting observation was that the point of maximum scour depth 
always occurs very closeto the pier wall, so that, if the depth of the foundation 
is not sufficient, the pier will easily be undermined, and collapse. 

As a conclusion, the fundamental importance of avoiding any angle of attack 
between the pier and the flow shouldagain be stressed. In order to obtain this, 
the river banks should be stabilized with adequate measures for a sufficiently 
long stretch upstream of every bridge crossing. In braided rivers, however, 
the possibility that an angle of attack occurs should always be taken into ac- 
count, and the pier foundations designed accordingly. 
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CONSERVANCY DISTRICTS AS FLOOD CONTROL ORGANIZATIONS# 
Discussion by Samuel A. Greeley 


SAMUEL A. GREELEY,” F. ASCE.—This is an excellent and timely paper. 
In my opinion, it has important bearing not only on conservancy districts, but 
also on other forms of municipal organizations such as authorities, sanitary 
districts, water districts, and, in Michigan, refuse disposal districts. The im- 
portance of organizations of this kind, which may be generally referred to as 
Authorities, was emphasized in finefashion by Mr. Robert Moses in his address 
before the American Public Works Congress in New York, on August 17, 1960. 
It is pertinent, therefore, to consider Mr. Chambers’ paper as applying to more 
kinds of public works activities than flood control. 

It is of interest to refer to the Progress Report dated January, 1945 of a 
Committee of the Society on Organization, Financing, and Administration of 
Sanitary Districts, of which the writer was Chairman. Relative to the organ- 
ization and creation of sanitary districts (or authorities) , the following is ab- 
stracted: 


In general, the Committee feels that the steps required tocreate, form, 
and organize a sanitary district should be taken within the locality which 
will comprise the district. 

The Committee favors the origination and creationof districts by local 
petition followed by ageneral vote by qualified persons withinthe desig- 
nated boundaries. 

The Committee favors the determination of the extent of the district 
by a competent engineer or board of engineers after an adequate survey 
and investigation of the locality. 

The Committee favors, in general, the election locally of the members 
of the administrative board. The Committee recognizes, however, that 
many appointed boards have been eminently successful and suggests that 
the appointment of board members by some locally elected official may 
be an acceptable method. 


In Item 3 of the Summary, Mr. Chambers refers to the needs of “such de- 
velopments.” The writer suggests that not only the need, but also the relative 
need as compared to other interests should be included in the original studies. 
It is further suggested that the paper does not give sufficient weight to fi- 
nancing procedures. In two elections under the Illinois Conservation District 
Act, organizations of the districts have failed. The writer thinks one reason 
was that the enabling act limits all financing to taxes derived from assessed 
valuations. It seems likely that properties in different parts of the district 
will have different views relative to the need and validity of the proposed dis- 
‘trict or authority to be organized and financed. The obtaining of funds by the 


@ April, 1960, by Cloyde C. Chambers. 
5 Partner, Greeley and Hansen, Engrs., Chicago, Ill. 
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Miami Conservancy District which was organized shortly after 1913 provided 
for charges or assessments by zones carefully related tothe resulting benefits 
or uses of the flood control works. This is a good example of an approach to 
fair financing. Having in mind the conception that the paper as a whole is use- 
ful in the consideration of other kinds of municipal organizations such as au- 
thorities and districts, reference should be made to the so-called two-part 
rate as developed by a Joint Committee of the Society and the American Bar 
Association with six other National organizations. The report of this group 
stated the following fundamental principle relativeto afair procedure in finan- 
cing water, sewage, and other public works: 


“The needed total annual revenue of a water or sewage works shall be 
contributed by users and non-users (or by users and properties) for 
whose use, need and benefit the facilities of the works are provided ap- 
proximately in proportion to the cost of providing the sue and the bene- 
fits of the works.” 


This statement means that the needed total annual revenue (operating and 
debt service costs) shall be derived from the users and properties in propor- 
tion to the amount which each causes to be spent. 


» 
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COMPARISON OF STREAM VELOCITY METERS* 


Discussion by Harold G. Golden and I. L. Trotter, Jr. 


HAROLD G. GOLDEN, 1° A. M. ASCE and I. L. TROTTER, JR.14—1n this 
most interesting paper, the authors concluded that the screw type and cup type 
current meters gave identical results in U.S. Lake Survey flow measurements. 


TABLE 4.—COMPARISON OF DISCHARGE MEASUREMENTS 
OBTAINED BY PRICE AND OTT METERS 


ee ae es 


Measured discharge, in Percentage varia- 
cubie feet per second tion of discharge 
obtained by Ott 
Date meter from 
; Ott meter discharge obtain- 
be eg (adjusted) mahi 
Price meter 
(1) (2) (3) (4) 
3-28-58 700,000 692,000 1.14 
4-16-58 796,000 774,000 -2.76 
4-28-58 795,000 796,000 +0.13 
5-23-58 1,178,000 1,184,000 +0.51 
6-— 9-58 475,000 474,000 -0.21 
6-30-58 688,000 688,000 0 
7-21-58 667,000 672,000 +0.75 
8- 4-58 1,032,000 1,005,000 -2,62 
8-25-58 624,000 624,000 0 
9- 8-58 335,000 326,000 -2.69 
9-26-58 399,000 397,000 -0.50 
10-20-58 291,000 291,000 0 
11-10-58 206,000 205,000 -0,49: 
Average percentage variation of 1958 measurements -0.69 


662,000 672,000 +1,51 

1,047,000 1,059,000 +115 

5-11-60 746,000 751,000 +0.67 
6- 3-60 755,000 763,000 +1.06 
6-17-60 546,000 547,000 +0.18 
7- 6-60 720,000 731,000 +1.53 


‘Average percentage variation of ‘1960 measurements 
= 

Average percentage variation of all measurements -0.15 
—E 
"8 April, 1960, by F. Wayne Townsend and F. A. Blust, 

_ 10 Hydr, Engr., U. S. Geol, Survey, Jackson, Miss. 

_ 11 Engr. Tech., Hydr., U. S. Geol. Survey, Jackson, Miss. 
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Of particular interest was the comparison between the Price (cup-type) meter 
and the Ott (screw-type) meter. 

The purpose of this discussion is to present additional data that indicate 
there are no Significant differences between independent discharge determin- 
ations using Price and Ott meters. During 1958 and 1960 the U. S. Geological 
Survey obtained nineteen pairs of discharge measurements using the Price and 
Ott current meters at the stream-gaging station on the Mississippi River near 
Vicksburg, Miss. The results of these measurements are shown in Table 4. 
All measurements were made from the downstream side of the highway bridge 
near Vicksburg. In all cases, the Ott meter discharge measurements were 
made after the completion of the Price meter discharge measurements, with 
an average time lag of 44 hr. The cross sectionalareas used forthe Ott meter 
measurements were based on the soundings obtained during the Price meter 
measurements. For comparative purposes, The Ott meter discharge meas- 
urements were adjusted when a significant change in river stage occurred be- 
tween the two measurements. 

An equipment change was made between the measurements obtained in 1958 
and in 1960. The measurements in 1958 were obtained using a 300 lb Vicks- 
burg-type lead sounding weight and a single pin attaching the Ott meter to the 
hanging-bar. The measurements in 1960 were obtained a 150 lb Columbus- 
type lead sounding weight and two pins attaching the Ott meter to the hanging 
bar. The effect of this equipment change is not evaluated butis notedas a mat- 
ter of record. 
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FRICTION LOSSES IN LINES WITH SERVICE CONNECTIONS#@ 


Discussion by M. H. Diskin, Annabel L. Tong, M. B. McPherson 
and J. V. Radziul, and K. M. Yao 


M. H. DISKIN.°—The author has presented a solution, based on the Hazen- 
Williams formula, for the problem of head losses in a pipe with a non-constant 
discharge along its length. The solution is based on a number of assumptions 
that should be critically examined before the solution can be adopted for gen- 
eral use. The main assumptions on which the analysis is based are: 


1. The head losses per unit length of pipe, at any given discharge, are the 
same whether it is a constant discharge or whether it is variable along the 
pipe. 

2. The local head losses inthe main stream due to the side outlets are 
negligible. 

3. There is no redistribution of energy content, per unit weight of fluid, 
between the main flow and the side flow. 


The first two assumptions are mentioned or implied in the orginal paper 
and require no further explanations. The third assumption refers to the pos- 
sibility that water arriving at a side outlet (Fig. 8) with a total head (energy 
per unit weight) Hy and discharge Q, will be divided into two streams Q2 and 
Q3 in such a way that one of them will have a higher head (He 2 Hy) and the 
other a lower head (Hg < Hj) than the original head (Hy). The total energy 
balance will, of course, still be in accordance withthe principle of conservation 


of energy: 
vy Q1 Hy > vy Q2 Hp + v7 Q3 HZ «eee eee eee eee (21) 


in which y is the unit weight of the liquid. 

That sucha redistribution is possible can be demonstrated witha simplified 
analysis of the conditions at a side outlet using the momentum equation, It is 
assumed for simplicity that the side outlet is ina horizontal plane. The free 
pody to which the momentum equation will be applied is defined by sections 
1-1, 2-2, and 3-3 in Fig. 8. Section 1-1 is taken just before the side outlet, 
and the other two sections are taken sufficiently downstream for the velocities 
to be approximately uniform. The momentum equations for the free body pre- 
iously defined, in directions parallel and perpendicular to the main pipe axis, 
and based on the usual assumptions for the momentum equation, are: 


Py Ay - Pe Ay=5 (Qy-Vg-Q1V1) cee eee ee (22) 


 @ April, 1960, by D. L. Muss. 
‘ 5 Acting Head of Hydr. Lab., Technion, Israel Inst. of Tech., Haifa, Israel, 


He 
; 


fh 


158 November, 1960 HY@ 


and 
pg Ag - P3 Ag = & (Q3 V3 - 0) De OE be (23 


in which Ay and Ag are cross sectional areas of main pipe and side outlet re 
spectively, py, pg, pg are pressures; Vy, V2, V3 are velocities at the corre 
sponding sections and pq is the pressure at the main pipe wall opposite the sid 
outlet. 


oe ® 


FIG, 8 


Dividing Eq. 22 by y Ay, it is transformed to: 


PI = P2 = Vase e, vie (2 : 

=F % ‘ RTL msctittens pay ko athe 4 

or 9 P 
2 3 

Sue Tt let Betis (25 

y g y g eee eve . ° q 

which may also be written as: 7 
2 2 r 

Hy + uf = Va ; 

1 2 g = Ho ot 2 g © 10 “es @ 0 16 ©. a) ee © 638-6) (26 


Since the velocity at section 2-2 is smaller than that in section 1-1 (Vo < vy 
the total head at section 2-2 is higher than the total head at section 1-1: 


4 
i Mate 8 ae tra arte trade prt (27 
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Dividing the second momentum Eq. 23 by y Ag and assuming that the pres- 
ure at the pipe wall opposite the side outlet (p4) is approximately equal to 
he inlet pressure (Py) 


P4 xs Py Rivera tha> pit, Pall wsttusantcitenb ita Satie libs (28) 
he equation is transformed into: 
2 
PplyF Pye Mg 
mipitntag -08 Desi Sirs eoriye apn tsnieie gotiie (29) 
Y g 
r 
Vv 2 
Bibel wo (30) 
< - A a ea ree ee 
hat may also be written as: 
2 2 
i ee we 
Sa Wallner ee hay SE wie “eteniareruieda Ss reniae ete (31) 
rom which it is seen that: 
EES), Seen Irae tae earesar | er (32) 


Jt is not claimed that the preceding analysis is exact; it was put forward 
rather to show that a redistribution of energy at a side outlet is possible and 
hat it could influence the results. Evidence that such a distribution takes place 
san be seen in the experimental results reported by J aeger® for head loss co- 
sfficients at side outlets that, at certain combinations of flows, assume negative 
values indicating a gain of head. 

Consideration of the possibility of the energy redistribution, together with 
jossible errors introduced by the other assumptions given previously, leads to 
-he conclusion that more experimental work is desired on the subject of head 
losses in pipes with varying discharge, before the results presented in the pa- 
ser can be accepted as describing the actual conditions in such pipes. 


ANNABEL L. TONG,” M. ASCE.—The author’s methods in determining the 
‘riction losses in pipelines with service connections will give exact solutions 
‘or single pipelines of uniform size and constant coefficient of friction in the 
pipe. In Tables 2 and 3, equivalent discharges Q, are given for ideal cases of 
inifinite and finite numbers of connections versus a/Qt) andalsoversus num- 
ver of connections “n” for the case of finite number of connections. But, in 
actual practice, water network systems in service are mostly constructed of 
snclosed loops and sometimes combined with a number of dead end lines. It 
would be interesting to know to what extent the author’s methods are applicable 
to a network system. 

In order to make use of his tables, one may convert a single loop system 
consisting of one source of supply into a single line by taking an equivalent 
length Le to represent the complete system as a unit. When the pipelines in 
one system were built at the same period in an area of similar environment, 
their physical condition would be assumed alike and a constant c value would 


«6 «Engineering Fluid Mechanics,” by C. Jeager, Blackie &Son Ltd., London, England, 


s 


1956, pp. 477-480. 
_ 7 Hydr, and Structural Engr., Stearns and Wheler, Cazenovia, N. Y. 
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be used forthe entire system with sufficiently good accuracy. With these sim 
plifications, the author’s tables can be readily applied to a simple loop syster 
without much computation. 

An example is shown in Fig. 9 to illustrate the degree of accuracy of hi 
methods applied to the loop systems. A water distribution network of 12 pipe 
lines, forming four enclosed loops is shown and the size and lengthof each lin 
are marked in parenthesis adjacent to the line. A frictional coefficient C c 
100 is assumed for the entire system. Seven service connections and the cor 
responding discharge rates are indicated in junction points B, C, D, E, G, B 
and I. If a hydrant test flow of.300 gpm (Q1) was make at junction F, and th 


(Q;), = 1000 
(Q;),= 1700 


(12 in.,2.000 ft) (12 in.,2.000 ft) 


(12 in, 1500 ft) 


(8 in,, 1750 ft) 


106 


(12 in., 2000 ft) 


(8 in, 


100 


(12 in. 1500 ft) 


(8 in., 1 250 ft) 


(10 in., 1000 ft) 


(10 in., 2.000 ft) (10 in, 2000 ft) 


100 


100 


FIG, 9.—WATER DISTRIBUTION NETWORK 


head loss (Hy) through the system during this flow is measured to be 4 f 
what will be the corresponding head loss (He) when the hydrant flow increase 
to 1000 gpm (Qe) ? ; 
With the information obtained from the hydrant test the head loss (Ha) ma 
be computed as follows: ¥ 
Method 1.—The system is assumed to be a single line of infinite number ¢ 
connections where q = 700 gpm and (Q) 1 = 1000 gpm. From Table 2, wit 
q/Q = 0.7 and Qe/Q = 0.676, one may find: (Qe) 1 = 676 gpm. To determii 
Le of the system, the Hazen-Williams formula ina simplified form, may b 
4 


used. } 
8 « Nomograph For Equivalent Lengths of Water Pipelines,” by Annabel L. Tong, Wa 


ter and Sewage Works, February, 1960. 
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fy WAS Sameera? We te), Re ee LE (33) 


119 x 10° 


here Le, is the length (in feet) of the equivalent pipe for the network. The 
nit for Qe is gpm and the unit for H is linear feet. 
‘hen 


_ 119x10° x4 
(676) 1-85 


Then Q changes to 1700 gpm and q/Q to 0.41 from Table 2, (Qe) 2 is com- 
uted as 0.087 x 1700 = 1370 gpm, and 


_ 280 x (1370) 1-85 
119 x 10° 


Method 2.—If seven connections are assumed inthe equivalent pipe, Table 3 
an be used to compute Qo. Using n = 7, q/Qt = 0.7 then (Qe)1 = .734 x 1000 
734 gpm. The equivalent length is computed as 


_ 119 x 10° x4 
(734) 1.85 
imiliarly (Qe) = -843 x 1700 gpm = 1430 gpm, and 


= 280 ft 


H = 14.8 ft 


Ee = 240 ft 


1.85 
_ 240 x (1430) Saat 
119 x 10° 


Head loss Hg is then computed by the Hardy Cross Distribution Method, the 
nswer is found to be 14.0 ft. Using this value as a basis for comparison, the 
ercentages of error from the preceding methods are approximately equal to 
6% and -2%. If the number of service connections reduces to 3 by discharging 
00 gpm at junctions H and D and 300 gpm at junction G only, the result by 
fethod 1 is unchanged and by Method 2, Hg equals to 13.0 ft. Again using the 
lardy Cross Method, the head loss Hg corresponding to 3 service connections 
3 found to be 15 ft. The degrees of accuracy by Methods 1 and 2 are -2% and 
15%. A change of location and rate of service discharges shall effect the value 
f the equivalent length for the system slightly, and subsequently shall effect 
he head loss under the same fire flow also. But both the methods are satis- 
actory to use for practical design problems. 

For a complex system consisting of pumping facilities and two or more 
ources of water supply, the assumption of using an equivalent length to rep- 
sent the system will not be correct, therefore, the above analysis should be 


imited to the simple networks only. 


Hg 


—M.B. McPHERSON,? M. ASCE AND J. V. RADZIUL.!°—Before commenc-~ 
ing a discussion on the analyses presented in the paper, consideration of some 
verall design factors is necessary to place the limited scope of the subject 
n perspective. 

“In United States metropolitan municipal water works, it is common practice 
Co) supply the distribution system, or separated service districts, via an arte- 


"9 prof, of Hydr. Engrg., Civ. Engrg. Dept., Univ. of Illinois, Urbana, Tl. 
10 Water Distribution Design Engr., Water Dept., Philadelphia, Pa. ~ 
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rial (or feeder) network, The smaller-pipe grid network is connected to the 
arterial system at relatively few points (about 1000-1500 ft apart). In a strong 
system a greater number of grid connections is not necessary while a mini- 
mum number is desireable to facilitate operations. Service connections into 
any arterial main grossly complicate operation, and are seldom tolerated. 
Service connections are normally made into grid mains. In small systems or 
districts, the arterial system and grid system are often inseparable or indis- 
tinguishable, Analysis of large systems is normally concentrated on the arte- 
rial network. Some of the major design considerations for arterial network 
analysis have been presented previously, 11,12 The author’s paper is con- 
cerned principally with grid “water mains.” In closing, the author states that 
his procedures are “particularly suitable for use in network analysis... .” 
While some of his procedures might be adapted to the skeletonizing of a grid 
network, they have no place in arterial network analyses. 

Only the more progressive water works maintain a program of field meas- 
urements of pipe friction loss coefficients. Usually time and fund limitations 
restrict measurement programs to the arterial network; only scattered sam- 
plings of grid network mains can be made. In the city of Philadelphia, for ex- 
ample, there are over 2,600 miles of mains equal to or smaller than 12-in 
dia, but only about 360 miles of mains larger than 12-in. With few exceptions, 
the arterial mains are larger than 12-in. The author refers to the head loss 
“which actually exists” as a point of reference in defining a “range of error.” 
It is seldom that a head loss coefficient is known to any resonable degree of 
exactitude for the grid mains with which the paper deals. In fact, one might 
ask how a flow test can be performed properly in the presence of active ser- 
vice connections. 

Let m represent the power to which the flow rate is raised to compute head 
loss, The author has assumed m = 1.85 throughout his paper. Arguments ad- 
vanced for the use of m = 2.00 for “old” pipes in arterial networks are not 
necessarily sound,1!2 With the little information normally available on grid 
mains it would seem quixotic to argue whether m = 1.85 is more suitable than 
2.00 for grid mains, and this discussion includes both powers in relations to 
be presented. For the low Reynolds numbers characterizing many grid mains, 
use of either 1.85 or 2.00 and assuming constancy of friction coefficient, re- 
gardless of Reynolds number magnitude and variation, are both unrealistic. 
Further, there are more local losses (fittings, valves, and specials) in grid 
mains than in arterial. Field flow coefficient valuesfor arterial mains include 
the influence of all local losses over the reach of measurement. Although local 
losses in grid mains may be neglected, asin the paper, the consequent exacti- 
tude of the results is even further diminished. 

While the paper implies that grid design is based on domestic demands, in 
the United States, fire flow requirements are often more severe, It is also 
implied that the local domestic demands, q, are known as a precise fraction 
of the imput to the local main, Q;, which is not likely. 4 

Retiring from the preceding reservations, the remainder of this discussion 
will be directed towards the development of the equations in the paper, As in 
the paper, a constant pipe coefficient (embodying friction coefficient and dia- 
meter) and equally spaced, equal-valued drawoffs will be assumed, . 


eG 
11 “Generalized Distribution Network Head Loss Characteristics,” by M. B. McPher- 
son, Proceedings, ASCE, Vol. 86, No. HY 1, January, 1960. P 
12 Discussion of “Relationships Between Pipe Resistance Formulas,” by M. B. Me- 
Pherson, Proceedings, ASCE, Vol. 86, No, HY 9, September, 1959, pp, 143-155. ; 
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The development of the equivalent of Eq. 9, has appeared in perhaps a 
dozen American hydraulics and fluid mechanics texts, but not in such general 
terms. One of the earlier presentations was by Dougherty. 3 Fair and Geyer 
have used the same principle in describing the hydraulic losses in filter under- 
drain manifolds. These two references use m =2.00. The author refers to the 
development of Eq. 9 as being for an “infinite number of connections” which is 
identical with the idealized uniformly decreasing flow (or velocity) assum- 
ed +»4* along with constant diameter and friction coefficient. When the flow 
into the service connections is large and the number of services are few (for 
example, n less than 10, q/Qt more than one-half), performance is more com- 
plex. It cannot be described properly by means of the simple pipe friction re- 
lations to which the paper is restricted because the hydraulic behavior of a 
manifold system15 is being approached. It is not clear why the series com- 
prising Eqs. 4 and 5 were introduced because Eqs. 1 and 6 are readily inte- 
grable, as shown in Eqs. 7 and 8. 

In the section on a “finite number of connections” (n), there are typograph- 
ical omissions and anerror in Eq. 12. Eq. 12 has been corrected and expressed 
in terms of K"' in Eq, 34. 


Referring to the definition sketch of Fig. 10(a) it will be noted that for an 
jdealized case the number of pipe lengths between cross-mains must be one 
plus the number of service locations (two services served at each location; 
there are thus six service locations in Fig. 10(a)). In Fig. 3 the author has 
‘either one too many connections or one too few pipe lengths. Eqs 10 through 
14 and much of the remainder of the paper are based upon this error. Letting 
the number of service locations be n' (that is a differentiation from n of the 
paper), the corrected, general expression for K is given by Eq. 35. 


: 1 n' “2.,| 
eee rae (ee eee (35) 


Table 9(a) lists values of K' for Eq. 35 using n' and m = 1,85, Eq. 9 was, 
therefore, used for an infinite n'. Curves which should replace those in Fig. 6 
are given in Fig. 10(b). It should be noted that for a realistic number of ser- 
‘vice locations, K' is practically constant, as opposed to the large variation 
indicated in Fig. 6 or Table 4. The most divergent case is for a q/Q of one, 
‘which would be for adead-end main or a main in which flow was supplied from 
both ends. Eq. 35 is not proper for this case because one too many lengths 
‘would be used—this is the only instance in which Fig. 3 and Eq. 34, would 


apply. 


Rel tal? saranumas Oe eG Fe eee 
~ 18 «ydraulics,” by R. L, Dougherty, McGraw-Hill Book Co., Inc., New York, N. Y., 


950, pp. 244-245. ; 
F 1 e Water Supply and Waste-Water Disposal,” by G. M, Fair and J. C. Geyer, John 


Wiley and Sons, New York, N. Y., 1956, p. 687. 
¥ 18 “Engineering Hydraulics,” Hunter Rouse, ed., John Wiley and Sons, New York, N. 


Y., 1950, p. 436. 
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Values from Eq. 35 for m =2.00 appear in Table 9(b) for comparison, Note 
that there is little difference in the values of K from those in Table 9(a) for 
an m of 1,85, For m = 2.00 Eq. 9 becomes Eq. 36. 


3 
K' = ste; (t-(1-&) ge Me (36) 


The most partinent section of the paper deals with the “ratio removed at the 
beginning of the line,” r. Eq. 19isnot affected by the error in the development 
of Eq. 34, and has been expressed in terms of K' in Eq. 37. 


Values of r for Eq. 37using K from Table 9(a) appear in Table 10(a) and from 
Table 9(b). While the author’s error leads to a large dependence of r on n, 
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/ 


Y 9 DISCUSSION 165 


ig. 7 and Table 8, the values of r in Tables 10(a) and 10(b) indicate an ap- 
roximate value of about 45% for most of the realistic combinations. The 
ithor says that “European practice has been to assume...45%...”, which ap- 
2ars reasonable. On the few occasions in which this type of problem has 
risen, the writers have used a 50-50 division, partly to allow for flow from 
ther direction but mostly because further refinement could not be justified. 
his practice, to the best knowledge of the writers, is common in water and 
as works practice; the writers are not aware of any utility where “a single 
jantity at the end of the line” is used. 

The assumptions leading to Eq. 35 are little better than those leading to 
q. 34. In Fig. 10(a), it may be seen that the assumption of equally spaced 
ads involves a series of premises: building lots must be equal in frontage, 
2rvices must be locatedinthe center of each lot, distance from center of side 
reet to property line of corner lot must equal spacing of services, or lot 
ontage. In the opinion of the writers these are negligible assumptions com- 
ured with the assumption that the design peak q, equally distributed among 


TABLE 9.—VARIATION OF K' FOR EQ. 35 VERSUS NUMBER OF LOCATIONS n!'. 


', No. of Ratio of Quantity Withdrawn (q) to Total Inlet Flow (Q 
toeations [00 [02 [ 04 | 06 [ 08 | 10 
(a) m= 1,85 


(b) m = 2,00 


4 Via Eq. 9. » Via Eq. 36. 
Le 


srvices in the bargain, is known to an exactness consistent with the numer- 
‘al values which issue from its use. 

If elegance is desired, a paper by Brown and Markel!® should be consulted. 
hey have mathematically treated the case of varying diameter, varying flow 
efficient, unequal service flows, unequal lengths between services, illustrated 
1r an m of both 1.85 and 2.00. Their development leads to a two-load equiva- 
snt, The only assumption for an m of 1.85 is “that the sum of a MacLaurin’s 
ries (of a specified form) canbe nearly approximated by the sum of the first 
ree terms of the series.” No assumptions are involved for an m of 2,00. 
irther, their development leads to apair of equivalent loads which give iden- 
cal losses regardless of direction of flow from one end or the other. Of the 
sveral discussions of their paper, the one by Malcolm S. McIlroy, inventor of 


16 “A Fluid Network Reduction Method for Accurate Representation of the Nonlinear 
aracteristics,” by Wilmer T, Brown and Robert A, Markel, AIEE Miscellaneous Pa- 
r 52-334, September, 1952 (Presented at Toledo Meeting, October, 1952), 
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the McIlroy Network Analyzer, is germane to this discussion, Because his dis- 
cussion is not generally available, it is reproduced here without deletion 01 
change: 


“The authors are to be congratulated for deriving the method pre- 
sented in the subject paper and for conceiving that it might exist. The 
method should often permit studies of pipeline networks to be carried 
out accurately on linear or non-linear electric calculating boards that 
would otherwise have inadequate capacity, and it should be generally 
helpful in studies conducted on paper. The opportunity of representing 
all the loads and pipeline sections between two adjacent network junc- 
tions by means of two pipeline elements, one intermediate load, and one 
terminal load is a very significant contribution to the art of pipeline- 
network analysis. 

“The suitability of the method for flow in either direction, subject to 
the stated limitations, and for exponents of any value encountered in 
fluid-flow problems, gives it wide usefulness. 

“In some studies of pipeline networks, exact values cannot be as- 
signed to loads at definite locations along the pipelines. In these in- 
stances, the assumption is often made that the total load delivered to 


TABLE 10.—VALUES OF r 


n', No. of Ratio of Quantity Withdrawn (q) to Total Inlet Flow (Qt 
octions [OZ [04 [| 08 | 08 | a 


(a) Via Eq, 37 using K' from Table 9 (a), m= 1.85 


(b) Via Eq, 38 using K' from Table 9 (b), m = 2.00 


users along the pipeline is uniformly distributed between network junc- 
tions and may be represented by two half loads at the ends of the sec- 
tion, This rough equivalence is false, asthe authors point out, but under 
many common operating conditions causes an error no greater than that 
inherent in the nature of known data. The error is not significant when- 
ever the flow delivered to other network elements at the downstream — 
terminal of the pipeline equals or exceeds the total load distributed 
along it. Wherever reasonable doubt exists regarding the ratio of dis- — 
tributed load to terminal load, or wherever definite values are known for ~ 
loads at specific locations along a pipeline, the authors’ method is a » 
powerful and accurate tool for reducing the number of quantities to be 


considered in an analysis. The procedure is a technical contribution of 
lasting value.” 


While some of MclIlroy’s statements would appear to refute several ared 
ments in the present discussion, it should be noted that the development b; 
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3rown and Markel was phrased in terms of gas network analysis. The magni- 
udes and variation in local gas loads are normally known to a much closer 
olerance than with water services (namely, the major share of the design load 
yften comprises the winter heating load which fluctuates rather predictably 
with temperature). MclIlroy’s discussion was, therefore, directed principally 
Owards gas applications. 

In summary, the common practice of splitting the service connection loads 
svenly between the two ends of a grid main is entirely consistent with the de- 
sree of accuracy to which other hydraulic factors can be evaluated. If grid 
main connections to an arterial network cannot be consolidated in an obvious 
und equitable manner, the procedure developed by Brown and Markel is pre- 
erred, but can be applied only between arterial main junctions. 

Some of the major design considerations for grid mains and networks in 
modern residential developments have been reviewed by the writers,! 

In Fig. 2, the term Q, is incorrect; perhaps Q, was intended. 


K. M. yao.!8_The arrangement of service connections as shown in Fig. 3 
S a bit different from the actual pictures. Usually, pipelines in a pipenet run 
‘rom one street intersection to another and it is unlikely to have service con- 
nection at both ends with possible exception of dead end lines, generally not 
sonsidered in pipenet analysis. Therefore, one service connection on the pipe- 
ine would mean a service connection at the midpoint of the line instead of at the 
ond of the line as assumed by the author. If this is the case, Eq. 10 will be: 


1.85 1.85 
A L 13852 S42 qd L nq 
pe” ae (0) 9+ Fe (@e- 8) teeter (@- 34) | haloes 


and Eq. 12 will be: 
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“17 “Distribution Design Considerations,” by M. B. McPherson and J. V. Radziul, 
al, AWWA, Vol. 51, April, 1959, pp. 489-502, 
18 Lecturer, Civ. Engrg., Univ. of Canterbury, Christchurch, New Zealand, 
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SEDIMENT PROBLEMS OF THE LOWER COLORADO RIVER@ 
Discussion by T. Blench 


T. BLENCH,!4 F. ASCE.—It is suspected that the authors, like the writer, 
acquired most of their knowledge of river behavior without the aid of formal 
instruction, for the simple reason that, even now, practically no colleges pro- 
vide any for engineers. Accordingly, it is hoped that they will not misunder- 
stand when the writer expresses the opinion that the knowledge required to 
forsee and arrange for the general trendof the events they have described ex- 
isted before Hoover and Parker Dams were built, but was not available to civil 
engineers generally; and he hopes they will agree with him that the absence 
from North American college curricula of some form of instruction on regime 
behavior of rivers and canals is as deplorable as would be the absence of soil 
mechanics. Actually, there is a little college activity in the matter already. 
The University of Alberta has given river andcanal instruction at postgraduate 
level since 1949, and is now (1960) introducing it in the Hydraulics, Soils, and 
Municipal Engineering options at undergraduate level. The University of Sas- 
katchewan has introduced a postgraduate course. The writer is aware of a few 
colleges in the United States that have given postgraduate instruction for some 
time. A Canadian of Indian extraction “drew attention to the value of teaching 
river and canal engineering, as is done in India, and... recommended... in 
the July issue of the Engineering Institute of Canada Journal” at the 1959 Con- 
gress of the International Association for Hydraulic Research. 

That article!5is relevant tothe present paper because it (1) gives six hypo- 
thetical typical situations in which civil engineers could make major errors (in- 
cluding lack of foresight), because they lack the knowledge of specialized col- 

eagues and geologists, (2) draws attention to official figures of the cost of 
sedimentationtroubles in the United States, (3) discusses the magnitude of er- 
rors that canbe made with the aid of modern construction equipment, (4) gives 
the opinion that “ignorance of river behavior is not universal; the trouble is 
that it is the prerogative of the civil engineers who build and plan water pro- 
jects,” (5) argues that the only effective method of removing this ignorance is 
: “offer civil engineering students the opportunity to study appropriately de- 

igned courses in river engineering, just as they now study soil mechanics, ° 
6) discusses the need for suitable textbooks and refers to the writer’s own,~*° 
(7) emphasizes that the geologic nature of river phenomena prevents any single 


@ April, 1960, by Whitney M. Borland and Carl R. Miller. 

14 prof,, Civ. Engrg., Univ. of Alberta, and Pres., T. Blench and Assoc. Ltd., Cons. 
Hydr, Engrs., Alberta, Canada, 

15 “River Engineering as a College Course for Civil Engineers,” by T. Blench, En- 
neering Institute of Canada Journal, July, 1959. 
_ 16 “Regime Behavior of Canals and Rivers,” by T. Blench, Butterworth’s Scientific 
Publications, London, England and Toronto, Canada, 1957. 
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person from acquiring expertness in terms of nothing but the experiences of 
his own short life, and (8) discusses the feasibility of college instruction. 

The paper necessarily condenses theory and the authors’ extensive field in- 
formation without which, of course, even the most satisfying theory cannot be 
applied satisfactorily. The writer, therefore, would appreciate further eluci- 
dation of the following points: 

1, In view of the facts that (a) Fig. 5 shows the bed material (median size 
about 0.3 mm) to follow closely the logarithmic Gaussian distribution of parti- 
cle size that sedimentary petrologists seem to have found general for sandy 
sediments, 17 and that the writer himself has found in river bed sands; 18 and 
(b) the bed material is following this distribution even after several years of 
retrogression has occurred, what is the factual or theoretical justification for 
believing differentiation willnow occur as implied in the statement “It was de- 
termined that bed stability wouldoccur when a6-in. armor layer of 4 to 10 mm 
material was obtained.” Also, assuming that this differentiation of coarse ma- 
terial is possible, in general, can a computation and description of transport 
mechanism be produced to show how 6 in. of it could differentiate out of the 
particular sand in the depth estimated as yet to be degraded? Does the All- 
American canal show such differentiation, and is its behavior relevant? 

2. How were channel widths found in Fig. 7 from Leopold and Maddock’s 
work’ when they did not propose any design formulas there, but instead re- 
corded relationships with very different constants found for very different rivers 
whose bed-material sizes were generally not observed? How were widths found 
by any of the methods when none of them concerned erodibility of banks? Why 
is the width for Q = 15,000 cfs and 0.3 mm material about 2,500 ft when the 
All-American Canal and afigure shown elsewherel9 for 20,000 cfs in the river 
at Yuma show that nature is satisfied with about 450 ft? 

3. Why was canalization adoptedas asolutionto the difficulties in degrading 
reaches when the effect of straightening a meandering channel is to reduce its 
regime slope radically and shorten its length, thereby causing it to degrade 
more? 

4, Is the canalized channel expected to meander? If so, what amount of 
stone protection is foreseen for the sides? 

5. Is the proposed section in Fig. 9 related to the channel dimensions ob- 
served by Leopold and Wolman!9 to have developed at Yuma since about 19407 
Their work seems to show that the channel is incised at discharges as high as 
20,000 cfs and had acquired, for that discharge, a breadth of about 450 ft from 
1940 to 1950, and a somewhat fluctuating depth that changed from 11 ft in 1940, 
to 13.5 ft in 1943, and back to about 10.5 ft in 1949, 

6. Considering that the section just mentioned is for an appreciable bed- 
load, obtained by dumping the total sediment of the All-American Canal inte 
the river (whose discharge is probably less than that of the canal), would the 
section needed for vanishing bed-load not be considerably narrower and deeper ? 

7. By what computations was “dominant discharge” decided? : 


17 “Size Frequency Distribution of Sediments and the Normal phi Curve,” by W. C. 
Krumbein, Journal of Sedimentary Petrology 8, 1938, pp. 84-90. “Gq 


18 “Normal Size Distribution Found in Samples of River-bed Sand,” by T. Blench, 


Civil Engineering, February, 1952. (Answers to discussion, February, 1953), “4 
“The Hydraulic Geometry of Stream Channels and Some Physiographic Impli : 

tions, by Luna B, Leopold and Thomas Maddock, Jr., Geological Survey, related 

Paper 252, 1953, Fig, 28, 
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8. Dothe authors attribute any significance to the following results from re- 
rime theory16 taking a bed factor 1.05 for 0.3 mm sand with vanishingly small 
wharge, and a side-factor of 0.20 for moderately cohesive sides: 


(a) Regime slope of straight canal of 15,000 cfs = 1/11,000. (So regime 
slope of meandering canal approximately = 1/5,500.) 

(b) Regime breadth for 20,000 cfs = 310 ft 

(c) Regime depth for 20,000 cfs = 16 ft 
Note that slope is very insensitive to discharge and side factor but roughly pro- 
jortional to bed-factor. Breadth varies as the square root of the bed-factor 
livided by the side-factor, and depth as the cube root of side-factor and in- 
versely as the 2/3 power of bed-factor, so b and c are not intended to be more 
han a guide. 

9. Why was the problem not solved by model originally? The reasons fa- 
yoring a model include (a) there would be no bed-load injection, (b) field data 
of every kind are available, or easy to collect, (c) models using natural bed- 
sand have been used satisfactorily for decades in India, and (d) regime theory 
srovides the means for scaling such models rapidly. (e) Friedkin29 demon- 
strated at the U.S. Waterways Experiments Station that the meandering of riv- 
ars, including degradation and aggradation, could be reproduced qualitatively 
using even artificial sediment. (f) Leopold has reproduced the details of braid- 
ing by using a model with suitably chosen sand.2! The conditions seem remark- 
ably favorablefor a quick accurate model solution, provided regime principles 
were used to eliminate unnecessary trial-and-error with scales. 

In the present state of river engineering, the writer feels that the publication 
of quantitative field information in accessible and discussable form is of high 
value, and he congratulates the authors on presenting so much information on 
a major project within the space available for a single paper. 


ee ae ee in fg i ae 
20 «4 Laboratory Study of the Meandering of Alluvial Rivers,” by J. F. Friedkin, U. 


3. Waterways Experiment Station, Vicksburg, Miss., 1945. 
' 21 «River Channel Patterns; Braided, Meandering and Straight,” by Luna B, Leopold 


and M, Gordon Wolman, U.S.G.S. Professional Paper 282-B, 1957. 
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TOLKMITT’S BACKWATER AND DROPDOWN CURVE TABLES* 
Discussion by Praxitelis A. Argyropoulos 


PRAXITELIS A. ARGYROPOULOS,21 F. ASCE.—The writer will comment 
on the formulas and tables by G. Tolkmitt for the computation in the actual 
practice of backwater and dropdown curves and to give practical and theore- 
tical supplementary information on this method. 

The Tolkmitt’s method is well known in many countries of Europe as well 
as in Greece. On the other hand, the Tolmitt’s formulas and tables for the uni- 
form prismatic channel of broad parabolic cross section form are not only 
found in the German texts.22 

It is very usual and important in Hydraulic Engineering to determine the 
stage at which, in a given artificial or natural channel, the flow takes place. 
On the other hand, it is well known that when a certain point of the channel is 
obstructed by a dam, sluice gates, weir, or other construction the raising of 
the latter will cause the water to set back up the stream. In this case, the 
water depth is increasing, and a new surface is formed. 

Usually, in the case of a backwater curve, two questions arise: how much 
the water is raising at agivendistance upstream from the point of obstruction 
and what is the distance for which the influence of the curve is felt. The ques- 
tions, in the case of dropdown curve, are analogous, 

Formulas, based on the assumption of channels possessing regularity in 

shape of cross section and having a uniform and constant bed slope, have been 
formulated. Also, in these formulas the roughness of the lining is assumed 
to remain constant, On the other hand, to facilitate the practical computations, 
tables were prepared. 
: Generally, in the case of a prismatic channel, the methods of computing 
backwater and dropdown curves can be divided into the following categories: 
a. Analytical methods; 
b. Graphical or semigraphical methods; and 
c. Methods by nomograph. 
The methods of Rthmann, Bresse Schaffernak, Tolkmitt, Bakmeteff, Du- 
i ‘ 23 24 2: 
 puit, Jansen, Chow, Von Seggern, Mononobe,“* Ming Lee,“” and Francesco 


| a2 May, 1960, by R. D. Goodrich, 
21 pr. of Civ. Engrg., Special Scientific Researcher, Athens, Greece. 
22 «Hydraulics,” by P. A. Argyropoulos, Ed. of the Technical Chamber of Greece, 


1951, pp. 354-361. 
23 “Integrating the Equation of Non-Uniform Flow,” by M. E. Von Seggern, Proceed- 


ings, ASCE, Vol. 75, No. IR 1, January, 1949. 
4 «Backwater and Dropdown Curves for Uniform Channels,” by Nagaho Mononobe, 


Transactions, ASCE, Vol. 103, 1938. : 
‘ 25 “Steady Gradually Varied Flow in Uniform Channels on Mild Slopes,” Univ. of Il- 


_linois, Urbana, Ill. 
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Ramponi,26 are noted. The method of Flavio Scano2” and Escoffier-Raytchine- 
Chatelain2® are graphical methods, while the Silber’s method29 is a method 
using a nomograph. 

On the other hand, in the case of an irregular natural stream with varying 
slopes of bed and different cross sections along its length, the solution of the 
problem of backwater and dropdown curves is more difficult and delicate. In 
this case, except for the tedious trial methods, there are the Krivoshein 
method, 30 the Grimm31 method, the Erza method, a simple method of the 
writer,32 the Leach diagram33 method and others. 

Also, the methods for computing backwater and dropdown curves for the 
artificial channels of trapezoidal, parabolic, and other cross-sectional forms 
can be divided in two categories whether the effect of the change in kinetic 
energy owing the accelerated or delated flow of water is considered or not: 
Therefore, the methods of Bresse, Struckel, Kozény, and others consider the 
effect of velocity-head v2/2¢ (in which, V is the mean velocity of flow and g 
the gravitational acceleration) whereas, the Tolkmitt’s method, the Rthmann’s, 
the Schaffernak’s method, and others, do not consider it. 

Tolkmitt, who was particularly interested innatural watercourses, assumed 
in his method a constant value of friction factor C throughout the whole 
range of depths, Tolkmitt chose as an “idealized” profile, a parabolic channel 
(Fig. 3). Let To be the top width and Ao the cross sectional area. 

The raising z the width Ty will be 


(17) 


and the wetted area: 


2 2 d+z 
A=5T (d+2) = 3 To Woe" (+2) ......... G8) 


On the other hand, the hydraulic radius: 


R= $F = OTe) 2 ds 2) vere see ty 09) 


26 “Determination approchée de la grandeur du remous dans les canaux rectangu- — 


laires ‘a faible pente,” Energie Electrice, December, 1951, 


27 “Procédé graphique rapide pout tracer la courbe de remous de 1’ecoulement per- 3 


menent dans les canaux prismatiques,” Energie Electrice, July, 1952. 


28 «Determination graphique de la ligne d’eau et calcul des remous,” La Houille © 


Blanche, No. 3, 1950, 


29 “Etude et trace’ des ecoulements permenents en canaux et rivi¢res,” Ed., Dunod, © 


Sat France, 1954, 
“Méthode naturelle de calcul des courbes de remous,” Rev. Universidad, Nati 
de Plata, December, 1951, : as 


: “Blackwater Slopes Above Dams,” Engineering-News Record, Vol, 100, 1928, 
“Computation of the Backwater and Dropdown Curves in Natural Streams,” Tech- 
nical Records, Technical Chamber of Greece, No. 401-402, 1957, 


33 “New Methods of the Solution of Backwater Problems,” Engineering-News Records 
Vol. 82, 1919, i 
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Parabola 


Tg=p,d 


FIG, 3 
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where, P= ~T=T y yz, is the wetted perimeter. 


In this case, the mean velocity 
V=CVRS =CVS (d+2) 8 .....7..05.. (20) 


in which, S is the fall of the free water surface per unit length. 
On the other hand, the discharge: 


Before the raising of the water (owing the obstruction), the discharge 


QM =AoVo=FT acyl as=- Fy To 43/2 c VS .. (22) 


But (from Fig. 4) 
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The general differential equation of the backwater curve is: 


dz [(a+ z)* . a4] 
—— =S5 gr ey ema PPLE MPR RAE si) 7 
axe (d+ m2 zt 


a 9 DISCUSSION 177 
Eq. 27 may take the form 


S, dx R isle nak 
=|oot + Gots TANS oc 
= (d+z)4 - a4 ee 
or 
d 1 1 bans 
Sed =). = — Sage dz (29) 
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The integration of Eq. 29 leads to Eq. 1. the procedure to find the formula 
for the dropdown curve, whose correct writing is: 


d : : D : 
Un,z) “55 E (2) e rete ( “Se @): (*2) cone 630) 


is analogous. 
In actual practice, the Tolkmitt’s method can also be applied in the case 


of artificial channels or natural streams having other than parabolic form. In 
this case, the cross section of the channel must be deformated in equivalent 
parabolic profile of the same to the given channel top width T, (Fig. 5) and 
cross sectional area A,. For instance, if the given channel is of a regular 


trapezoidal cross section (Fig. 5), 


Tot+b 2 
‘e} = *5 to = 3 To OG s casn oe at ieomewente eke (31) 
or 
3 Ao 
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On the other hand, 
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in which 
a3 bt 
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Example 1,—A rectangular channel, 98 ft wide and 4.6 ft in depth, has a 
slope of 1 in 2,000 ( = .0005). The discharge is 1,483 cfs. By constructing a 
weir, the original water level has been raised 2.6 ft. Estimate (by the Tolk- 
mitt’s method): (1) at which distance upstream from the point of the obstruc- 
tion the water is raised 7 ft, and (2) what is the total length of the backwater 
curve? 

The depth of the equivalent parabola will be, 


tes nee Bot, a 
d=5 ee | t = 4.6 = 6.9 ft 


d+h =6,9 + 2,6 = 9.5 ft and, d+z =6,9 + 1,0 =7,9 ft 


+ 
a 


From the Tolkmitt’s table, for ash = 1,38 f (: “ £(X) = 1,235 


when, (83 ‘| f(1,14) = 0.818 
So, 
6,9 
(hz) 9995 (1,235 - 0.818) = 5,755 ft 


On the other hand, as 


d+z\ _, /6,9+007\ _ " 
f ( q ) =f (Sr) = £(1,001) = - 0,507 


and, 


6,9 
L(h,z) = 9905 (1,285 + .507) = 24,040 ft 


Summarizing the present discussion, we note that, in Europe, the Tolkmitt’s 


method for computation of backwater and dropdown curves, is considered 


valuable and satisfactory. European engineers feel a Special confidence in © 


this method, By Goodrich’s paper, the Tolkmitt’s formulas and tables, slightly 
enlarged by the addition of twenty items in each table (thus in many cases 


obtaining intermediate values by interpolation will be avoided), will be avail- 
able to American engineers, 
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HOOD INLETS FOR CLOSED CONDUIT SPILLAYS# 


Discussion by Bernard L. Golding 


BERNARD L. GOLDING, !1 M. ASCE.— Mr. Blaisdell and the engineers of 
the Agricultural Research Service are to be congratulated for recognizing the 
potentialities of, and fully developing, the hood inlet culvert discovered by 
Karr and Clayton at Oregon State College, in 1954. 

In the writer’s opinion, the very simple hood inlet is the answer to the 
problem of making culverts on steep (supercritical) slopes flow full (prime). 
This could result in the saving of many dollars on the cost of culverts of the 
yet-to-be built portions of the Federal Interstate Highway System as well as 
other highways. It is the intention of the writer to discuss this paper as far as 
its application to the design of culverts under highway fills. 

In the determination of a culvert size under a highway fill, the headwater 
depth upstream from the culvert entrance and the exit velocity of the flow from 
the culvert are the two principal factors that must be considered. The head- 
water depth upstream from the culvert entrance must be limited so that ex- 
tensive flooding of upstream property does not occur. The exit velocity of 
flow must not be so great as to cause excessive scour of the stream bed and 

banks which could conceivably cause a washout of the highway embankment. 
Generally, up to the present time, engineers inthe writer’s office have felt 
obligated to assume that a culvert would not flow full (an entrance control or 
orifice condition existed) in computing the headwater depth caused by culverts 
on hydraulically steep (supercritical) slopes, even though the normal depth of 
flow in the culvert was “greater” than the height of the culvert, and the actual 
headwater depth was greater than 1.5 times the height of the culvert. This 
design assumption was based on the many observations and reports of others 
as well as our own field and laboratory observations, 
Making the assumption that the culvert would not flow full, resulted in the 
headwater depth being greater than that which would have occurred if the same 
pipe was to flow full, or for a particular allowable headwater depth, that is 

the usual design situation. The full flow condition would result in a smaller cul- 
vert size than the entrance control condition. An example of the possible re- 
duction in culvert size that would result from the use of the hood inlet to insure 
full pipe flow will be shown subsequently. 

The use of a smaller culvert flowing full will not necessarily mean that the 
exit velocity will be any greater than in the case of the larger pipe size with 
entrance control. Generally, the cross sectional areas occupied by the flow in 
each case will probably be the same. In the past, engineers in the writer’s 
office have made the reasonable assumption that culverts with entrance con- 
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trol will flow two-thirds full at the exit and have placed rip-rap, stilling 
basins, and so forth, at culvert outlets based on this assumption. 

In computing the headwater depth for culverts acting as orifices (entrance 
control condition), nomographs of the Bureau of Public Roads!2 have been used. 
These nomographs have been reproduced in two publications. 13,14 

In order to compute the head H over the invert of the culvert entrance for 
the full flow condition, the writer has simplified the Eqs. 4, 5, and 6 as follows: 


Be Hp tiD = Sodio 2 vleinis 1 eee aupitercthe' iy aih (19a) 
in which 

i= (I= He) pe + ORE se A insula (19b) 
or 


2.52(1-Ke) | 466.18 n21 Q2 


Ht = p 16/3 10 


p? 
In deriving Eqs. 19a to c, the following simplifications and substitutions were 


made: 


1. The term cos (sin71 S) has been completely dropped from Eq. 5b as 
being a neglible quantity. 
2. Sol has been substituted for Z, where So is the slope of the culvert in 


ft per foot and 1 is the length of the culvert (not including the length of the 


hood), 


3. Manning’s roughness coefficient “n” has been substituted for the Weis- 
bach friction factor “f,” 


In order to further simplify the computation of the head H over the culvert en- 
trance, a nomograph is given in Fig. 16 that enables an immediate solution of 
Eq. 19c. In the solution of Eq. 19a, the reasonable assumption can also be 
made that the hydraulic gradient pierces the plane of the outlet end of the con- 
duit at a distance above the invert equal to 0.75 times the diameter of the cul- 


vert (B = 0.75) that enables the complete solution of the head H over the invert 
of the culvert entrance. 


Generally, the standard state highway department headwall is equal to the © 


Headwall Type Anti-Vortex Wall shown in Fig. 8(e), so that no other special 
device to inhibit vortices seems necessary. However, the writer would appreci- 


ate more discussion by the author concerning the effectiveness of this type of 


anti-vortex device. 


When a conduit flows full, itis usually a good idea for the designer to com- _ 


pute the subatmospheric pressure at the entrance of the conduit to make sure 


that it is not less than the vapor pressure or cavation that will occur at the — 


12 “Highway Drainage Manual,” Bur. of Public Roads, Div, Two, Hagerstown, Md., 
July, 1953 (Preliminary). 


13 “Highway Engineering Handbook,” K, B, Wood, ed., McGraw-Hill Book Co, Inc., 
New York, N, Y., 1960, 


14 “Design of Small Dams,” Bur, of Reclam., Gov, Printing Off,, Washington, D. C. 
1960, 
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H,=Head in feet 
Ke= Entrance loss coef ficient 


D= Diameter of pipe in feet 
n= Manning's roughness coefficient 
L= Length of culvert in feet 
Q=Design discharge rate in cfs 
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pipe entrance. When the conduit entrance is of the re-entrant type, such as the 
hood inlet, the computation of the subatmospheric at the entrance should never 
be ignored. The actual computation of the pressure at the conduit entrance is 
quite simple once the head H atthe entrance is known, From Fig. 17, it can be 
seen that 


hyp + he + hy = hga + (H-D).......-.--. (20) 


v..2 
wi 
in which he= entrance loss in feet = K, ae ; hr = reduction in pressure head 
g—-¥ 2 


caused by the contraction at the entrance = Ky a in which K;= pressure 
g 


BD 


FIG, 17 


reduction coefficient; and hsa = subatmospheric pressure head at the crown of 
entrance, Eq. 20 may be rewritten as follows: 


V2 

- (Koi Ky # Ky) he, (H-D)\}...0.1.-+. 2a 

hii 7 
or 

Vp i 

BPP (Ret Kyi Ky) (=D) a “7 


For streamlined entrances, the pressure reduction coefficient K,is quite low 
but for the re-entrant type of entrance, the pressure reduction coefficient ze 


se 
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ypproach unity. Generally, a subatmospheric pressure of 25 ft at mean sea 
evel is considered allowable in hydraulic design. 

For purposes of comparsion between the sizes and costs of culverts ar- 
eived at under the two design assumptions, the typical field conditions are set 
orth as follows: 


Design flow 175 cfs 

Length 200 ft 

Slope of culvert 6% 

Pipe Type Reinforced concrete (n = 0.015) 


Maximum Allowable Headwater 8 ft 


Full Flow Condition.— (Hood Inlet, Ke = 0.8) Substituting in Eqs. 19a and c, 
1 42 in. pipe flowing full will cause a head H = 7.9 ft over the invert of the cul- 
vert entrance. The exit velocity from the culvert will be 18.2 fps, Substituting 
in Eq. 22, the subatmospheric pressure at the crown of the culvert entrance is 
9.0 ft that is well withinthe minimum subatmospheric of 25 ft allowed at mean 
sea level, if a pressure reduction coefficient Ky = 0.8 is used in design. 

Entrance Control Condition.— (Square edged entrance) Eq. 23 may be used 
n lieu of a nomograph for orifice flow square edged concrete pipe entrances 
under low heads (H > 1.0 D, H< 3.0 D) 


Substituting in Eq. 23, a 54 in. pipe will cause a head of 7.2 ft over the invert 
of the culvert entrance. If the depth of flow is assummed to be two-thirds of 
the pipe diameter at the outlet, the exit velocity of the flow will be 20.0 fps. 
In the northeastern states, the cost differential between a 42 in. reinforced 
concrete pipe and a 54 in. reinforced concrete pipe in place is about ten (10) 
dollars per ft, that would result ina savings of $2,000 in the preceding typical 
culvert installation. If the reasonable assumption of three such culverts in 
each mile of roadway is made, a cost saving of $6,000 per mile or $600,000 
for each 100 miles of roadway would result. 

The writer earnestly hopes that state highway departments will take 
cognizance of the advantages of the hood inlet as a culvert entrance, will 
sponsor further tests of it, and rapidly incorporate it into their standards. 
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ERRATA 


Journal of the Hydraulics Division 


Proceedings of the American Society of Civil Engineers 


May, 1959 
p. 31, 6 lines from bottom ‘Varwich” should be “ Varwick” 
p. 45, 12 lines from the bottom, the word “to” is omitted after “formula” 


p. 45, 11 lines from the bottom, the symbol for Reynolds number should be 
R, not IR 


p. 45, bottom line and page 64, IR should be R. 

p. 46, Table I in headings, R.10-3 should be R. 1079 

p. 50, 6 lines above the Table, whole line should be omitted. 
p. 66, “Chintu Lai” should be ‘Lai, Chintu” 
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